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Abstract 

Lanczos-Lovelock models of gravity represent a natural and elegant generalization 
of Einstein's theory of gravity to higher dimensions. They are characterized by the fact 
that the field equations only contain up to second derivatives of the metric even though 
the action functional can be a quadratic or higher degree polynomial in the curvature 
tensor. Because these models share several key properties of Einstein's theory they serve 
as a useful set of candidate models for testing the emergent paradigm for gravity. This 
review highlights several geometrical and thermodynamical aspects of Lanczos-Lovelock 
models which have attracted recent attention. 
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1 Introduction 

The principle of equivalence, along with the principle of general covariance, suggests that 
gravity has a geometrical interpretation in terms of the spacetime structure. These princi- 
ples are adequate to determine the kinematics of gravity (viz., how curvature of spacetime 
influences the dynamics of spacetime) in an elegant manner. Much of the beauty attributed 
to general theory of relativity arises from this geometrical interpretation and the kinematics 
alone. 

To complete the picture, one also needs to determine the dynamics, viz., how matter 
determines the spacetime structure. Based on the structure of other physical theories, it ap- 
pears reasonable to expect that this should arise through a second order differential equation 
relating the metric tensor to suitably defined source built from matter variables. Unfortu- 
nately, we do not have any equally elegant principle to determine these field equations of 
gravity. 
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It is therefore reasonable to study the most general system of field equations which are 
second order in the independent variables, which are the components of the metric tensor. 
Under very reasonable assumptions, we are then led to a class of models first discussed by 
Lanczos [1] and Lovelock [2]. These field equations follow from an action functional which 
is a polynomial in the curvature tensor and reduces — uniquely — to Einstein's general 
relativity when D — 4. In higher dimensions (D > 4) the Lanczos-Lovclock models differ 
from Einstein's theory and have a richer structure. 

The purpose of this review is to describe several features of the Lanczos-Lovelock models. 
Given the fact that they are more general than Einstein's theory and have been the subject 
of intensive study in the recent years, it is not possible to cover all the varied aspects of 
these models. In order to retain focus in the review, we will concentrate on the aspects of 
Lanczos-Lovelock models closely related to the emergent paradigm for gravity [3, 4, 5]. It 
turns out that several features related to the thermodynamics of horizons carry over from 
Einstein's general relativity to Lanczos-Lovelock models. Since the emergent paradigm of 
gravity uses the thermodynamic features of horizons as a key input, this allows a natural 
extension of this paradigm to all Lanczos-Lovelock models. The discovery that several 
features of the emergent paradigm transcends Einstein's theory of gravity strongly suggests 
that the relationship between gravitational dynamics and horizon thermodynamics is more 
general than was first recognized. In fact, this has been a strong motivation to study Lanczos- 
Lovelock models of gravity in the recent years. Such studies have also unravelled several 
peculiar geometrical features of Lanczos-Lovelock models which deserve attention on their 
own. We will describe these features as well, as we go along. 

Unless otherwise specified, we work with a D dimensional spacetime with signature 
(—,+,+,...). Latin alphabets run over spacetime indices, whereas Greek symbols 

li,v, . . . run over indices representing a (Z? — 1) dimensional hypersurface. This would be usu- 
ally a spacelike hypersurface, except in Sections 3.2, 3.3 where it is timelike. Finally, we use 
upper case Latin letters A,B, . . . for indices representing a (D — 2) dimensional spacelike hy- 
persurface. The convention for Ricmann tensor we follow is: i?"^^^ = '^'^[■^"^ d]b~^'^^'^ muF"^ d]b- 



Wc use the symbol '=' to denote the fact that the equation defines a variable or a notation. 

This review is structured in terms of five parts consisting of Sections 2 to 5. Section 2 
introduces Lanczos-Lovelock models and their major structural and geometrical features. 
Section 3 begins by introducing the concept of horizon entropy in theories more general than 
Einstein gravity and uses it to study several aspects of horizon entropy in these models. 
Sections 4 and 5 describe the connction between Lanczos-Lovelock model and the emergent 
gravity paradigm. 

2 Lanczos-Lovelock models: Structural Aspects 

As emphasized above, one distinguishing feature of Lanczos-Lovelock models is that, although 
the action functional of the theory could be an arbitrary higher order polynomial in the 
curvature tensor, it leads to equations of motion which do not contain derivatives of the 
metric tensor of order higher than two. Originally, this class of theories were arrived at from 
a somewhat indirect route by Lovelock [2]. He wrote down the field equations in the form 
Eab = {'^/'2)Tab where Tab is the energy momentum tensor of matter fields and demanded 
that Eab satisfies the following requirements: (i) Eab is symmetric in a and 6; (ii) it is a 
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function of metric and its first two derivatives; (iii) it is divergence-less. Having determined 
the most general form of Eab satisfying these criteria, he obtained a Lagrangian from which 
one can obtain these field equations. (Note that we have not demanded Eab to be linear in 
the second derivatives of the metric. If such a restriction is also added, one obtains Einstein's 
theory with a cosmological constant.) It is, however, possible to obtain Lanczos-Lovelock 
models in a more streamlined and transparent manner along different lines (see e.g.. Chapter 
15 of Ref. [6]) which is what wc will pursue here. 



2.1 Lagrangians from curvature and metric: general results 

Let us consider a general Lagrangian built out of Ricmann tensor and the metric, but not 
containing derivatives of the Riemann tensor, and ask the question: What is the condition 
on such a Lagrangian so as to ensure that equations of motion should contain no derivatives 
of the metric tensor higher than second order? The answer, as it turns out. uniquely picks 
out the Lanczos-Lovelock Lagrangians. 

Wc will now describe this result highlighting the main conceptual points, many of which 
are of interest in their own right. Wc start with the action functional in _D— dimensional 
spacctime: 

A = [ d^Xy/^L (1) 

Jv 

where L is a scalar built from the metric and the curvature tensor. We can think of the 
Lagrangian L as a functional of any one of {Rated, R"" bed' ^cd---) combination with either 
gab g^^^ The variation of the Lagrangian with respect to curvature and metric will bring 
in tensors defined by the derivative of Lagrangian with respect to different kinds of curvature 
components and the metric, like, for e.g., 

ibed — ( 1 . nab _ / 



pabed _ / 1 . pab _ / 1 (t^X 

~\dRabedJg^: ~\dgabjR^^^, 

It is obvious from the above definition that P"'' is symmetric and P"-^'^'^ has the algebraic 
properties of the curvature tensor: 

pabed pbaed pabdc. pabed pedab. pa[bed] q ^<j^ 

Curiously enough, the mere fact that L is a scalar imposes several further relations among 
these derivatives which play a crucial role in deriving the Lanczos-Lovelock field equations 
[7] . The most important ones among these relations can be expressed in terms of the tensor 

^ab ^ pa^Jkpb^^^ 

which plays the role analogous to Ricci tensor of Einstein's theory. It can be shown (see [7] ; 
Appendix A summarises these results with proof) that the following relations hold: To begin 
with, P'^'' and W'' arc closely related to each other: 

dgabJn,,,, '^'"K.dRatjk, 



2^%fc U7^ ; i.e., P"" = -27e'^'' (5) 
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It immediately follows that TZ°'^ is symmetric; this is by no means an obvious result and does 
not follow from the algebraic properties of curvature tensor. Using Eq. (5), one can obtain 
several other relations connecting the partial derivatives when {g°'^, Rated) or {g'^^ , R"'hcd) 
or are used as independent variables to describe the system. If we use the pair 

(g'^^jRabcd) as independent variables, we have 



f dL \ f dL 



= 27^,™ (6) 



On the other hand, if we use the pair (g"'', -R^^cd) ^'^ independent variables, one can show 
that 



dL 



dg 



ini 



7^^™ (7) 



and, more interestingly, if we use the pair (5°^,-Rc^) as independent variables, we get 

dL 



dg' 



= (8) 



Finally, it can be proved that the tensor V nP°'^'^^^ is symmetric in a, h. 

An intuitive way of understanding these results is as follows: To construct scalar poly- 
nomials of arbitrary degree in the curvature tensor, using only the curvature tensor Rabcd 
and metric g'^^ (with index placements as chosen) , we need two factors of the metric tensors 
to contract out the four indices of each curvature tensor in each term. So Eq. (6) is clearly 
true if L is a polynomial in Rabcd- Any arbitrary scalar function of curvature and metric 
can be represented as a (possibly infinite) series expansion in powers of the curvature tensor. 
Because Eq. (6) is linear in L, if it holds for a polynomial of arbitrary degree in curvature 
tensor, then it should hold for arbitrary scalar functions which possess a series expansion. 
This (intuitive) argument also gives the result in Eq. (8) which states that in scalar polyno- 
mials constructed from the metric cannot appear and all the contractions must be with 
Kronecker deltas! This is obvious because, if we use 5°^ to contract any two lower indices, 
that will leave two upper indices hanging loose which cannot be contracted out because we 
do not have gab available as an independent variable. Once we have Eq. (8) the other results 
can be obtained from it quite easily. Of course, the argument based on polynomials does not 
capture the full generality of our results, which hold for arbitrary Lagrangians that are not 
necessarily polynomials. 

Having obtained these useful identities, let us now proceed with the variation of the action 
treating the Lagrangian as a function of Rabcd and g°'^ with the index placements as shown. 
This leads to the result: 



/ d^xy/^Vj6v^ (9) 
Jv 



where 



Eab = , — ^ — T^ab — 2V™V"Pamnb 



1_ fdy/^g£ 

= [ Qgab J 

nab - \gabL - 2V™V"Parn„6 (10) 
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and 

SyJ = [2P"'J"'^(V,<5grf,) - 26gd^iVcP''"')] (11) 

In arriving at the second equality in Eq. (10) we have used Eq. (6). Note that each of the 
three terms in Eab in Eq. (10) is individually symmetric in a, b. 

In order to obtain sensible field equations from SA = 0, we need to ignore the boundary 
term involving dv^. From Eq. (11) it is clear that both the variation of the metric and its 
first derivative normal to the boundary contributes at the boundary. In other words, we do 
not get a sensible variational principle if we only assume Sgat = at the boundary. (This 
problem, of course, is well known even in Einstein's theory.) To obtain a well defined action 
principle, we have to add some boundary terms to our action such that the variations of the 
normal derivatives of the metric are cancelled out. We will say more about it in Sec. 2.6 
but will assume, for the moment, that the boundary term involving Sv^ can be ignored or 
cancelled. In that case, taking the total action to be Atot = A + Amatter and defining the 
energy momentum tensor for matter from the variation of Amatter with respect to the metric 
in the usual manner, our field equations will reduce to Eab = (l/2)Ta6 where Tab is the matter 
energy momentum tensor. 

2.2 Lanczos-Lovelock action functional 

Since P"'"^'' involves second derivatives of the metric, the term V'"V"Pamn6 in Eq. (10) 
will contain up to fourth derivatives of the metric. Hence, to satisfy our requirement that 
equations of motion should not have derivatives higher that second, we need to impose the 
condition: 

(It turns out that the seemingly more general condition, V'^V'^Pamnb = does not lead 
to anything different.) Due to symmetries of Pabcd, this makes P"-'"^'^ divergence- free in all 
indices. In this case, the equations of motion become 

Tlab IgabL = (^^^) ^ = iTab (13) 

In arriving at the first equality we have used Eq. (7). Clearly, the field equations do not 
involve more than second derivatives of the metric but are, in general, nonlinear functions 
of the second derivatives (the only exception being the Einstein- Hilbert Lagrangian). Note 
that the structure of Eab is very similar to that in Einstein's theory with TZab playing the 
role analogous to Ricci tensor. 

Our original task of finding an action functional, leading to equations of motion which 
contain no derivatives higher than the second order, now reduces to finding all scalar functions 
of curvature and metric such that Eq. (12) is satisfied. We shall now describe how one such 
class of Lagrangians can be constructed but will not bother to prove that our construction is 
actually unique — for which we refer the reader to original literature [1, 2]. 

For this purpose, let us consider a general class of Lagrangians of polynomial form in 
curvature: 

M 

L^J2 [(Riem)™] (14) 

m=0 
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where Lm contains m factors of the Riemann tensor, chosen with index placements i?"^. Our 
condition in Eq. (8) then imphes that the Lagrangian is independent of the metric tensor and 
aU the indices of the curvature tensor must be contracted out using Kronecker deltas. All 
we have to do is to find the correct manner of contracting the indices so as to satisfy Eq. (12). 
Ignoring a constant term corresponding to to = 0, the next simplest case corresponds to to = 1 
which should be a Lagrangian linear in i?"^. The only way of contracting the indices of 
using Kronecker deltas leads to 

Leh = SllRlt; Sll = SlSl - S^Sl (15) 

where S24 is the determinant tensor involving products of two Kronecker deltas and we are 
using an obvious notation as regards contraction of indices where the numeral n actually 
stands for an index a„ etc. This clearly leads to Einstein's theory with Leh R. 

The next order Lagrangian (for m = 2) will be quadratic in This implies that P'^'"^'^ 
(defined by Eq. (2)) will be linear in the curvature tensor. The corresponding Lagrangian can 
be found by obtaining the most general fourth rank tensor built from metric and curvature 
which is (i) linear in curvature and (ii) possesses symmetries of curvature tensor and (iii) 
divergence-free in all indices. The choice is unique and given by 

jjabcd j^abcd ^ac ^bd ^bc ^ad _|_ JZ^^g"^^ (-^^) 

(In four dimensions, this tensor is essentially the double-dual of Riemann.) Using this we can 
find the to ~ 2 Lagrangian to be 

L = l {g^ag'^g-'g"' - ^g^ag'^g^'g'^ + S^S^g'-'g^') R',,iR\,, = \ [R""'" Rabcd - iR^'Rab + R'] 

(17) 

This is the first non-trivial example of a Lagrangian leading to second order equations of 
motion from a Lagrangian which is nonlinear in the curvature tensor. The Lagrangian above 
is called Gauss-Bonnet Lagrangian, and a more explicit construction of it is given in Appendix 
B. The Gauss-Bonnet Lagrangian can also be written in the form: 

LgB = (52468^13^57 (18) 

where (S2468 again the determinant tensor involving four factors of Kronecker deltas. 

The examples of Einstein-Hilbert and Gauss-Bonnet actions suggest an obvious general- 
ization for higher orders that indeed leads to the Lanczos-Lovelock Lagrangian, which — in 
D— dimensions — is given by the sum: 

L = ^c„iL„i (19) 

rn 

t{D) „ _j_ 1_ i-ai6i...a„,fc„, pcirfi uc^dm (OfW 

~ i67r 2™ '^I'^i - '^'"''™ "i''! " ^ ' 

where the tensor appearing in the right hand side of Eq. (20) is a completely antisymmetric 
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determinant tensor defined as: 



cmi hi . 



dot 



" S} 


SI 


sL ' 




ci di . ..Cmdm 




Sj 







(21) 



for m > 0. The lowest order terms, to = 0, 1 correspond to cosmological constant and 
the Einstcin-Hilbcrt action respectively, as can be easily seen by expanding the alternating 
determinant. For m = 1, Li = (167r)~-'^i?, and the factor of IGtt in the definition of Lm 
essentially changes the right hand side of equations of motion from the conventional 87rTab 
to {l/2)Tab- (Getting an explicit expression for L^s remains a complicated task since it 
requires evaluating the determinant of a 2m x 2m matrix; fortunately, we shall never require 
explicit form of Lm's for the results which we wish to cover in this review!) 

The Cm's in Eq. (19) are arbitrary coupling constants which are a priori unconstrained. 
Some restrictions on Cm's might arise while considering specific solutions in Lanczos-Lovelock 
theory (such as black holes), or these might be related to certain parameters of a fundamental 
(quantum) theory when we consider low energy effective action that contain such Lanczos- 
Lovelock terms. This second situation arises in string theory where C2, for e.g., is related to 
the string scale, C2 (x l^. For most of our discussion in this review, we shall assume c„i's to 
be arbitrary. 

To complete our argument wc now need to prove that, for this Lagrangian. Eq. (12) is 
indeed satisfied. To do this, note that, being homogeneous in Rabcd, the L^s can be written 
as 



Lr. 



1 dLm 

m dRabcd 



Rabcd = (1/m) P(tf i^abcd 



(22) 



It is obvious that this tensor is same as the tensor Pabcd defined earlier, restricted to a 
particular value of m. The condition in Eq. (12) is therefore ensured if 

VaP(tf = (23) 
This is easily proved using the standard Bianchi identity for the curvature tensor: 



^a^cd (m) ^aOcd**..... ^ 



(771— l)factors 



(>cd**...;^[aR 







Rll 

ctors 

Rll 



(24) 

(25) 
(26) 



where in the second step we have used complete antisymmetry of determinant tensor, and 
the final equality follows from the Bianchi identity V[ai?^^j = 0. We have also introduced a 
convenient symbolic notation of using similar symbols (such as • etc.) to denote contracted 
indices. This helps us display relevant expressions (especially in later sections) in a form which 
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conveys same information in a much more transparent manner. Therefore, the Lanczos- 
Lovefock Lagrangians certainly would lead to desired equations of motion (not containing 
higher derivatives of the Riemann tensor). That they are unique in doing so requires some 
more work, but the demonstration for Gauss-Bonnet given in Appendix B should at least 
make it plausible. 

Just as in the above result, almost all of the unique features of the Lanczos-Lovclock 
models arise due to (i) the combinatorial properties of the determinant tensor and (ii) the 
fact that the Lm's are homogeneous polynomials of order m in the Riemann tensor. In 
particular, the Lanczos-Lovelock Lagrangians are finite order polynomials; i.e., although the 
sum in Eq. (19) is defined for all m € Z, only a finite number of Lm's are non-zero in a 
given spacetime dimension D. This follows essentially from complete antisymmetry of the 
determinant tensor and can be seen as follows. For even D, the determinant tensor vanishes 
if m > D/2, whereas for odd D, it vanishes if m > {D — l)/2. Therefore, for a given D, there 
is a critical value m = rric such that the action is non-zero only for Lm 's with m < mc, where 
rric = [{D — l)/2] with [x] denoting the smallest integer not smaller than x. The dimension 
D = 2m is called the critical dimension for a given L^, since L^*^^™ vanishes identically. 

The structure of the Lanczos-Lovelock action in D = 2m is interesting in its own right, 
and highlights certain topological significance of the Lanczos-Lovelock actions. This will be 
discussed in detail in Sec. 2.5. In particular, the critical dimension for Einstein-Hilbert action 
(m = 1) is D = 2 in which Gab identically vanishes; similarly the critical dimension for the 
Gauss-Bonnet Lagrangian (m = 2) is D = 4 and it does not contribute to the equations of 
motion in D = 4. We will say more about critical dimensions in Sec. 2.5. 

It is possible to write down an alternate form of the Lanczos-Lovelock action, which 
will be shown (in Sec. 2.7) to be related to the 'holographic' structure of these actions. 
This alternate form is a direct generalization of the so called FF (or Dirac) Lagrangian for 
Einstein gravity. We will sketch the basic structure symbolically, and will describe the details 
while discussing holographic structure of Lanczos-Lovelock actions in Sec. 2.5. Consider any 
Lagrangian which can be written as L = Q abcdR'^^'^'^ such that (i) WaQ'^^'^'^ = and (ii) Qabcd 
has all the algebraic symmetries of Riemann tensor. We know from the above discussion that 
Lanczos-Lovelock models belong to this class. We can then (symbolically) decompose L as 
follows 

L^QR - Q{dT + FF) 

- QFF + d (QT) + FVQ 

- QFF + c>(gF) (27) 

One immediately sees, in the final expression above, (a generalization of) the structure of the 
FF Lagrangian familiar from Einstein theory. The surface term in the above split is important 
in studying the thermodynamical properties of horizons in Lanczos-Lovelock models, as we 
shall show in a later section. 

Before concluding this section we will also comment on the Palatini form of the variational 
principle for the general Lagrangian. In Einstein theory, the conventional metric variation 
of the action yield equations of motion which are equivalent to those obtained upon treating 
metric and connection as independent quantities and varying g"'^ and F^^. independently. The 
latter variational principle (the so called Palatini variation) actually forces the connection to 
be metric compatible. For the general Lagrangian L[g'^^, Rabcd], it is easy to show that, when 
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the connection and metric are varied independently, we obtain 

6{L„,V^) = V, g'^^PklSr';^ + [T^at - l9abL]5g^' - V, (2P^|5«'V^) St% (28) 

where T*^^^ is an independent connection. If we ignore the total divergence term and set the 
coefficients of Sg"'^ and ST'^j^ independently to zero, then we get the conditions 

Tlab - \9abL = ^Tab; V, (2P;p5")'V^) = (29) 

It is difficult to make sense of the second equation in general and obtain the correct field 
equations. But if we restrict ourselves to Lagrangians which satisfy the condition V^P^.; = 0, 
then the vanishing of the coefficient of '^r'"'^-^ requires Vi (3"'-^— g) = 0. This in turn can 
be shown to lead to V;^— g = = Vigah, leading to the standard relation between metric 
and connection. The vanishing of the coefficient of Sg'^^ now leads to the correct equations 
of motion. Thus for the Lanczos-Lovclock Lagrangians, the Palatini form of the variation 
leads to the same field equations, just as in Einstein's theory. 

The converse is also true [8]. If we assume that V^-y/— g = = V igab (which makes 
connection and metric compatible), then the second condition in Eq. (29) leads to '^iPlj = 
0. That is, the compatibility of metric and connection (which is sometimes taken as a 
manifestation of the equivalence principle) for the generic Lagrangian of the form L[gab, Rabcd] 
picks out Lanczos-Lovelock Lagrangians uniquely. (Same conclusion was also reached in [9] . 
See also [10].) 

2.3 Explicit form of the equations of motion 

The equations of motion for the Lanczos-Lovclock model, given by Eq. (13), can be explicitly 
computed by first determining the form of P'^^'^d. fYom the Lagrangian in Eq. (20). This gives, 
for a generic Lanczos-Lovelock action L = J2m ^mLmi the equations: 



where 



b(m) 



Tia 
b 



E 



IGtT 2'" d,i---<^rn.dm aibi 



2 167r 2™ ici<il---Cmrfm aihl 



R 



ambr, 



(30) 



The equivalence of the two expressions in Eq. (30) can be established through the identity; 



H <:aibi...a„tb„ 
"j"cidi...c„,d,r 



lin ()^ ^aifci...a„6„ 
ZA« "ci"jdi...c„,d,„ 



^aibi 



R 



iaibi...ambm pCi^i 



rial 01 
"jcidi...c 



■R 



Cmdn 



(31) 
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which, in turn, is easily verified by noting that the alternating tensor in the right hand side 
of Eq. (31) can be written as a determinant: 



jcidi 



dot 



caibi...ar 
ci c/i ...c„ 



(32) 



The first term on the left-hand-side of Eq. (31) therefore comes from the multiplication of Sj 
with the lower right block of the above matrix. The remaining terms in the determinant can 
be grouped as 



"ck"3cidi---dk 



-ldfc---c„ 



ri raibi ■■■ a,„6,„ 
dk jcidi---CkCk + i- 



(33) 



where I < k < m. Since this whole determinant is multiplied by the product of curvature 
tensors, only the piece antisymmetric in the pair {cfc, dk} will be picked up, producing a factor 
of 2 for each pair {ck,dk}- Further, each of the m such pairs contribute the same amount 
due to the symmetries of the alternating tensor and the curvature tensor. This gives another 
factor of m, so that the contribution of the remaining terms in the above determinant becomes 
equal to 2m times the contribution of any particular term, say, the term corresponding to 
the pair {ci,di}. Noting the overall minus sign in (33), we obtain the second term on the 
left- hand-side of Eq. (31), thereby proving the desired equality. 

There are three more results which are easily obtained from the explicit form of the 
equations of motion: 

1. The field equations resulting from a m-th order Lanczos-Lovelock Lagrangian, which 
involves in Eq. (30), become vacuous when D = 2m, the critical dimension for a 
given Lanczos-Lovelock term. This is immediately obvious from the second form of 
Eq. (30), on recalling that the determinant tensor is completely anti-symmetric, and 
hence the number of values that the 2m -I- 1 indices of the tensor can take (which of 
course is the same as the number of spacetime dimensions Z?), must be greater than or 
equal to 2m -I- 1. For D < 2m the vanishes identically. 



2. Since Eab = Ti-ab — \gabL arises from the variation of a generally covariant scalar with 
respect to it follows that Vai5^ = which is the analogue of standard Bianclii 
identity. (The proof proceeds exactly as in the case of standard Bianchi identity; see 
e.g., page 255 of Ref. [6]). This in turn implies that 



kij 



1 



dbL 



(34) 



One consequence of this result is the following: Suppose we have a theory in which the 
field equations given by are the projection of the Lanczos-Lovelock field equations on 
to null vectors. That is, suppose the field equations of the theory are given by: 



2P, 



kij 



(35) 
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for all null vectors Ua with being a covariantly conserved symmetric tensor. This 
implies the equations 



TI + K51 (36) 



where A is a constant. To sec this, note that Eq. (35) implies that 2Pj^^-'^ R'^^-j^ ^ ^6 = 
— A(S^ where A(a:;) is some scalar. Rewriting this equation as 2£'^' + (L + X)5^ = T^, 
taking the divergence and using V = = aT^ we get db{L + A) = 0. Therefore, 
A = — _L+ constant. In other words, a theory with null projection of the equations of 
motion is equivalent to the original theory with a cosmological constant added to it. 

3. It can be shown that [11] for each Lanczos-Lovelock term of order to, one can construct 
a fourth rank tensor T?.^"^-*^ from P^^^ and i?a6cd' ^'^^^ all the algebraic symmetries of 
the Riemann tensor and having the property that, the trace of it's Bianchi derivative 
V[a7?.|^-|j^^j vanishes, yielding a divergence-free tensor which is the analogue of Einstein 
tensor for Lanczos-Lovelock models. Moreover, it can be argued [11] that Lanczos- 
Lovelock Lagrangians are the only ones which admit such a description in terms of a 
Bianchi derivative. It is also worth mentioning that, in [12], it has been proved, using 
the above mentioned generalized Riemann tensor for pure Lanczos-Lovelock actions, 
that for any pure Lanczos-Lovelock term of order to, the vacuum is trivial with respect 
to Ti^j^i ^ in all odd dimensions D ~ 2m + 1; i.e., it corresponds to 'R-^j^i '' = 0, which 
one might call as Lovelock-flat. This generalizes the well known fact in Einstein gravity 
(to = 1); that the vacuum solutions in Z? = 3 are Riemann flat. 

The additional characterizations of Lanczos-Lovelock models like the one mentioned 
above (as well as the one arising from validity of the Palatini variational principle) 
are not connected (at least in any obvious way) to the original requirement of quasi- 
linearity of equations of motion from which these models were first obtained. While it is 
possible that all such "uniqueness" results for Lanczos-Lovelock theories are eventually 
connected to (and derivable from) the topological significance of the Lanczos-Lovelock 
action in the critical dimension, we do not have at present in the literature a proof of 
the same. Such a proof might have important implications for several of the results we 
shall discuss in this review, which hinge on very special algebraic properties of certain 
tensors encountered in the Lanczos-Lovelock theory. 

4. It also follows from Eq. (13) that, for the mth Lanczos-Lovelock Lagrangian, the 
trace of the equations of motion is proportional to the Lagrangian: 



E ^ g'^'E^b = -fl ^^i^"" = -[{DI2) - TO]L(„) (37) 



This off-shell relation (that is, a relation valid even field equations are not imposed) is 
easy to prove from the fact that we need to introduce m factors of to proceed from 
R'^hcd to Red ^^'^ ^^^^ V~9 is homogeneous function of 5"'' of degree --D/2. Using this 



bed cd 

we get = —2E/{D — 2m) = —T/{D — 2m) which allows us to rewrite the field 

equations of the TO-th order Lanczos-Lovelock theory in the form: 

2^ =Tl- S^bj^^f (38) 
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This is similar in structure to Rl = 87r(T^ - i5^(l/2)r), to which it reduces to in D = 4 
for the m = 1 case. Note that this result does not hold for the theory obtained by 
adding together Lanczos-Lovelock terms of different order. 

2.4 Diffeomorphism invariance and the Noether current 

In any generally covariant theory, the infinitesimal coordinate transformations a;" — >■ + 
^"(x) lead to the conservation of a current (called the Noether current) that can be obtained 
as follows: The variation of the gravitational Lagrangian resulting from arbitrary variations 
of 5g°^^ can be expressed as form, 

5{L^) = ^ {EabSg'^' + VaSv'') . (39) 

where the (5u° term leads to the boundary contribution. When the variations in Sg""^ arise 
due to the diffeomorphism x'^ — > .t" + we have, S{Ly/—g) = —^—gVaiLS,""), with 5g'^^ = 
(V"C'' + VT). Substituting these in Eq. (39) and using the (generalized) Bianchi identity 
Va-E°^ = 0, we obtain the off-shell conservation law Vq J° = 0, for the current, 

ja ^ (2£;°^^b + L^'' + S^v") (40) 

where S^v"" represents the boundary term which arises for the specific variation of the metric 
in the form 5g°''^ ~ (V"^'' + V''^"). It is also convenient to introduce the antisymmetric tensor 
(called the Noether potential) J"'' by J° = Vt J"**. Using the known expression in Eq. (11) 
for in Eq. (40), it is possible to write an explicit expression for the current J" for any 
diffeomorphism invariant theory. For the general class of theories we are considering, the J*^^ 
and J" can be expressed in the form 

jab ^ 2P'''"'''^VcU - {VcP'''""^) (41) 

J" = -2Vh {P"'^'"' + P""'"') VcCd + 2P'''"'''Vb'^cU - 4^<iVfcVcP°'"='' (42) 

where Pabcd = {dL / dR'^'"^'^). The Lagrangians we are considering here depend, upto this 
point, on metric and Riemann tensor, but not on higher derivatives of the curvature tensor. 
However, the expressions in this section hold even in the latter case, the only difference 
being that the tensor Pabcd will now be defined with additional terms involving terms such 
as dL/d {S/R) etc. Algebraically, Pabcd simply gives the "equation of motion" with Rabcd 
treated as the dynamical variable. 

In the case of Lanczos-Lovelock models, V aP'^^'^'^ — and one can obtain from the 
Eq. (11) the expression for the boundary term 5v^ to be: 

<5«' = -2P/'^'<5«r\^. (43) 

where the Lie derivative of the connection (which is a generally covariant tensor) is defined 
through: 

C^V\, = VaVbC - R\bmr (44) 

This leads to the expressions for Noether potential and current given by 

jab ^ 2P'''""^\/c£.d; J" = 2P'''""^Vb'^c£.d ^ (27^°^Cfc + f^e^") (45) 
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In this case, we can also write the current as: 

J' = 27^'''6 + S^v' = 27^°''a - 2P^^ S^T\^ (46) 

This resuh shows that, near any event V where behaves hke an (approximate) KiUing 
vector and satisfies the conditions C^gab ~ ^ (a^b) = and C^V^^f^ = (which a true KiUing 
vector wiU satisfy everywhere), we get 

ja ^ 27^°''6■ (47) 

Incidentally, the expression for the Noether current can also be obtained from the Palatini 
variation in Eq. (28). Assuming that the metric is compatible with the connection and 
considering variations arising from diffeomorphisms for which S{Lm\/—g ) = —\f--g Vi (-^C): 
we obtain an expression for the conserved current as: 

V. {l^i^ + + 2.g"'P'\, ,5^r^-^„) = (48) 

where FjJ^ is now the metric compatible connection. This yields the same expression as in 
Eq. (46). 

2.5 Action functional in the critical dimension 

As we said before, the field equations resulting from a m-th order Lanczos-Lovelock La- 
grangian become trivial when D = 2m, known as the critical dimension for a given Lanczos- 
Lovelock term. The Lanczos-Lovelock action itself also vanishes identically for D < 2m. 
But for D = 2m, the action itself does not vanish. This raises the question: What is the 
structure of Lanczos-Lovelock action itself in the critical dimension? The vacuous nature of 
field equations requires that the variation of the action will be a pure surface term in the 
critical dimension. But it is not obvious whether the action itself is a pure surface term. We 
discuss these and related issues in this section, beginning by stressing why such issues are 
interesting and non-trivial. 

Let us begin with D = 2, which is the critical dimension for m — 1 Lanczos-Lovelock 
theory, the Lagrangian of which is simply R leading to Einstein-Hilbert action. It is trivial 
to show, using the language of differential forms, that the integrand of the Einstein-Hilbert 
action R^—g vn D = 2 can be written as an exterior derivative of a quantity which can 
be expressed in terms of the tetrads (see Eq. (50) below). But what we are interested in is 
whether or not Ry^—g can be expressed as djR^ for a doublet of functions (i?°, R^) built from 
the metric and its derivatives without using tetrads in a general coordinate system. That is, 
the line element is: 

ds^ = A{t, x)dt^ -f 2C{t, x)dxdt + B{t, x)dx^ , (49) 

(in which we have not imposed any gauge or coordinate conditions) we want to find an 
explicit expression for {R^ , R^) in terms of A, B, C and their derivatives. This is again trivial 
to accomplish in the conformally flat gauge with C = 0, A = £? in which Ry/—g — djd^ In A. 
This has the required form but it will not survive if we transform to an arbitrary coordinate 
system. 

Since we already know the answer in terms of differential forms, one might think that there 
is no room for controversy about obtaining an explicit expression for R^ . However, a paper 
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by Deser and Jackiw [13] seems to suggest that this is impossible, which led Kiriushcheva 
and Kuzmin [14] to conclude that there are fundamental differences between the vierbein 
formalism and the standard formalism and, because of those differences, the Einstein-Hilbert 
Lagrangian being a total derivative in the vierbein formalism does not imply that it is a 
total derivative in the standard formalism. Obviously, this issue deserves clarification: on 
one hand, several papers claim (at least implicitly) that the action itself is a total derivative 
for all Lanczos-Lovelock models in critical dimensions but without concrete proof; on the 
other hand, there are arguments, summarized in the articles by Kiriushcheva and Kuzmin 
[15, 14], that these claims are incorrect. 

These contradictory conclusions originate largely due to lack of an explicit expression for 

, as a function of the metric and its derivatives, satisfying Lm^/—g = djR^ . We will first 
summarise some recent results [16] establishing that R^ can indeed be found in terms of the 
metric and its derivatives and giving an explicit expression for R^ for the case of D = 2. 
This R^ is not unique since one can always add to it any set of functions (x) which satisfy 
djf^ = 0. There is, in fact, one natural i?* associated with each unit normalized vector field 
in D = 2 and any two of them differ by an with djf^ ~ 0. 

We will also summarize similar results for = 4. The Lanczos-Lovelock term, for 
which D = 4 is the critical dimension, is the m — 2 Gauss-Bonnet term. The question 
again is whether one can write Lqb V~5 = ^j^"*- expression for R^ was given in [17] 
based on an expression originally given in [18]. Unfortunately, this expression was obtained 
through an invalid identification of frame indices with spacetime indices and turns out to 
be incorrect. We will describe how correct expressions can be obtained and outline a proof 
that any Lanczos-Lovelock Lagrangian in the critical-dimension can be written as a total 
derivative of functions of the metric and its derivatives. 

Let us begin with D = 2. In reference [16], three distinct methods for obtaining an 
explicit expression for R^ in this case are given, of which we discuss the one based on Cartan 
formalism because it is easy to generalise to an arbitrary Lanczos-Lovelock term, Lm- For 
this purpose, we introduce an orthonormal basis w° = tufdx^ (satisfying gij ~ ujfuj^r]ab) as 
a basis for the cotangent space. The spin connection, defined by the first Cartan structural 
equation (see, for e.g.. Appendix A of [16]) must be antisymmetric: r^oi = — i^io- The 
Einstein-Hilbert two-form can be written (up to an overall sign) as: 



(50) 



where 8"'' {dVL''^ + VL"-^ A fi*'') is the curvature two-form. So, to find R^ such that 
Ry/—g = djR^ we only need to calculate explicitly the spin connection fioi and take its 
exterior derivative. This gives [16] the result: 



1 



-A— ffooa + (1 - A)— 3114 - 2(1 - A)5oi,i + .911,0 



500 



511 



-.9 



/T \\501 , \501 r,\ , 

(i — A) — .911,0 + A — 500.0 ~ ^A9oi,o + .9oo,i 
511 .900 



(51) 



which is an one-parameter family of solutions labelled by a constant A. This parameter can 
be thought of as describing the gauge freedom available in defining an orthonormal basis for 
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our cotangent space. The R^s for two different values of A are an example of non-uniqueness 
mentioned earlier. One can directly verify that = [R^Xi) — R^{X2)] satisfies djf^ = 0, as 
it should. This gauge freedom is entirely unphysical, and one can fix it without any extra 
condition on the spacctime or the metric. The gauge-fixed Cartan formalism is equivalent to 
the standard formalism for General Relativity. An explicit expression for R^ for the m ~ 2 
(the Gauss-Bonnet term) can also be found in Ref. [16]. 

We will now outline a general method for finding i^-^ 's for a generic Lanczos-Lovelock 
term, referring the reader to [16] for technical details. There exists a simple geometrical 
interpretation for the Lanczos-Lovelock Lagrangians being total derivatives of expressions 
involving tetrads. The Lagrangian L,„ is simply the Euler density of a 2m-dimensional 
manifold (without a boundary), as can be seen, for example, using the results in the section 
on index theorems in ref. [19]. In 2m dimensions, L„i is a topologically invariant scalar 
density because it is a closed form and therefore locally exact [20]. The Chern-Simons form 
Qm is defined to quantify this equality: evaluating it between two connections 51 and fi', 
one obtains Lm{^) — L,„(51') = dQm{^, ^')- In particular, if we restrict ourselves to a local 
coordinate patch, then we can introduce a fiat connection fJ' = and can therefore write 
Lm = dQm{^,0)- Hence, the existence of the Chern-Simons form is a valid proof that Lm 
is an exact differential. Moreover, this means that R^ can be thought of as providing the 
local coordinate description of the Chern-Simons form. Once it is clear that we can express 
Lm as an exact differential in the form language, we can proceed as in the two-dimensional 
case, and compute the set of (non- unique) functions R^ such that Lm\/~-g = djR^. One 
simply has to define an orthonormal basis uj"', calculate the spin connections fi"'', and use 
them in the expression already known in terms of differential forms. Therefore, LmyJ—g can 
be written as a total derivative of a set of functions of the metric and its derivatives in any 
critical dimension. 

We conclude this section by showing how R^ takes a very simple form in any spacetime 
possessing a Killing vector. This is a direct consequence of the results obtained in ref. [21] and 
the explicit proof goes as follows: Suppose the spacetime metric admits a timelikc Killing 
vector which we can choose, without loss of generality, to be: = (1,0, •• • ,0). Since 
Eab = 0, in the critical dimension, it follows from Eq. (10) in the case of Lanczos-Lovelock 
models with V aP"'^'"^ = that 

Lm = '^P^"^'' Rocde (52) 

We now note that 

Lm = 2P°'='^''i?0cd6 = P''"'^RacdbC = 2P°^'"'\/ cV d^b = -^d, (2 ^^pO^'^^ Vrf^) (53) 

where we used the facts that P"-^^'^ is divergence-free and has the same symmetries as the 
Riemann tensor. It follows that Lm\/—g = dcR"^ with i?^ = 2^—gP^'^'^^V d£,b- Using V d^b = 
r&odj we find that 

w = -2V^p/^ors, (54) 

In this case R^ is actually a component of the Noether potential corresponding to diffeomor- 
phism invariance, giving a direct physical interpretation to the local version of the Chern- 
Simons form in spacetimes which have symmetries. (For such spacetimes, the Noether current 
also gives the horizon entropy provided a Killing horizon exists; we shall discuss this further 
in Sec. 3.1.) 
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Obviously, similar results hold for spacetimes which have a spacelike Killing vector; if the 
Killing vector is taken to be, say, = in a coordinate system, we can obtain a similar 
expression starting from the identity El ~ 0. When the spacetime has more than one Killing 
vector, we obtain, in general, more than one choice of W which is again an example of the 
non-uniqueness mentioned earlier. 

2.6 Boundary terms for the action principle 

Part of the interest in the Lanczos-Lovelock models arises from the fact that they admit a 
standard initial value problem, because the equations of motion contain only up to second 
derivatives of the metric. At the level of the action, this means that one can hope to derive 
the equations of motion from a well defined variational principle by fixing the metric on 
given boundaries, without having to fix also the canonical momenta. For this to happen, 
one must be able to find suitable boundary terms which depend on the intrinsic geometry 
of the boundary, the variation of which can cancel the variation of normal derivatives of 
the metric on the boundary. Such a prescription is needed even for EH action, and the 
appropriate term is the well known York-Gibbons-Hawking term, given by the integral of the 
trace of the extrinsic curvature, K = h^'^Kfj^^, on the boundary [22]. The question arises as to 
whether one can construct a similar term for a generic Lanczos-Lovelock action. The answer 
is yes, and the boundary term itself is most easily obtained in the Hamiltonian formalism, 
which gives the canonical momenta the trace of which is proportional to the boundary term. 
For want of a better name, we shall call the Lanczos-Lovelock boundary terms as GYGH 
(Generalized York-Gibbons-Hawking). The first such analysis, to the best of our knowledge, 
seems to have been provided by Teitelboim and Zanelli [23]. 

Although we shall not discuss the details of the Hamiltonian formalism for Lanczos- 
Lovelock , it is relevant to mention one aspect of the same. For Lanczos-Lovelock action, 
the canonical momenta tt^j/ can indeed be found as a function of extrinsic and intrinsic 
curvature of the boundaries, but the corresponding expression cannot be inverted in closed 
form to find Kf^i, in terms of tt^i/ (i.e., one cannot express "velocities" in terms of "conjugate 
momentum"). More importantly, K^i, would generically be multiple valued functions of tt^i, 
for arbitrary Lanczos-Lovelock couplings (the Cm's). Hence the Hamiltonian evolution is 
not unique, and the generators of diffeomorphisms cannot be expressed in terms of canonical 
variables. Wc refer the interested reader to Rcf. [23] and Ref. [24] for further details on the 
Hamiltonian formulation of Lanczos-Lovelock models. Since wc will be working with the 
Lagrangian formalism all along in this review, the expression for boundary terms in terms of 
K^i, is sufficient for the specific topics we wish to cover in later parts of this review. Hence 
we shall focus on the same. 

In what follows, we give the boundary term for a timclike boundary. The corresponding 
expression for a spacelike boundary, which is the case dealt with in Ref. [23] , is similar except 
for a few (non-trivial) sign changes which one must be careful about. The timclike case turns 
out to be more relevant for study of quantum black holes since one usually describes the 
near-horizon physics by replacing the horizon with a timelike surface - the stretched horizon 
- in the spirit of the membrane paradigm. In fact, the boundary terms given in this section 
will be used later in Sees. 3.2, 3.3 where we will obtain some important conceptual insights 
into horizon entropy essentially by studying the near-horizon structure of these boundary 
terms. 
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For the Lanczos-Lovelock models, the final form of the relevant boundary term can be 
written as [25] 



r(m) 



dM 



2m 



1 r4l42i3-- 
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The above form can be further simplified, and the complete action can be written as 



M dM 



where 



(m)fj. om+l / y * 



(s)/i 
s=0 

— s-[l'fl---l'2777-l] 'TpA'lA12 . . . '7?A'2s-lAi2s ffp,2s + l . . . I('M2777- 1 /CyA 

'^(s)m - ''[/i/ii-. .,,2777-11 'Vi^2 ''^i^2»-11^2» ■'^I^2» + l -"■l'2777-l ^ ' 

and the coefficients are given by 

^m — s 

^ s!(2m-2s-l)!! ^^^^ 

As the notation already anticipates, the object ijyj-)^ defined above can be shown to be the 
surface stress tensor in Lanczos-Lovelock theory; that is, 

,(m) _ 2 (5^sur(m) . 

'^'^ - ^\ ^^^^ 

In the above expressions, TZp^pav denotes the projection of full spacetime curvature tensor on 
the {D — 1) hypersurface dM. It is more convenient to express the stress tensor in terms 
of the intrinsic curvature R^pau of the hypersurface (i.e., the one defined using the induced 
metric h^^j), through the Gauss-Codazzi relation 

Rfipav ~ T^^pau I^p,(7^pv ^ i^pvi^pa (60) 

After some simple combinatorial manipulations, the surface stress energy tensor can be rewrit- 
ten as 

I m— 1 

^^^{m)p - 2^ ^{s)p 

s=0 

-V _ r[!'I^l---l'2777-l] r,pip2 . . . 0^23-1^23 T^p2s + 1 . . . I<'M277i-l ( a-\\ 

^(s)p - °[ppi...p2^_i]^yiU2 ^1^2,-1^2. ^1^2s + l -'^V2^-l l^-^i 



18 



with 



m — 1 




(62) 



q=s 



One can check that for m = 1 and m ~ 2, the above expression reduces to that of the surface 
stress tensor in Einstein [26] and Gauss-Bonnet [27] theories respectively. 

Ahhough we have quoted the structure of the surface terms, we have not formaUy shown 
that these terms indeed make the variational principle well defined. In fact, this is a remark- 
able feature of the Lanczos-Lovelock actions and the reason why things work out again lies 
in the topological origin of the Lanczos-Lovelock actions. Wc refer the reader to [19, 28] for 
a detailed description of this point. It is however worth giving a sketch of the main argument 
here. 

As we mentioned in the previous section, in the critical dimension D ~ 2m, the Lanczos- 
Lovelock action Lm"^'^ gives Euler character of a manifold M without a boundary; in presence 
of a boundary dM, one needs to add a specific boundary term, say — (3m™\ to obtain 
the correct Euler characteristic. Since Euler characteristic can not change under smooth 
variations of a metric with a fixed boundary, the variation of im™'' — JgM '3™'"^ with 
respect to normal derivatives of the metric on the boundary must vanish; in other words, the 
corresponding variations in im"'^ and Qm™'' must cancel. This cancellation is independent 
of spacetime dimensions, and hence the above result is valid even away from the critical 
dimension. So the boundary term we need is essentially the term which — when added to 
the bulk term — will make the sum equal to the Euler characteristic in the critical dimension. 
From this one can show that the form of the surface term is given by the Chern-Simons form 
defined on dAi. (The above points are easily understood in Einstein gravity, where the 
standard York- Gibbons-Hawking term is precisely what is needed to make the full action 
Euler characteristic in _D = 2, which is the critical dimension for Einstcin-Hilbcrt term 
corresponding to m = 1.) This leads to the boundary terms quoted above. 

Wc stress a few important characteristics of the above surface term, some of which will 
be relevant from the point of view of topics to be discussed in Sec. 3.3. 

1. The surface term for a given Lanczos-Lovelock term of order m has an apparent pole 
at the critical dimension D = 2m. There is, however, no real pole since f^s itself has 
a zero at D = 2m which precisely cancels the pole. This will also become evident when 
we give explicit form for t'^^^-^^ near a static horizon in Sec. 3.3. In fact, what happens 
at D = 2m (which is the critical dimension) is straightforward to understand from the 
fact that, in Z? = 2m the corresponding Lanczos-Lovelock term is the Euler density 
of the manifold, and the above surface term must be included if the manifold has a 
boundary. 

2. For a given m, it is easy to see that the boundary term is a polynomial containing odd 
powers of K^i, of degree 1 to (2m — 1). For e.g., for the Gauss-Bonnet term (m = 2), 
the boundary term will be a linear combination of terms linear and cubic in K^^. 

3. The notation for the boundary tensor anticipates the fact that the jump in t^^i, across 
a shell, say, can be related to stress-energy tensor of the material comprising the shell, 
to yield the (generalization of) Israel junction conditions for Lanczos-Lovelock models. 
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This fact is useful in the study of collapse and black hole formation in Lanczos-Lovelock 
models. Further, the same junction conditions can be used to set-up the membrane 
paradigm for black hole solutions in Lanczos-Lovelock theory, in the sense that one 
can replace the details of region hidden by the horizon by endowing the stretched 
horizon with the above stress tensor. We expect the membrane paradigm to work in 
this case just as it does in Einstein theory, due to the initial value problem being well 
defined in Lanczos-Lovelock models. We will discuss this further in section 3.3. 

2.7 Holographic structure of Lanczos-Lovelock action 

In a classical theory, action functional is just a device to obtain equations of motion and 
to encode the symmetries of the system. The actual value of the action functional is of no 
importance. In contrast, action functional plays a more crucial role in quantum theory: in 
the scmiclassical limit, it gives the phase of the wave functional and in the full description it 
appears in the path integral. So the off-shell structure of the action functional is important 
in quantum theory. In Lanczos-Lovelock models we have already seen that the action has a 
peculiar structure and requires a special treatment if it has to lead to a well-defined variational 
principle. By and large, this requires careful handling of the surface term in the action 
functional. 

The existence of the surface term brings in another aspect in the Lanczos-Lovelock models 
viz. that the bulk and surface terms of the action are closely related to each other and 
duplicates the information to certain extent. This is reminiscent of the notion of holography, 
which represents — in the most general sense — a correspondence between bulk and boundary 
degrees of freedom. For lack of better terminology we will call the correspondence between 
bulk and boundary terms in the action (rather than in the degrees of freedom) as " holography 
of the action functional" . (In this review, we shall use the term 'holography' as explained 
above and it should not be confused with the usage of this term in e.g., string theory.) In 
this section, we shall explore this feature and show that Lanczos-Lovelock actions possess a 
holographic structure which eventually gets reflected in the horizon thermodynamics in these 
models. 

The basic idea we wish to build upon was already expressed, symbolically, in the decompo- 
sition of the Lagrangian in Eq. (27). We will first describe a toy example of higher derivative 
action in ordinary mechanics to demonstrate the fact that such an action must have a very 
special structure if it has to yield equations of motion which are still second order in time 
derivatives [29] . It is then not a surprise that the generalization to gravitational actions picks 
out the Lanczos-Lovelock action uniquely as the one exhibiting holography. 

2.7.1 Higher derivative actions yielding second order equations of motion 

Consider a dynamical variable q(t) in classical mechanics described by a Lagrangian Lq(q, q). 
Varying the action obtained from integrating this Lagrangian in the interval (^1,^2) while 
keeping q fixed at the endpoints, we obtain the Euler Lagrange equations for the system 
(dLq/dq) = dp/dt, where p{q,q) = {dLq/dq). (The subscript q on Lq denotes the variable 
that is kept fixed at the end points.) The Lagrangian depends only up to first derivatives of 
the dynamical variable and the equations of motion are — in general — second degree in the 
time derivative. 
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Interestingly, there exists a wide class of Lagrangians L{q, q, q) which depend on q as well 
but still lead to equations of motion which are only second order in time. We will now analyse 
this class which will generalize to the holographic actions in field theory. 

To do this, let consider the following question: We want to modify the Lagrangian Lq 
such that the same equations of motion are obtained when — instead of fixing q at the end 
points — we fix some other (given) function C{q, q) at the end points. This is easily done by 
modifying the Lagrangian by adding a term —df{q,q)/dt which depends on q as well. (The 
minus sign is for future convenience.) The new Lagrangian is: 



Lc{q,q,q) = Lq{q,q) - 



dJiq, q) 
dt 



(63) 



We want Lc to lead to the same equations of motion as ig, when some given function C{q^ q) 
is held fixed at the end points. We assume Lq and C are given and that we need to find /. 
The standard variation now gives 
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At the end-points, we invert the relation C = C{q,q) to determine q = q{q,C) and express 
p{q,q) in terms of {q,C) obtaining the function p = p[q,C). We then treat / as a function 
of q and C, so that the variation of the action becomes: 
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dL 
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(65) 



since SC = at the end points by assumption. To obtain the same Eulcr-Lagrange equations, 
the second term in Eq. (65) should vanish for all Sq. This fixes the form of / to be: 



f{q.C)= p{q,C)dq + F{G) 



(66) 



where the integration is with constant C and F is an arbitrary function. 

Thus, given a Lagrangian Lq{q,q) which leads to certain equations of motion when q is 
held fixed, one can construct a family of Lagrangians Lc{q,q,q) that will lead to the same 
equations of motion when an arbitrary function C{q, q) is kept fixed at the end points. This 
family is remarkable in the sense that though Lc will now also involve q, the equations of 
motion are still of second order in q because of the structure of the Lagrangian. (The results 
obtained above can also be interpreted in terms of canonical transformations etc. which we 
purposely avoid since we want to stay within the Lagrangian framework). So, even though 
a general Lagrangian which depends on q will lead to equations of higher order, there are 
special Lagrangians which will not. The analysis extends directly to a multicomponent field 
g^(x°) in a spacetime with coordinates cc" where A collectively denotes the tensor indices 
[29]. 

When one considers the action just as a tool to obtain the field equations, the above 
procedure can be used with any C . But once the dynamical variables in the theory have been 
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identified, there are two natural boundary conditions which can be imposed on the system. 
The first one holds q fixed at the boundary and the second one keeps the canonical momenta 
p fixed at the boundary. In the second case, C = p and Eq.(66) gives f — qp — q{dL/dq) 
leading to the Lagrangian: 

This Lagrangian Lp will lead to the same equations of motion when p{q, q) ~ (dL/dq) is held 
fixed at the end points. In field theory, with multicomponcnt field QAix"") with canonical 
momenta Tr\ = [dL/d{diq^)], when we fix = tt^ at the boundary, the new Lagrangian 
we need to use is: 

= -^bulk + Lsur (68) 

We will now summarise [29] the implications of the above results: 

1. These results have a natural interpretation in quantum theory: A path integral de- 
fined with Lp will give the transition amplitude in momentum space G{p2,t2',Pi,ti), 
just as a path integral with Lq leads to the transition amplitude in coordinate space 
K {q2,t2;qi,ti). 

2. When the momenta are held fixed, the surface term has, what we shall call, the "^(gp)" 
structure. It is obvious from this relation that the surface and bulk terms in the 
action are closely related and the form of the surface term will put some constraints on 

ibulk = Lq. 

3. It turns out that the Einstein- Hilbert action (and more generally the Lanczos-Lovelock 
action) has precisely this form, which is evident when the action is written in the FT 
form. This is the key to the holography in Lanczos-Lovelock action Junctionals which 
we now describe. 

2.7.2 Holography of gravitational action functionals 

Let us first discuss the Einstein-Hilbert action, rewriting the well known FF form of the 
action in a manner which generalizes to the Lanczos-Lovelock actions. The Einstein-Hilbert 
action can be written using a tensor Qabcd = (1/2) (gacgbd - Qadgbc) as: 

Qahcd r> 
nabcd 

25c [^/^gQa^'^nd] + 2^gQ,'"--^Tl^Tl 

isur + ibulk (69) 

with 

Lbulk = 2^gQj'''^Tl^Tl ; L.^r = 2d. [V~gQa''nd] ^ dc [^gV^] (70) 

where the last equality defines the D-component non-tensor V. As is well known, one can 
obtain Einstein's equations varying only keeping gab fixed at the boundary. 



Lp{d'^qA,dqA,qA) = Lq{qA,dqA) ~ 
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The non-trivial aspect relevant for holography is a specific relation [30] between Lbuik and 
Lsur allowing Lgur to be determined completely by ibuik- Using gab and dcgab as independent 
variables in ibuik one can prove that: 



1 



[{D/2)-l 



gab 



d{digab) 



(71) 



"d(ijp)" 



The connection with the discussion in the last subsection is obvious: the "d{qpy' structure 
of Lsur suggests that Leh is obtained from Lbuik by a transformation from coordinate space 
to momentum space. The action based on Lb^ik is the co-ordinate space action and leads 
to a well-defined variational principle if the metric is fixed at the boundary. On the other 
hand, the usual Einstein-Hilbert action {R^—g = Lbuik + Lsur) is actually a momentum- 
space action and leads to well-defined variational principle if the momenta are fixed at the 
boundary. This is very transparent if we use f^^ = 
the momenta being 



gg as the dynamical variables with 



\{5]ria + 5iT%)] 



(72) 



if we fix the momenta A'jj, on the 



The equations of motion will follow from SAeh 

boundary. In fact, if we decide not to add any surface term to Einstein-Hilbert action, then 
we can still obtain [5] the field equations if we demand 



SAeh = - 



av 



Sxy/hn,g^''5N\ 



(73) 



instead of the usual 5Aeh = 0. 

Similar results continue to hold for Lanczos-Lovelock actions, which we will now describe. 
Taking a cue from Einstein's theory, consider a Lagrangian of the form: 



Ik 



(74) 



where Q"'"^'^ is a tensor having the symmetries of curvature tensor. By straight forward 
algebra, one can prove [29] the following identity: 

-gTt.V.Qj'^'' 
(75) 



-gL = 29, [V~9Qa'"''nd] + '2^9Qa'''''^aAc 



The general covariance of L requires the condition VcQa = 0- Then we get 



(76) 



We have, therefore, once again arrived at the Lanczos-Lovelock models in a very different 
context, both mathematically and conceptually. The final step of establishing holography of 
such actions requires some additional care particularly since, unlike in the case of Einstein 
theory, Lbuik will in general depend on second derivatives of the metric. For Lanczos-Lovelock 
models, it is possible to show that [29] the following generalization of holography exists: 



[{D/2) ~ m]L, 
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bulk 



(77) 
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where the Eulerian derivative 5 is defined as 

<5i^[0,a,(/),...] ^ dK[(i),d,4>,- 

5(j> d(j) 

Note that the second term on right hand side of Eq. (77) vanishes for m = 1 Lanczos- 
Lovelock term (which is the Einstein theory). The proof of above relation, as well as some 
other forms of holography for the Lanczos-Lovelock class of actions (which arise by choosing 
the Christoffel symbols in place of first derivatives of the metric as independent variables), 
can be found in Ref. [29] . 



dK[cj>,d,cb,. 



d{da<P) 



(78) 



2.8 Black hole and cosmological solutions of Lanczos-Lovelock mod- 
els 

There is a vast literature on known solutions for gravitational theories based on the Lanczos- 
Lovelock Lagrangians, and even a brief discussion of these would form a topic for a new 
review. We shall therefore use this section just to redirect the interested reader to relevant 
literature and reviews where more details about specific solutions, and more importantly, 
additional references, can be found. 

As far as black hole solutions of the theory are considered, one of the very first solution to 
have appeared (to the best of our knowledge), is the so called Boulware-Descr solution [31] 
which is essentially a spherically symmetric black hole solution of 5 dimensional Einstein- 
Gauss-Bonnet (EGB) action. The metric for this solution is given by 



ds2 = _/(r)dt2 + J-^dr^ + rMrj^ 



16am AaK 
1±A 1 + ^ + — 



(79) 



where dVl\ is the metric on a unit 3-sphcrc, A is the cosmological constant, and a is the 
Gauss-Bonnet coupling constant normalised so that the total lagrangian takes the form: 
L (X —2A + R + a (i?^ — iRatR'^'' + RabcdR'^^'''^) ■ The parameter m is an integration constant 
which is proportional to the "mass" of the solution; its exact form with all the constants, can 
be found, for example, in [31]. The ± is a consequence of the fact that field equations are 
quadratic in f{r). There exists an extensive literature on properties of this solution, as well 
as its modification upon addition of the electric charge Q. 

There also exists considerable literature on study of Birkhoff 's theorem in EGB theory, 
the most relevant of which arc listed in [32]. Although the analysis leading to the Boulware- 
Deser solution presented in [31] was motivated by appearance of EGB term in low energy 
expansion of super-symmetric string theory, it was followed by extensive work on finding 
such solutions for generic Lanczos-Lovelock models treated as classical models of higher 
dimensional gravity in their own right. A nice review of such solutions, and their various 
asymptotic properties (which arc of some interest also for the AdS-CFT correspondence), 
can be found in Refs. [33] and [34] (see also [35] for work generalizing known solutions in 
Einstein theory to Lanczos-Lovelock ). A systematic study of cosmological solutions in the 
Lanczos-Lovelock theory can be found in [36], which deals in particular with maximally 
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symmetric space-times, Robertson- Walker universes, and product manifolds of symmetric 
subspaces. 

3 Lanczos-Lovelock models: Horizon thermodynamics 

One fascinating aspect of general relativity, which is still not completely understood concep- 
tually, is the association of thermodynamic variables with horizons. Originally, these results 
were obtained in the context of black hole physics [37, 38] but it is now clear that all null 
surfaces possess these attributes in an observer-dependent manner [39] . It is important to un- 
derstand whether this is a special feature of Einstein's theory or whether all Lanczos-Lovelock 
models share this intriguing connection. It turns out, as we shall describe in this and ensuing 
sections, that the latter is the case. 

As we said before, one key utility of the Lanczos-Lovelock models is in their power to 
discriminate between widely differing conceptual explanations for known results in Einstein 
gravity. So the study of horizon thermodynamics in Lanczos-Lovelock models will allow us to 
distinguish those aspects which transcend Einstein's theory and will lead to a more promising 
and general route to understanding these aspects. In this spirit, there has been an extensive 
study of thermodynamics of black hole solutions in particular and horizons in general in the 
Lanczos-Lovelock theory which we now turn to. Since it is impossible to do justice to all the 
literature in this review, we will be selective in our discussion. A good introduction to this 
topic can be found in the (possibly the very first) paper on the topic by Myers and Simon 
[40]. 

3.1 Horizon entropy in Lanczos-Lovelock theory 

The results from the study of quantum field theory in a given background spacetime does 
not depend on the action governing the gravitational dynamics of the background. In par- 
ticular, results [39] such as the Davies-Unruh effect etc continue to hold in a given metric, 
independent of the field equations of the theory. In fact, the notion of horizon temperature 
only involves the near horizon geometric properties of the spacetime metric; therefore, the 
notion of temperature for a horizon remains the same independent of the theory of gravity. 

The same is not true for the horizon entropy, though. Algebraically, this is easy to under- 
stand in the context of a simple black hole solution in Lanczos-Lovelock theory characterized 
by a single parameter corresponding to mass M of the black hole (identified in a suitable 
way, for example, by looking at the asymptotics of the solution). Then, the entropy can be 
identified by integrating dM/T{AI), for which one would need an explicit form for T{M) 
This depends on the equations of motion of the theory. 

There exists a general way of identifying horizon entropy in higher derivative theories, 
given by Wald [41] and Visser [42, 43], which makes the dependence of entropy on gravitational 
action very evident. In this approach, horizon entropy is identified as a surface integral over 
the horizon of a (£>— 2)-form corresponding to the Noether charge arising from diffeomorphism 
invariancc of the theory. The basic idea is to find the Noether current J of diffeomorphism 
invariance (which is a closed {D — 1) form), and use the associated Noether potential Q 
defined by = dQ to define black hole entropy. By considering variations (not necessarily 
stationary) of a stationary background spacetime possessing a Killing horizon, Wald argued 
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that the first law of black hole dynamics is consistent with the identification of the following 
expression as the horizon entropy: 



^Noethcr ^ l3 J d^-^S^ J'^ = J d^-^j.^^ jafc (gp) 

where k is the surface gravity, (3~^ = 27r/K is the temperature of the horizon and J", J"'' are 
the Noether current and the Noether potential. In the final expression in Eq. (80) the integral 
is over any (D — 2) dimensional spacelike cross-section of the Killing horizon on which the 
norm of vanishes. For example, in the case of Einstein gravity, Eq. (41) reduces to 

where is the Killing vector generating the horizon. Writing the area element as dJ^ab = 
{la^b — lb£,a)\/^d^'^'^x whcrc la is an auxiliary vector field satisfying la^"" = —1, the integral 
in Eq. (80) reduces to 

^Noether = -(3^ J d'^'^X V^(Za6)VT = /^^'^ J d^^^xV^ = ^Ah (82) 

where Ah is the horizon area. (We have used the relations ^°Va^'' = k£^^ and /3k = 27r.) 
This, of course, agrees with the Bekenstein-Hawking entropy. 

In what follows, we first give a quick derivation (in component notation) of the argu- 
ments of Wald as to why the definition of entropy based on Noether potential makes sense 
for stationary Killing horizons. It turns out that the Noether current J° and it's associated 
potential J"*" (or more accurately, their off-shell versions) characterize most of the results 
related to thermodynamic aspects of gravity, and provide the most natural substitute for the 
equation of motion tensor Eab in attempts to interpret gravitational dynamics as an emer- 
gent phenomenon. So this description is of importance in our later discussions on emergent 
paradigm. 

Let us now sketch the argument leading to an expression for black hole entropy in a 
general gravitational theory with the Lagrangian L[gab, Rabcd]- The main input needed is 
the definition of covariant version of Hamiltonian equations of motion, which can be given in 
terms of a symplectic (D — 1) form cj°, which is defined, in component notation, as 

lo^eeS (5^?;") - (V^ dv") (83) 

where is the Lie variation with respect to a vector field ^, and 5^v°' is defined from the 
variation of the Lagrangian under a diffeomorphism: 

5{L./^ = V^EabSg'^' + V^Va{Siv'') (84) 

In terms of the above symplectic form, covariant Hamilton's equations take the form (note 
that w° is a tensor density of weight unity) 

5H[^] = J d^a ooyV^ (85) 
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where C is a Cauchy surface. 

A rigorous justification for this can be found in [44], which studies the covariant phase 
space formulation for field theories. We here take a digression to simply motivate the above 
expression in analogy with point particle mechanics. In this case, the symplcctic form is 
© = dp A d(7, in terms of which Hamilton's equations can be written as 



dH = ©(..., 9t) 

= q{q,p) dp- p{q,p) dq 



(86) 



The Hamiltonian H here generates time translations (hence the vector dt in the 2nd slot of 0), 
while the "d"'s are arbitrary variations whose coefficients characterize the {q,p) dependence 
of the Hamiltonian. One may re-write the above expression as follows. 



dH = q dp — p dq 

~ d{pq) — pdq — pdq 
= d(pdf(7) - dt (pdq) 



(87) 



Now, the surface term in point particle action is v{dq) ~ pdq. Hence we can write 

dH = d[v{q)]-dt[v{dq)] (88) 

which has the same structure as Eq. (85) with uj"- given in Eq. (83). (Similar conclusion 
would be reached if the surface term is qdp, with the roles of q and p interchanged.) 
Returning now to the main discussion, we sec that a;° can be written as: 



S {V^ 5^v'') - (V^ Sv'^) 



(89) 



This allows us to relate to the Nocther current J", obtained from diffeomorphism invari- 
ance of the Lagrangian, as given by Eq. (40) which is conserved ofF-shell. If we omit from it 
the term involving Eat we get the expression: 



which is conserved on-shell, i.e., when the equations of motion arc satisfied: 



(90) 



/-9 



(91) 



since {Ly^) = y/^Vam"). Therefore, Va J" = when Eab = 0. We wih use the 
definition of J" in the expression for a;° and also utilize the facts: (i) S^v"' = J°[^] + L^"", (ii) 
S{y/^ L) = \/~^ VafJw" on-shell, (iii) (J^-/— g = V^^9^a£."', and (iv) S represents arbitrary 
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field variations, and therefore do not act on ^. Straightforward algebra then leads to the 
result: 



= S {rmV^) + 2V^ Vb [^^-SvO^j (92) 
which can be substituted in Eq. (85) to obtain 

SH[^] ^5 j dllar + j d^a 2Vb (e['^<5«''l) (93) 
c c 

(The steps leading to the above relation are best understood by doing a 3 + 1 decomposition 
w.r.t. normal Ua = ±Ndan to the hypcrsurfacc C, with ± corresponding to spacelikc/timelike 
surface respectively, and using y/—g = N^y\h\ with being the induced metric on C.) From 
the above, one can already note that the flux of Noether current is almost the same as the 
Hamiltonian, except for a surface term. Further, if the variations Sgat satisfy the linearised 
equations of motion, then the perturbed current is also conserved, and one therefore has 
S[J'^y^—g\ = db[S{J"-^y/^)]. Repeating the steps leading to the preceding expression, this 
time in terms of J"'', leads to 

SH[^] - s J dj:aVbr' + J dj:ab&Sv''^ 

C Soo 

= -^S j d^abr^' + ^S Jd^abr''+ Jd^ab&Sv'^ (94) 

5bif <Soo 5oo 

where iSoo is the {D—2) dimensional boundary of C at asymptotic infinity. The inner boundary 
of C is taken to be the bifurcation surface S, where = (hence this surface does not 
contribute as far as the S}-°-5v^'^ term is concerned). It is clear from the second expression 
above that, for a meaningful Hamiltonian H to exist, the last integral must be expressible as 
a variation of a [D—l) form, say B. Assuming this is so, then, as argued in [41], the combined 
contribution to H of this term and the second term involving J""^ above, at asymptotic infinity 
5*00, gives the numerical value of the ADM Hamiltonian when the vector field ^ is chosen 
appropriately to correspond to asymptotic time translation and rotation(s). Now, when ^ is 
a Killing field, from its very construction, = 0, although H[^] itself need not vanish. 

Using all these, one may identify the Soo contributions above as the "work terms" 5£ — VLh5J ^ 
with £ and J7 being the energy and angular momentum respectively (sec [41] for a rigorous 
justification of this). This in turn allows us to identify the black hole entropy S'Noether from 
the iSbif contribution via the relation {K/2Tr)SS ~ 6£ — VLrSJ^. Symbolically, 

- -i (5 j dllab J"' + 1^ J d^ab + J d^ab S.^" 5v''^ (95) 

5bif Soo Soo 



-(K/27r)<5S SE-HhSJ 

Using all the above arguments, Wald eventually identified the entropy as 

•SNocthcr = 9/^ / dSab (96) 
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The above expression for entropy can be rewritten in a simplified (and appealing) form for 
stationary horizons in Lanczos-Lovelock models. However, for sake of simplicity, we sketch 
how this form is arrived at for static horizons. Using Eq. (45) for Jab in Lanczos-Lovelock 
models, we can rewrite Eq. (96) as 

Wher =\P j (2P'^'"='^ V[,Cd]) {'ii[ah\) V^d^^-^)^ (97) 

where la is an auxiliary null vector satisfying and I ■ £, = —1, in terms of which 

dSab = {'^^[ah]) \f5dL^^~'^^x because ^ is the horizon generating Killing vector field. One 
can now use the near horizon static structure of the metric (discussed in detail in section 
Sec. 4.L1; see Eq. (197)) to show that Vf^Cd] = {-2t^n) <5p^(5°] and /'^ = (l/TV^, \/N,Q,{).. .) 
near the horizon defined by n = whence the lapse behaves as iV = Kn + 0{iv'). Using these, 
it is easily seen that 

^Nocthor =\P j (S^Pot) V^d(^-2)^ (98) 

From the total anti-symmetry of the alternating determinant tensor, and the fact that 
= R^B ^ oin"^), it follows that 

P'.: - (y) L^^:? (99) 

which eventually gives 

^Noether = 4^™ j 4f- 1'^ V^d^^^^) ^ (^gO) 

Thus the entropy of m-th order Lanczos-Lovelock theory in D-dimensions can be related 
to the Lagrangian of the (m — l)-th order Lanczos-Lovelock theory in (D — 2)-dimensions. 
The above expression (which holds for stationary horizons as well) provides a description 
of horizon entropy in Lanczos-Lovelock models purely in terms of intrinsic geometry of the 
horizon, without any reference to any additional vector fields etc. For a more general (non 
Lanczos-Lovelock ) Lagrangian, this, of course, will not be the case. In particular, the fact 
that entropy is purely a functional of intrinsic horizon geometry can presumably be taken 
as another indication of 'holographic' properties of Lanczos-Lovelock models discussed in 
Sec. 2.7. 

Two other aspects of Noether charge entropy worth mentioning in this context are the 
following. First, in Einstein gravity, we know that the York-Gibbons-Hawking boundary 
term, evaluated on the stretched horizon Hg, also gives the correct horizon entropy in the 
limit e 0. A natural question, therefore, is whether the Lanczos-Lovelock surface term 
given in Eq. (56) also yields the same entropy. This definition of entropy was first given by 
Zanelli and Teitelboim [23], and an argument by Iyer and Wald in [41] shows that it matches 
with the Noether charge entropy. However, an explicit correspondence between these very 
different forms of the two expressions for entropy is not known in general. Perhaps the 
structural properties of the surface term described in detail in Sec. 3.2 and Sec. 3.3 might be 
useful to establish a direct correspondence between the two entropy expressions (at least for 
maximally symmetric horizons) . (It is also possible to obtain the entropy of the horizon by 
a procedure (called 'replica procedure') related to the conical singularity and deficit angle in 
the Euclidean sector; this is discussed, e.g., in Ref. [45].) 
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Second, there has been extensive work, based on complex path formahsm [48], to interpret 
the Hawking radiation as a tunnelhng phenomenon (see e.g., [49, 50] and also [51]) and the 
question arises as to the status of this approach in Lanczos-Lovelock models. In this context, 
S. Sarkar and D. Kothawala were able to prove the following result for spherically symmetric 
horizons [52]: For any theory of gravity in which field equations near the horizon can be 
written in the form TdS — dE ~ PdV, the tunnelling formalism of [51] yields an entropy 
which matches the Wald entropy once a specific choice is made for the Hamiltonian used in 
the tunnelling formalism. As we shall show in section Sec. 4.1.2, the above criterion is satisfied 
near any static horizon (not necessarily spherically symmetric) in Lanczos-Lovelock models, 
with S being the Wald entropy (see Sec. 4.1.2 for definition of E etc.), and hence the result of 
[52] implies that the tunnelling mechanism can reproduce the correct expression for entropy 
in Lanczos-Lovelock models as a simple consequence of the first law of thermodynamics. 

We conclude this section with a brief mention of an approach to the entropy of horizons, 
which does not seem to generalise from Einstein's theory to Lanczos-Lovelock models easily. 
This approach is based on entanglement of quantum fields. The modes of the quantum fields 
outside the horizon are entangled with modes inside and when the latter are traced out, 
the physical processes outside the horizon is governed by a density matrix p oc exp(— /3i?). 
A natural entropy associated with this density matrix is 5 = —Tr{p\ogp). This quantity 
has a divergent behaviour of the form lime_^o(^/4e^) where A is the horizon area and e has 
dimensions of length. If one cuts-off the limit at Planck length, one can match the expression 
with the standard result in Einstein gravity with entropy coming out as proportional to area. 
(This cut-off can be introduced in a covariant manner, using, for e.g., the techniques discussed 
in [46].) But since the entropy in Lanczos-Lovelock models is not, in general, proportional 
to area it is not clear whether the idea of entanglement entropy generalizes beyond Einstein's 
theory. (If it does not, it raises the possibility that it may not be the correct approach even in 
Einstein's theory.) The fact that it is divergent makes it difficult to draw firm conclusions. It 
is, for example, possible that the regularization procedure to introduce the cut-off somehow 
leads to the correct expression [47] but these ideas have not been worked out completely. So 
it is fair to say that whether entanglement entropy approach reproduces the correct Lanczos- 
Lovelock entropy is at present unclear. 

3.2 Entropy of self-gravitating systems and Wald entropy 

It is well known that statistical mechanics of systems with long range interactions exhibit 
various peculiar features (see, for a review, Ref. [53]) which could lead to subtleties in the 
thermodynamic description of such systems. It is therefore natural to ask whether effects of 
self-gravity can alter the conventional (extensive) nature of entropy of a matter configuration 
that is undergoing gravitational collapse, when the matter configuration is on the verge of 
forming a black hole. 

This issue was first raised by Israel et. al. [54] in an attempt to give an operational 
definition of black hole entropy in Einstein gravity. Their results reproduced the well known 
area scaling of entropy in the limit of a matter shell collapsing to form a black hole (probably 
not a surprising result for a single shell), but more importantly, gave the correct numerical 
coefficient of 1/4 in the near horizon limit. This then motivated Oppenheim [55] to study a 
configuration of closely packed shells and ask how exactly the conventional volume scaling of 
entropy transmutes to an area scaling as the horizon is approached. The results obtained in 
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these works (also see [56]) suggested that entropy scahng approaches the Bekenstein- Hawking 
area law for self-gravitating systems provided the shells are assumed to be instantaneously in 
equilibrium with the local acceleration temperature in a manner made precise in the discussion 
below. 

This result seems too good to be true at a fundamental level, since it seemingly provides 
an explanation for origin of black hole entropy without any appeal to quantum gravity! 
Perhaps this is just what happens, but one must be cautious since there are several ways 
of obtaining the correct entropy scaling using very different arguments. The contribution of 
effects of self-gravity to horizon entropy can be more firmly established if we could test the 
above analysis for theories other than Einstein gravity, since that would help isolate dynamics 
from kinematics in the final result. In this section, we shall describe the re- analysis of this 
approach for Lanczos-Lovelock models, and show that several non-trivial differences arise. 
The conclusion is that horizon entropy and entropy due to effects of self-gravity are different 
from each other in general, although the difference itself turns out to have an interesting 
form. This is one example in which the study of Lanczos-Lovelock models allows us to 
clearly isolate certain numerical coincidences which occur in Einstein's theory and do not 
carry over to Lanczos-Lovelock models. 

To calculate the standard thermodynamic entropy of a self gravitating system, we will 
basically follow the set up suggested by Oppenheim [55], but in a slightly different manner 
to bring out the key features more prominently. We consider a system of n densely packed 
spherically symmetric shells in D spacetime dimensions, assumed to be in thermal equilib- 
rium and supporting itself against its own gravity. We shall be interested in the entropy of 
this system when the outermost shell is close to the event horizon of the system. We first 
describe a general set-up to study the system described above for any spherically symmetric 
spacetime. This is important for two reasons. First, as we shall see, it will highlight the 
key mathematical feature responsible for entropy density S/A of the system being 1/4 for 
Einstein's theory in arbitrary dimensions. In [55], this result was obtained for D = 4 by 
explicitly using the Schwarzschild form of metric functions and the expression for Hawking 
temperature. However, we shall not require any explicit expression for metric functions or 
Hawking temperature to obtain the area scaling. Second, the set-up we describe can be used, 
without any modification, to obtain entropy density for higher curvature gravity actions. 

Let us denote the variables describing the ith. shell with a subscript i. Since the system 
is spherically symmetric, we can write the metric outside the zth shell as 

ds^ = -Cifi{r)dt^ -f h{r)-^dr^ + r^dn^ (101) 

where dfl^ is the metric of a unit {D — 2) sphere. At the location of the shells (given by 
Ti =:const), the above metric must satisfy the Israel junction condition which states that 
the induced metric on the hypersurfacc should be continuous. This leads to the following 
conditions on the constants cf. 

Cifiin+i) = Q+i/i+i(ri+i) (102) 

and we shall choose Cn = I (which fixes the interpretation of t as proper time for asymptotic 
observers). Then Eq. (102) can be solved to give 

fk+i{rk+i) /n-i(r„-i) ./„(r-„) 
Ck = — TT ^ 7 7 77 7 — 7 (103) 

Jk{rk+i) /n-2(r„-i) /„_i(r„) 
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Note that when fn{rn) = — i.e., when the outermost sheU is exactly on the horizon — we 
have c„ = for all (i ^ n) which implies that 

500 = V i (104) 

Here the condition that the goo's vanish even for the inner-shells indicates that our assumption 
regarding the staticity of the inner shells is not valid when /„(r„) is exactly zero. (This is 
because we know that a particle cannot be kept at a fixed position inside a black hole.) 
However, for the purpose of our discussion, we only need to consider the limit in which the 
outer shell is very near to the horizon, that is, /«(?"„) — >■ 0, then the result that the of 
Ci's (other than c„) and goo vanish is just the leading order result in this approximation. 
Henceforth, we shall assume that wc arc working in this limit and will not state it explicitly 
unless otherwise needed. 

We shall next use the second junction condition which relates jumps in geometric quan- 
tities across a shell, to the matter stress tensor of the shell, Let the normal to a 
r =const. surface (dropping the subscript i for convenience) be Ua = b^^^^dar, so that 
the induced metric on r = constant surface becomes /i^^ = {gab — naUb) 51^5^ in coordinates 
(t, 9a)- Here, we will use Greek letters to define components in the r = constant surface with 
coordinates {t^d^). Within the hypersurface, we can further define il^ = y/cf d^t, and the 
induced metric on the level surfaces of g^j/ = h^^, + Uf^Ui,. It then follows from spherical 
symmetry that tf^i, for each shell has the general form 

= PiU(.i)^,u^i)^ + Piq(i)^,j (105) 

with Ei ~ Aurfpi is the energy of the shell. The physical interpretation of Ei and Pi is that 
of energy and pressure as measured by a local observer at rest on the shell. Further, the 
condition for thermal equilibrium in curved spacetime implies [57] 

Ti\J- glair i) = Too 

where T^o is the temperature of the system as measured by a static observer at infinity, and 
/Xi's denote the chemical potential. We can thus express all the T^'s and /x^'s in terms of just 
two unknown parameters Too and fin respectively. In thermodynamic equilibrium, each shell 
satisfies the Gibbs-Duhem relation 

E,=T,S.,- P^A, + fi,N, (107) 

where Ni is the number of particles composing the ith shell and Ai is the area of the ith 
shell. Using Eq. (106) in the above expression, we can write the total entropy S'mattcr of the 
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system in the form 



E 

i 

E 



rp y 900\ri) rp 



Tf, \/c^fi{r,) - fn(r 



(108) 



Now assuming that /i„-/V is a finite quantity, the last term in the above expression vanishes 
in the near horizon hmit, since 5oo(''") vanishes (see Eq. (104)). Hence the only non-trivial 
contribution to Smatter can come from the first term, which we wish to evaluate in the limit 
when outermost shell is close to the horizon of the system, r„ — > vh (where fnirn) = 0); we 
shall refer to this as the near-horizon limit. We can read-off this contribution as follows. 

First, we note that (Ei + PiAi)/T,^ has no dependence on since does not depend 

on Ci (see below). Further, in the near horizon limit, Ci = Vi 7^ n, c„ = 1, from which it 
is easy to see that the only non-zero contribution to S'matter will come from the i — n term, 
and in particular from those terms in (En + Pn-An) which diverge as l/-\/ fn{rn)- Therefore, 
our strategy to calculate the entropy of the system will be to analyze the form of En and 
Pn and look for such divergent factors. As we shall now demonstrate, this strategy provides 
an extremely quick way to obtain the result in Einstein theory (compared to the explicit 
computations in [55]), and also facilitates easy generalization to higher derivative theories. 

In Einstein theory, the surface stress tensor is given by 

8^4S^ = {^KK,-K,}^ (109) 

where (. . .); = (. . — {■ ■ ■)i evaluated at r = r^. The extrinsic curvature of r =const 
surface, for the metrics described by Eq. (101)), is given by (after dropping the subscript i 
for convenience): 



K^, = -[Vbf'/2f^u^u,+ (^Vb/r^q^, 



K = hf^^K,,, 

= [Vbr/2f)+{D-2){V~b/r) (110) 

Note that, as stated above, there is no dependence on c,;. Based on our strategy outlined 
in the previous section, we must now find terms which, for i = n, diverge as 1/ \/ fni^n) as 
fn — >■ fH- We now specialize to the case bi{r) — fi{r) Vz. (This condition is automatically 
satisfied, in the near-horizon limit, for all the cases in which the horizon has no curvature 
singularity; so this choice entails no loss of generality.) Using Eqs. (109), (110) and (108), we 
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can immediately write, for the divergent parts, 



87r[p„]div = Stt 



(E) 
(n)05 



8^[Pn]div 5^ = 8^ 



{E)A 
(n)B 



in the hmit 



div 
(n) _ 



(111) 



r-n th, fnirn) (here = K^^u'-'u"). We have recalled the definition 
of Pi's and P^'s from Eq. (105). Hence, we see that the main divergence comes from 

l-\ hm 

167r. 



"J div 



1 



lini — , 

4 ^-^0+ Vfn{rH+e) 



(112) 



where we have used the definition of Hawking temperature for a system which has collapsed 
to a black hole. Using the above result in Eq. (108), wc finally get 



lim 

e->0+ 

Ah 

4 



(l/4)To./V/„(rg + e) ^ /7^— T 
Ah X x/fnirn + e) 



(113) 



Hence, we see that the entropy of the configuration, as the outermost shell approaches the 
event horizon of the system, is dominated by the Bekenstein-Hawing entropy. The crucial 
factor of 1/4 comes from the last line of Eq. (112), which can be recast as 



[Pn]dw/Tn = 1/4 



(114) 



Interestingly, this same relation arises in the context of a derivation of Navier-Stokes equation 
and its relation to emergent gravity paradigm [58] . 

In what follows, we shall repeat the above analysis for Lanczos-Lovelock gravity, and 
show that the correspondence between shell entropy as calculated above, and the black hole 
entropy as given by Wald formula does not always hold. We will also point out that our result, 
nevertheless, has a curious interpretation in terms of the "equation of state" of the horizon 
degrees of freedom, and also has a very direct mathematical connection with the membrane 
paradigm for black hole horizons. 

As discussed in Sec. 3, the surface stress tensor in the m"^ order Lanczos-Lovelock theory 
is given by 



iTT t 



2m+l 



s{m) \^s(m) J 



■M2,7i-l] •^l'^2 f2s-ll'2s f2s+l f2™ 



(115) 
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where the coefficients Cs(m) ^iro given by 



~ g!(2m-2,-f)!! ^^^^^ 

Here R^ipau etc. indicates that these quantities arc to be evaluated for the induced metric 
h^iu- To calculate the matter entropy of the system of n shells, we proceed in the manner 
similar as in the case of Einstein's gravity described in the previous section, that is, by finding 
. We note that, in the near horizon limit (and assuming that the intrinsic Riemann 

div 

tensor of horizon surface is finite), the only divergence in the stress tensor will come from the 
A'P component of the extrinsic curvature of the r = r„th shell, which diverges as (l/\/7n) 
as can be seen from Eq. (110). The transverse components Kg vanish as {VTn/r). Further, 
note that the determinant tensor has the property 

te::;: = ^if'k x {szs^^ . . . szs^^) (117) 

That is, the presence of in each row of the determinant tensor forces all the other indices to 
take the values 2, 3, • • • (D — 1). Thus, at most, we can have only one Kq present in any term 
of Eq. (61) which shows that the maximum possible divergence is 0(l/-\/7n) even though the 
surface stress tensor is a (2m — l)th degree polynomial in K. From the structure of Eq. (61), 
it is easy to see that only the s = (to — 1) term gives this divergent contribution while all the 
other terms being of 0(v7) or higher vanish while from Eq. (117), it is again obvious that 
only the transverse component of the surface stress tensor will contribute. Thus we get 

where in the second line wc have used Cm-i — 4/(to— 1)! from its definition in Eq. (58), and 

replaced with the intrinsic curvature Rq^ defined completely in terms of the induced 
metric qab of the constant t, constant r surface. Such a replacement is valid since the 
corresponding extrinsic curvature vanishes. We can now obtain the pressure easily using the 
fact that QAB is maximally symmetric (see Eq. (105)), and hence 

rp 1 _ 1 „AB r, ] _ 2TOQ;m 1 1 ^AAi---A2^-2 £,32,^-382^-2 7^6 

2TOa™L,(f-) 4 (119) 

In arriving at the second equality, we have used the following two relations concerning 
Lanczos-Lovclock actions of order to in Z? dimensions derived earlier in Eqs. (30), (37) 

pi _ riaitil-ambm DCidi JDCmdm 

^kni) = 5 (120) 
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where -E'j(m) is the equation of motion tensor for Lanczos-Lovelock action. Hence, we once 
again find that the only non-zero contribution to matter entropy comes from the pressure 
term as in the case of Einstein gravity (see Eq. (111)). No further computation is needed 
since we can simply obtain the entropy by comparing expression in Eq. (119) with the first 
of Eqs. (112) in Einstein theory. This immediately allows us to read-off the matter entropy 
as 

^i™Lr = [^^) 4™™4f-^)A« (121) 

where Ah is the D — 2 dimensional hypersurface area of the horizon. On the other hand, 
the Wald entropy for general Lanczos-Lovelock theory obtained earlier in Eq. (100), when 
evaluated for the spherically symmetric case is given by: 

^Waid - 47rma™i|,f_f)AH (122) 
which leads to the following relation 

'-'matter — I JJ _ 2 ) Wald l^^'JJ 

So the matter entropy in Lanczos-Lovelock theory is proportional to the corresponding Wald 
entropy of the black hole. We also see that ^matter < 'S'waid with the equality holding only in 
Einstein's gravity m = 1. This inequality strongly suggests that, in general, the microscopic 
degrees of freedom responsible for the entropy of the black hole are certainly not the extrinsic 
(matter) degrees of freedom forming the black hole. Further, since a general Lanczos-Lovelock 
theory will be described by Lagrangian of the form L = am,Lm, the total matter entropy 
in such theories would be 



lD-l)/2] , 

D — 2m 



E 



D 



^tid (124) 



which has no simple relation, or proportionality to the Wald entropy 

[(D-l)/2] 

•S'wald = E ^wld (125) 

m=0 

We highlight another point which is relevant to the membrane paradigm for Lanczos- 
Lovelock models, to be discussed in Sec. 3.3. As is evident from above analysis, it is only 
the pressure term which leads to a non-zero contribution to the matter entropy and hence, 
in the near horizon limit, we have the result: 

'5'matter [Pn] div 



Ah Tn 

P , ^ 

= 7^ (126) 

where for the sake of brevity, we have replaced [P,i]div by P. Thus the scaling of entropy 
with the horizon radius rn, as well as the proportionality constant, is completely determined 
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by the relation between P and T, which is to say that the equation of state P = P{T) 
determines the matter entropy. In the case of Lanczos-Lovelock gravity, we can read off the 
equation of state P = P(T) = [t^^^^Jdiv from Eq. (119) by using the relation Kqq = —A'? = 
-I'h/M) = -27rr„ to get 



P /D-2m 



4™,„L|f_i') (127) 



Tn \ D-2 

Using Eq. (123), the above relation can also be written in terms of Wald entropy as: 

PAh f D-2m \ (,„) 

In fact, the pressure P obtained above is exactly what one would call as surface tension of 
the stretched horizon in the membrane paradigm. This is mathematically obvious from the 
way we arrived at the near-horizon for P, but it must be emphasized that the equivalence 
with membrane pressure is strictly valid only in the near-horizon limit. 

The main implication of the above result is that matter entropy of a system is not neces- 
sarily the only contribution to entropy of a black hole formed from collapse of that system. 
This unambiguously establishes the fact that there might be more to horizon entropy than 
just the effects of self-gravitation, at least to the extent that the approximations made in 
the above analysis are valid. (We have not studied the system in full generality, since the 
dynamics of the collapse has not been taken into account. We only considered the shells at 
various stages of collapse, imposing consistency conditions in the form of junction conditions 
at each instant of the collapse. While this may be sufficient, it will be interesting to study 
what happens in a fully dynamical situation.) However, whatever be the nature of additional 
corrections introduced by a more general analysis, it is unlikely to establish the complete 
equivalence of the resultant entropy with the Wald entropy. 

It is also interesting to view the above result in the light of several entropy bounds 
on matter entropy which exist in the literature. The model calculation above gives precise 
expressions for matter entropy and its relation to Wald entropy, rather than just an inequality. 
Finally, the analysis above also highlights a possible connection between horizon entropy 
and the "equation-of-state" P(T) satisfied by the degrees of freedom residing on the horizon 
surface. This would become more evident when we present a study of the membrane paradigm 
for the Lanczos-Lovelock models in the next section, deriving all the transport coefficients 
to linear order in the perturbation theory. 



3.3 Membrane paradigm in Lanczos-Lovelock models 

We now use the mathematical tools of the previous section to derive several properties of the 
so called stretched horizon in a theory described by the Lanczos-Lovelock action. That is, 
we set up the so called membrane paradigm for the Lanczos-Lovelock models. After a brief 
introduction to the membrane paradigm in Einstein gravity, we shall discuss the analysis 
for the Lanczos-Lovelock models, keeping in focus the new features introduced by the non- 
triviality of the Lanczos-Lovelock action. 

Specifically, we shall present a derivation of the various transport coefficients associated 
with the membrane to first order in perturbation theory, and demonstrate how the structure 
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of these coefficients brings into focus several fundamental features of the membrane paradigm 
itself, which are otherwise not evident when attention is restricted solely to Einstein gravity. 
This reinforces our original point of view about the utility of Lanczos-Lovelock models in 
gaining insight into aspects of Einstein gravity itself. 

The membrane paradigm of black holes [59] takes the point of view that, as far as the 
static observers outside the horizon are concerned, the black hole horizon can be replaced 
by a stretched horizon, a membrane, endowed with specific physical properties which encode 
the presence of the inaccessible black hole region. Although originally developed to facilitate 
inclusion of the physics of black hole background in astrophysical applications, work over the 
last decade or so has indicated that physical properties of the membrane (such as viscosity 
etc.) might have a deeper relevance, from the point of view of holographic dualities which 
map gravitational systems to non-gravitational systems in one lower dimension. This has led 
to a renewed interest in the study of membrane paradigm. 

The membrane paradigm in Einstein gravity has also been analysed from a view point of 
the emergent gravity paradigm. The Damour-Navier Stokes equation governing the dynam- 
ics of the black hole membrane could be obtained [58] starting from an action which could 
be given a thermodynamic interpretation as an entropy production rate when expressed in 
terms of thermodynamic variables such as temperature, entropy, pressure, etc of the hori- 
zon. This comes very close to quantifying a precise connection of properties of the horizon 
membrane with the thermodynamics of horizons, but as we shall show in this section, a more 
direct connection between the transport coefficients of the membrane and its thermodynamic 
properties can be established by studying the Lanczos-Lovelock models. 

In Einstein gravity (as well as with the Gauss-Bonnet term added to it), the geometric 
structure of the fluid transport coefficients of the membrane is not evident from their algebraic 
expressions. For example, in Einstein's gravity, it is known that the shear viscosity 77 is equal 
to l/(167r). But the origin of this factor is not well understood in terms of the geometric 
properties of the horizon such as the intrinsic curvature of the horizon etc. As we shall show, 
setting up the membrane paradigm [60] for the Lanczos-Lovelock models once again helps 
us in not only understanding such geometric aspects of the transport coefficients in a more 
general manner, but even explains the origin of factors such as l/(167r) for shear viscosity in 
Einstein theory! 

Let us briefly describe the geometric set-up and the procedure to construct the membrane 
paradigm. More details can be found in [59], and we only refresh the essential expressions 
below to set up the notation etc. We shall work in a D dimensional spacetime containing an 
event horizon H generated by the null geodesies I. The surface gravity k is then determined 
through Vii = kI. We shall deal with a stationary background spacetime, in which case I 
corresponds to the timelike Killing vector ^ restricted to H where = 0. For asymptotically 
flat spacetimes, the norm of I can be fixed at infinity and there is no ambiguity in the 
definition of k. The stretched horizon Hs is then defined as a timelike surface infinitesimally 
close to H, and separated from H along spacelike geodesies n. This surface is spanned by a 
unit timelike vector field u and (D — 2) vectors e^^-j's. The correspondence between points 
on H and "Hs is made via ingoing null geodesies k normalized to have unit Killing energy, 
k • I = —1. The horizon H is then located at A = 0, where A is the affine parameter along fc, 
related to the norm of ^ by A = —{1/2k)^^. We shall define N = \/—^ to simplify notation 
and also to facilitate comparison with known results in the literature. In the limit A — > 0, 
H and u ^ N-'^l and n N-'^l. 
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The induced metric hf^^ on Hs is given by h^^i, = hate'^^-^e'l^^y where hat = 9ab ~ naUi, 
and e(^)'s are basis vectors spanning Jig (e(^) • n = 0). Similarly, the induced metric jab 
on the (Z? — 2)-dimensional space- like cross-section of Hs orthogonal to u is given by 'jab = 
7(j{,e°^^ e^^^ where -fab = hab + UaUf, and e^^-j's are basis vectors spanning 7 (e(^) • u = = 
£{A) ' fi)- The extrinsic curvature of Hs is defined as /•sT^jy — e^^^e^^^ Van.f, and it is easy to 
verify that in the limit N ^ 0, we have: 



Kai 



^tJ.v^-[Ay-(B) 

N~'kAB, (129) 



where kAB is the extrinsic curvature of the {D — 2)-dimensional space-like cross-section of 
the true horizon H. This can be decomposed with respect to its trace- free part as 

kAB ^ CTAB + JJ^^^^OjAB, (130) 

where j^^ctab = 0, and 9 and aAB arc the expansion scalar and shear respectively of the 
I congruence generating the horizon. For a Killing horizon, both of these vanish as N'^ and 
therefore Kqq remains the only divergent contribution to K^^. However, since we wish to 
investigate dissipative properties associated with the horizon, we shall perturb the background 
away from staticity and then obtain the membrane stress tensor to first order in perturbations 
to determine the transport coefficients. The extrinsic curvature can be written, upto linear 
order in perturbations, as 

(0) (1) 

K^^K^+K^ (131) 

(0) 

where is the unperturbed extrinsic curvature. Since we assumed the background geometry 
to be static, we have 

(0) (0) (0) 

i^oo = -^/N, K^A^O^ KAB^{N/r)6AB (132) 
To study the additional divergences brought in when the horizon is perturbed, it is convenient 

(0) 

to choose perturbations which do not affect A'qq. (As we shall see, our final relations will 
only have divergence of 0{N~^), and perturbing k will not lead to any qualitative difference.) 
Formally, since we assume the background k to be non-zero, the perturbation in k is not gauge 
invariant and our choice corresponds [59] to a gauge wherein Sk = 0.) Therefore, we have 

(1) (1) (1) 

= 0, A'a^^O, KAB^{l/N)kAB (133) 

where we have not bothered to put a "(1)" over kAB since its unperturbed part is strictly zero. 
The expansion and shear parameters of the perturbed horizon will now introduce additional 
divergences in the stress tensor via Kab- 
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We are now ready to use the above facts to analyze the membrane stress tensor in Lanczos- 
Lovelock theories of gravity. For further details about the membrane paradigm and aspects 
of the perturbation scheme, we refer the reader to earlier literature on the subject [59, 61], 
and also to the introductory sections in the recent paper [62] which focuses on Einstein- 
Gauss-Bonnct theories. Since the derivation for Lanczos-Lovelock involves considerable 
combinatorics, we present a summary of the final results first and then provide a sketch of 
the derivation [60]. 



3.3.1 Membrane paradigm: summary of the results 

We shall be mainly interested in membranes which have isotropic tangential stresses (pres- 
sure), which is only possible when all directions, everywhere within the horizon surface, are 
equivalent. (That is, the horizon surface is maximally symmetric). Therefore, the cross- 
section of the horizon of the background geometry has the curvature 



iiABCD 



(/C/r^) {"{AClBD - lADlBc) 



(134) 



where /C = ±1 and r is a constant of dimension length which is related to the intrinsic Ricci 

(0) 

scalar of the horizon cross-section as K,{D — 3){D — 2)r^^ = ^^^^^ R. For simplicity of 
notation, we will assume /C = +1 in the intermediate steps, and only restore it in the final 
expressions (which is easy to do). Doing so will come in handy while discussing the case of 
planar horizons, for which we will simply put JC = 0. With this assumption, it can be shown 
that the horizon membrane has the stress tensor 



with p 



tap 



Pressure : p 
Energy density : 

Shear Viscosity : 

Bulk Viscosity : 



^ PsUaUfi + e\f^e-p {PslAB - 2r]s(TsAB - CsSs^ab) 

etc., and the different variables are given by: 
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where we have introduced the notation Dk = {D — k) to avoid clutter, and a„i is the coupling 
constant of the m*^ order Lanczos-Lovelock Lagrangian 
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The transport coefRcients above can be connected with the thermodynamic properties of 
the horizon in the Lanczos-Lovelock gravity, by noting that the Wald entropy and quasi-local 



40 



energy of the horizon in these theories are given by 



J H 



^Wald 

This immediately imphes that 



^waid = l\xj^d^L^-^ (138) 



with 



m— 1 



Too w^o ^\D~0 

m— ^ 

Similarly, it is obvious from the dependence on L^z\ and L^~^ factors appearing in the 
third and fourth equalities in Eqs. (136) for the shear and bulk viscosities, that both 77 and 
C are also expressible solely in terms of and E^"'^\ 

3.3.2 Membrane paradigm: details of the analysis 

As discussed in Sec. 2.6, for the w}^ order Lanczos-Lovclock Lagrangian Lm in Eq. (20), the 
surface stress tensor is given by (see Eq. (61)): 

, m— 1 

^ 6=0 

[We have replaced Cm in Eq. (61) by for notational clarity.] One can check that for to = 1 
and m — 2, the above expression reduces to that of the surface stress tensor in Einstein [26] 
and Gauss-Bonnet [27] theories respectively. 

As in the case of the membrane paradigm for Einstein gravity, we will interpret the tap as 
due to a fictitious matter source residing on the stretched horizon . Using the assumptions 
and expressions given above we can show that the surface stress energy tensor can be written 
in a form analogous to that of a viscous fluid, namely, 

tal3 = PsUaUf3 + 6^^' 6^^'' {pslAB ~ 27]s(TsAB ^ Csf^sjAs) (142) 

and thus read off the corresponding energy density Ps, pressure p^, shear viscosity co-efhcient rj 
and co-efficient bulk viscosity Cs in terms of the geometrical quantities describing the horizon. 

From the fact that Tf(s)'^ in Eq. (61) is a polynomial of odd degree [2(to — s) — 1] in i^T^jy, it 
is easy to see that, to first order in perturbations, the only divergences in t^p are of 0{N~^), 
and they arise from s = (m — 1) and s = (m — 2) terms in the series. This follows (upon 
simple counting) from Eqs. (132), (133), (61) and (117). All the other terms in the series are 
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of 0{N) and hence vanish in the limit N 0. We give below an outline of the contribution 
of divergent terms, stating only the relevant expressions and skipping the intermediate steps 
which arc algebraically straightforward. We will also define 
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(143) 



A. Divergence due to the zeroth order term: 

Only the s = (m — 1) term will contribute to the zeroth perturbative order. The 
corresponding contribution [60] turns out to be 
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B. Divergence due to the first order term: 

At first order in the perturbations, we get divergences from s ~ (m— 1) and s = (m — 2). 
The contributions from these can also be found in a straightforward way. and are given 

by 
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where the symmetry factor of 6 in the last equation comes because we need to choose, 

(1) 

from 3 factors of Kap, one factor of K ab (3 ways) and one factor of Kuu from the 
remaining 2 Kap^s (2 ways), which can be done in a total of 3 x 2 = 6 ways. 

The above equations can be further simplified by using the following identities (the proofs 
for which are given in Appendix C) 
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where we have again used the notation Dk = {D — k) introduced earher. We can now collect 

/ -IN • 

the 0[N ) terms in the total stress tensor tap = t ap+ t upto first order in perturbation 
to get 
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We have thus accomplished our main task of finding the divergent contributions to sur- 
face stress tensor i^i, of the membrane to linear order in horizon perturbations. The only 
divergence is of 0{N~^), which is precisely the same as in Einstein theory, and can therefore 
be regulated in the same way as is done in conventional membrane paradigm for Einstein 
gravity. That is, we simply multiply by to obtain the finite part of t^^, with the rcgular- 
ization having the interpretation of canceling the infinite redshift at the horizon. Writing fc^ 
in terms of the traceless part and trace 9 as in Eq. (130), we can express the stress tensor 
in the required form 



and read off the pressure, transport co-efficients, etc as 



(151) 



Pressure : p^™-* 
Energy density : p["^^ 



Shear Viscosity : 77 
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where, again, Dk = {D — k). One can check that for Einstein's gravity m = 1 and Gauss- 
Bonnet gravity m = 2, the above expressions reduce to the corresponding expressions in the 
literature [59, 62] . But of course, as is evident, the final expressions yield much more than just 
a generalization of conventional results from Einstein theory to Lanczos-Lovelock models. 
As we already indicated earlier the above expressions highlight a connection between the fluid 
transport coefficients and thermal properties of the horizon and calls for a more thorough 
investigation into possible physical implications of the results. 

In particular, one quantity which has gained significance from the point of view of the 
AdS-CFT correspondence is the ratio of shear viscosity to entropy density, rj/s. Using the 
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expression of Wald entropy given in Eq. (138) in above expressions, we arrive at 
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where we have used the relation L^^ ^ = (/C/r^) D2rnD2m+iL^i-\ which holds for a (D — 2) 
dimensional maximally symmetric spacetime. We can therefore write 
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47r 



1 ^ C^ii'^rn 
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1 + I] Sm 

m=2 
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where nic = [{D — l)/2] and :— Sm/si — Asm- 

One can now investigate the ratio ry/s on a case-by-case basis, considering specific solutions 
with horizons in Lanczos-Lovelock theory. Specifically, for planar ^ horizons, /C = 0, and we 
note that in this case, C,„ — > (• ■ ■)IC^^ ■ Therefore, since s„i — > (. . .)/C™^^, only the m ^ 2 
term in the sum in the numerator survives for /C = 0. One is therefore led to the conclusion 
that no Lovelock term other than Gauss-Bonnet contributes to the rj/s ratio for the planar 
(/C = 0) case; this contribution is given by 



For /C = 



1 + [C2S2] 



K=0 



1 

47r 



4tt 

1-4^2(^-4) I — 



(155) 



Although only the Gauss-Bonnet contributes for the planar case, it must be noted that rj/s 
will, in general, depend on other Lovelock coupling constants as well, through the ratio {k/th) 
which must be calculated using a specific /C = solution of the full Lovelock action. 

A similar result has been noted in the literature [63, 64] in the context of AdS-CFT 
correspondence with the difference being in the boundary used — which is the AdS infinity 
whereas we have considered the boundary to be the stretched horizon. However, it has been 
argued in [65] that the ry/s ratio for the two cases are related by a renormalization-group 
flow equation in Einstein's gravity. It would indeed be interesting to check whether such 
a flow exist in Lanczos-Lovelock theory, since our result for the membrane paradigm, in 
conjunction with the AdS-CFT results for Lanczos-Lovelock gravity, does indicate such a 



couple of clarifying points while dealing with planar horizons arc in order: (i) although {ktu) would in 
general depend on K, we expect it to be finite for K = 0, and (ii) the most general form of a planar horizon 
metric, consistent with the assumption of staticity, is: ds|j = fl(r/Li) (^dx^ + dy^ -|- . . .), where fl(r/Li) is a 
dimensionless function constructed out of ratios of the coordinate r and any other length scale(s) Li appearing 
in the metric. For known solutions which arc asymptotically AdS, Q(r/Li) = r/£, £ being the AdS scale. 
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possible connection in the /C = case. Further, in the special case of a black brane solution in 
Einstein-Gauss-Bonnet theory, we have n/rn = {D — l)/2£'^, in which case the result above 
matches exactly with the one obtained in the context of AdS-CFT (see [62] and references 
within). Using the expression for rj/s, it would be interesting to check whether it obeys or 
violates the KSS bound [66], rj/s > l/(47r), by explicitly computing it for black hole solutions 
known for Lanczos-Lovelock gravity. For JC ~ case, the KSS bound is found to be violated 
in the AdS-CFT context in Einstein-Gauss Bonnet gravity [67, 62] and also in a general 
Lanczos-Lovelock theory of gravity [63, 68]. 

We can now summarize the broad implications of the results discussed in this (and the 
previous) section. 

1. The connection between various transport coefficients and the thermodynamic prop- 
erties of the horizon, such as temperature, entropy and quasi-local energy, becomes 
immediately obvious in our analysis, and it would be interesting if such a connection 
holds for horizons in more general class of gravity theories. 

2. The calculation highlights the connection between Wald entropy S'waW of the horizon 
and the pressure Ps of the horizon membrane. This relation is encoded in the form of 
an equation of state. 



The algebraic steps leading to the above result are precisely the same as those leading 
to the evaluation of entropy of a self-gravitating configuration of densely packed shells 
on the verge of becoming a black hole, discussed in the last section. 

3. The relation between bulk and shear viscosities for Lanczos-Lovelock order m, C^™^ /j/™) 
— 2(D — 3)/(D — 2) is also noteworthy. This ratio is independent oi m, which immediately 
implies 

for an action which is a sum of Lanczos-Lovelock terms. This corroborates the results 
already known for Einstein gravity, and recently also established for Einstein- Gauss- 
Bonnet case [62] . It would be worth investigating whether there is a deeper reason for 
the robustness of this ratio. 



3.4 Second law for black hole entropy in Lanczos-Lovelock models 

Having analyzed the Noether charge entropy in Lanczos-Lovelock models from various dif- 
ferent perspectives in the previous sections, we now discuss the issue of whether the change 
in this entropy can be shown to be positive definite; that is, whether the second law can be 
shown to hold during a physical process for black hole entropy in Lanczos-Lovelock theory. 

In general relativity, the "area theorem" implies that area of a black hole cannot decrease 
in any physical process, so long as the matter stress tensor satisfies the null energy condition 
[69]. However, the equilibrium state version of first law for black holes, which has been 
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established by Wald and collaborators, says nothing in general as to whether the Wald entropy 
always increases under physical processes in arbitrary higher curvature theories. In fact, black 
hole entropy is no longer proportional to area in arbitrary higher curvature theories in general, 
and Lanczos-Lovelock models in particular, and hence the question of whether the second law 
is valid in such cases is important and not yet solved (except for f{R) models, see [70]). The 
reference [70] also gives an argument for second law for arbitrary diffeomorphism invariant 
theories, but the argument relies on (a rather strong) assumption that stationary comparison 
version of the first law implies the physical process version for quasi-stationary processes. 

We shall now sketch a recent proof of the validity of the second law for quasi-stationary 
physical processes involving black holes in Lanczos-Lovelock gravity [71] . Such a result was 
first proved for the Einstein-Gauss-Bonnet theory in [72], but since the result in [72] is a 
special case of the one proved in [71] , we shall essentially outline the steps given in the latter 
reference. The main result of the analysis is that, Wald entropy of stationary black holes in 
Lanczos-Lovelock theories increases monotonically for physical processes which are quasi- 
stationary, at the end of which the black hole settles down to a quasi-stationary state. One 
of the reasons why this result is important is that it reinforces the physical fact that Wald 
entropy indeed behaves like ordinary thermodynamic entropy. 

We shall concentrate on stationary, non-extremal. Killing horizons. In D spacctime di- 
mensions, the event horizon is a null hypersurface Ji with tangent vectors k = d\, where 
A is an affine parameter and k satisfies geodesic equation. The constant A slices foliate the 
horizon and are spacelike. An arbitrary point p on such slices can be assigned coordinates 
{A, x^} where x^, {A = 2, - ■ ■ ,D) are coordinates of the point on A = shce connected to p 
by a horizon generator. To complete the basis, one needs to introduce a second null vector 
I such that I • k = —1, thereby getting {fc", e^}. Finally, we would require the induced 
metric on any slice, 7ab = gab + 2k(^ah), with k°-^ab = = l"'jab- The change in induced 
metric from one slice to another is given by [73] , 



where aab and 9 are shear and expansion, respectively, of the horizon. We will further assume 
that the event horizon is also a Killing horizon. This fact is true for general relativity due to 
strong rigidity theorem [74]. However, no such proof seems to exist for the Lanczos-Lovelock 
models to the best of our knowledge. With this assumption, the horizon is then generated by 
orbits of a Killing field ^ = dy which becomes null on the horizon, and surface gravity k is 
defined in the standard manner as ^°VaC'' = i^^^- For stationary spacetimes with a Killing 
horizon, aab = 0',9 = 0. Using evolution equation for 6 (Raychaudhuri equation) and shear, 
we obtain [73], 



Now consider a physical process wherein a weak stress-energy tensor perturbs a stationary 
black hole for certain duration, and after which the black hole settles to a stationary state 
in asymptotic future. This last point implies that ^ is an exact Killing vector in asymptotic 
future. Further, the flow of stress-energy is assumed to be slow, so that changes in dynamical 
fields are first order in some suitable (bookeeping) parameter e and viscous effects can be 
neglected. To be more specific, 6 ~ C(e), (Tab ^ C>{e). 




(158) 



= 0. 



(159) 
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In general relativity, a concrete example of such a physical process is slow accretion of 
matter by a black hole of mass M for a finite time, after which the black hole settles down 
to a stationary state. Linearized Raychaudhuri equation then gives, 

do 

— ^-Rabk''k''^-&T: Tabk''k\ (160) 

where Einstein's equation have been used to obtain the second equality. If Tab satisfies 
null energy condition Tabk°'k^ > 0, the rate of change of 6 is negative on any slice prior to 
asymptotic future. Since = in the future, the generators must have positive expansion 
during the accretion process. This suffices to prove the monotonic increase of area in the 
physical process. It is evident that the result depends crucially on field equations, and hence 
monotonicity of horizon area is only valid in case of general relativity. 

For Lanczos-Lovelock gravity, the starting point, of course, is the Wald entropy of a 
stationary Killing horizon [42, 75, 43], 

S=\JpVldA (161) 

where the entropy density (after separating out the m — 1 term) is given by (see Sec. 3.1): 

[D-l)/2] \ 

1 + J2 ^GT^rnarnL^^-J^^ . (162) 

m=2 / 

The aim is to prove that this entropy always increases when a black hole is perturbed by a 
weak matter stress energy tensor oiO{e) obeying null energy condition. To study the entropy 
change 

AS = \J^^^^+Op)dXV^dA. (163) 
it is useful to define a quantity Q as, 

e=(^ + <-.). 

In general relativity, is the same as the expansion parameter 9 of null generators, but for 
a general Lanczos-Lovelock model, & is the rate of change of the entropy associated with 
a infinitesimal portion of horizon ([70] discusses similar construction in f{R) gravity). Our 
aim is to prove that, given our assumptions, Q is positive on any slice in a physical process. 

To proceed further, we use the fact already noted in Sec. 4.1.2, that the change of the 
(£> — 2)-dimensional scalar L^^^^l can be thought of due to the change in the intrinsic metric, 



(™-i) 
I simpl 

th order Lanczos-Lovelock term in [D — 2) dimensions. So, for a general Lanczos-Lovelock 



and this variation of simply gives the corresponding equations of motion of (m — 1)- 
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gravity, we can write 



dp 
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in=2 
[D-l)/2] 

- E 
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(165) 



where the surface term has been ignored since horizon sections are assumed to be compact 
surfaces without boundaries. The tensor ^^~^''TZab, which is the generahzation of Ricci tensor 
for (to — l)-Lanczos-Lovelock term, is given by 
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Using Eq. (158), one now obtains. 
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To study the rate of change of along the congruence, we use the evolution equations 
for expansion (the Raychaudhuri equation) and shear [73]. Since we work to first order in 
perturbation over a background stationary spacetime, we can always extract the part linear 
in perturbation from a product of quantities XY as 



(168) 



where X'^^^ is the value of X evaluated on the stationary background, and X'^^^ is the 
linear part of its perturbation. Note that, for stationary background, Raychaudhuri equation 
demands ui^h'^k'' = 0. Since T^f ^fc'^fc'' = 0, this implies E^'^\ ^k^k^ = 0. We shall further 
simplify the calculation by using diffeomorphism freedom to make the null geodesic generators 
of the event horizon of the perturbed black hole coincide with the null geodesic generators of 
the stationary black hole [76]. 

Using the evolution equation of 9 and Gab to linear order: d6/dX f 
doab/dX ~ C^cJ^i^k'^k'^, and Eq. (159) for the background, the evolution of Q to linear order 
in perturbation can be written as 
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where 



Vabk'^k'' 



[D-l)/2] 

E 

m=2 



167ra„ 



{m)ao 



2TO(^-2)^(S)a6 (P) 
(m— 1) acoa 



(170) 



and we have used the expression of perturbed Weyl tensor in terms of curvature and Ricci 
tensors, along with the relation '^^'^^ E^^-^ab = '-^~'^^T^{rn)ab - (l/2)7ab^^~^'i(m)- 
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The next, and the most crucial, step is to prove that the first order part of T>abk°'k^ 
vanishes identically. We omit the proof of this, which can be found in [71]. It is, however, 
worth mentioning that the proof depends crucially on various algebraic properties of Lanczos- 
Lovelock models that have been discussed in earlier parts of this review. Whether there exists 
another way of proving the result which is adaptable to gravity theories more general than 
Lanczos-Lovclock remains an open issue. 

Once the first order part of Dabk^'k^ is shown to vanish identically, we finally arrive at, 

^ = -8^T,6fc"fc^ + 0(e2). (171) 
ctA 

The above equation shows that if null energy condition holds, the rate of change of Q is 
always negative during a slow classical dynamical process which perturbs the black hole to 
a new stationary state. Since, the final state by assumption is stationary, both 9 and u, and 
as a consequence 0, vanishes in the asymptotic future. We can use therefore use the same 
argument as in the case of general relativity to conclude that Q must be positive at every 
slice during the process, thereby proving that the horizon entropy of black holes in Lanczos- 
Lovclock gravity is indeed a monotonically increasing function during any quasi-stationary 
physical process, i.e. 

dS 

- > 0. (172) 

Further comments on the implications of this result, and suggestion for further investigations, 
can be found in [711. 



3.5 Entropy from horizon Virasoro algebra 

There is now considerable evidence to show that features like entropy, temperature etc — 
which were originally thought of as properties of black holes [37, 38] — arises in a much more 
general context. Several investigations suggest that there exist certain level of universality in 
the thermodynamic properties of null surfaces. Hence approaches which compute the entropy 
density of a null surface in a generic, universal, manner might provide us with deeper insight 
into the quantum structure of spacetimc. This aspect was emphasized, among others, by 
Carlip who has developed such an approach to compute the entropy using horizon Virasoro 
algebra. 

The basic idea is based on an approach first presented by Brown and Henneaux [77] , which 
was originally used in the context of asymptotically AdS space-time in (1 -I- 2) dimensional 
gravity. They found that the Fourier modes Qj of the charges corresponding to the asymptotic 
diffeomorphism symmetry generators obey a Virasoro algebra with central extension: 

C 

«[Qrn, Qn] = (m - ri)(9,„+„ + Y^m^S„,+n,o , (173) 

where C is known as the central charge. The work by Strominger [78] and others showed 
that if one uses the above central charge in the Cardy formula [79, 80], the resulting entropy 
comes out to be the Bekenstein-Hawking entropy for the (1-1-2) dimensional black hole. 
This was further developed by Carlip [81, 82] using the diffeomorphism symmetry generators 
near the horizon to obtain the black hole entropy. In this approach, one begins with the 
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diffeomorphism generators ^n{x) which preserve a set of boundary conditions near the horizon. 
The Fourier modes of these generators obey one sub-algebra isomorphic to Diff given by 

^{e™,^nr-("i-n)C+„ (174) 

where {, } is the Lie bracket. One can then construct the Fourier modes Q„ of the charges 
corresponding to each either by Hamiltonian [83] or a covariant Lagrangian formahsm 
[84, 41, 85, 86, 87, 88] and evaluate the Lie brackets among them. A comparison between 
this algebra and Eq. (173) allows us to identify the central charge. Finally, one finds the zero 
mode eigenvalue Qq and computes the entropy of the black hole using Cardy formula. 

Several related approaches have been developed using these ideas with the hope that 
diffeomorphism symmetry generators may shed some light towards the microscopic degrees 
of freedom responsible for entropy of the horizon [89, 90, 91, 92, 93, 94, 95, 96, 97, 98, 99, 100, 
101, 102, 103]. All these approaches developed in the literature have the following ingredients: 

• The Noether current used in the approaches are defined on-shell, usually by ignoring a 
term which vanishes when equations of motion are used. 

• The calculation of Lie brackets is complicated and different approaches lead to slightly 
different results and it is not clear how to interpret these differences in the calculations. 

• To obtain the correct result, one has to further impose specific boundary conditions on 
the horizon in order to set certain terms to zero. Again, the physical meaning of these 
boundary conditions is often not clear. 

• Most of the analysis (except the one in [99]) is confined to general relativity and it is 
not clear how to generalize the results for a wider class of theories. 

We will now describe an approach to horizon entropy based on the existence of a Virasoro 
algebra and its central charge as developed in Ref. [104] in which it was shown that there 
is a relatively simple way of obtaining the central charge and the horizon entropy using the 
off-shell Noether current in any generally covariant theory of gravity, and in particularly for 
the Lanczos-Lovelock model. The method does not require any boundary conditions to 
make unwanted terms to vanish. The essential idea is to use diffeomorphism invariancc of the 
Lagrangian under a;* — > a;* to define the Noether current Ja[Ci] ^'Hd then use its variation 

Ja[^i] under a second diffeomorphism a;' — > x'^ + ^l to define the Lie bracket structure. This 
can be done without using the explicit form of Noether current, and the calculation of the 
resulting Lie bracket is algebraically quite simple and leads to the standard results. One can 
then identify the resulting Virasoro algebra, the central charge and the zero-mode eigenvalue; 
plugging these into the Cardy formula, one obtains entropy density of a null surface which 
turns out to be the same as Wald entropy for the Lanczos-Lovelock models. 

The key new features of this approach are as follows: 

• The current is now defined and conserved off-shell i.e. we do not use the equations of 
motion in its definition and in its conservation. A derivation of Noether current and 
its conservation for on-shell condition is given in [105]. In this case a specific boundary 
condition is required to obtain the necessary results. But the approach in Ref. [104] 
deals with off-shell situation and does not require any ad-hoc boundary condition. This 
is in contrast to — and an improvement on — the previous approaches which use 
equations of motion extensively. 
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The definition of the bracket among the charges is completely general in the sense that 
one does not need the explicit expression for the Noether current to obtain its structure. 

This approach also provides a derivation of the central term (and zero mode energy) 
from the surface term contribution Jw" in the Noether current. This is new and does 
not seem to have been noticed in the previous works. 

• In the derivation of bracket we do not use the equations of motion or specific boundary 
conditions (like Dirichlct or the Neumann boundary condition) to make non-covariant 
terms vanish; therefore this approach is general enough to include any covariant theory 
of gravity. (This feature was again missing in the earlier literature which used the 
equations of motion for the evaluation of the relevant brackets as an on-shcU construct.) 
As a direct consequence of the above, the central term obtained is also an off-shell 
construct. 

Consider any generally covariant Lagrangian, for which the conserved Noether current J° 
(with a corresponding current density P°) that can be expressed as the covariant derivative 
of an anti-symmetric Noether potential J°^. These satisfy the standard conservation laws 
which arc valid ofF-shell: 

J" = '^bJ'''', P° = V5^"; VaJ" = 0, daP^^O. (175) 

Let us now consider the variation of the current density itself for an arbitrary diffeomorphism 
x° — >■ x° -I- and use it to obtain the variation of the Noether charge by straight forward 
calculation: 

S^Q = J d^aS^P" = J dY^ab^i^'r , (176) 

where h is the determinant of the induced metric of [d — 2)-dimensional boundary in a d- 
dimensional spacetime. We next define a (Lie) bracket structure for the charges which leads 
to the usual Virasoro algebra with the central extension for a general Lagrangian of the kind 
L = L(gab, Rated)- The relevant bracket among the charges is defined as: 

[Qi,Q2] = (%QK2]-%QKi]) (177) 

Then using Eq. (176), we obtain: 

[Ql, Q2] I VhdJ^ab [^U'i^l] - eiJ'[^2]] = I VhdJ^ab [^Ul - eiJ2] (178) 

where we use the notation = •^''Ki] etc. We note that: (a) This definition is quite 
general and has not used any field equations, (b) The derivation of Eq. (178) is simple and 
straightforward, (c) If we use the form J° = Vf,(V°^'' — V''^'^) for general relativity one 
can easily obtain the result obtained earlier in the literature, like e.g., in Ref. [82]. We will 
hereafter concentrate on Lanczos-Lovelock gravity and evaluate Eq. (178) on the {d — 2)- 
dimensional null surface which is a Killing horizon. 

To evaluate Eq. (178) over the Killing horizon, we will follow the stretched horizon ap- 
proach of Carlip [82] . Let us first mention some of the key results needed for this computation. 
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The location of the horizon is defined by the vanishing of the norm of a timelike (approx- 
imate) Kilhng vector x°' ■ Near the horizon, one can define a vector orthogonal to the 
orbits of the Killing vector x° , by the following relation 



(179) 



where k is the surface gravity at the horizon, with x°'Pa = 0. Consider a class of diffcomor- 
phism generators given by: 



(180) 



where T and R are scalar functions chosen such that the generators obey the (near-horizon) 
condition [{x'^x^)/x'^]3i9ab which preserves the horizon structure. This condition leads 
to a relation among R and T given by: 



X 



R^i^] x'^^aT = ^DT , 



(181) 



where D = x"Va- Wc can now compute the Noethcr charge: 



Qm = - 
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and the central term of the Virasoro algebra defined by the relation, 

K[Ci,^2] = [Qi,Q2]-Q[Ui,6}] , 

where [Qi, Q2] is given by Eq. (177). We find that the central term is given by 
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32nG 



^d''-^xp-'"■■''^lab^icd- 



DTi D^T2 - DT2 D'^Ti 



(182) 



(183) 



(184) 



This was obtained earlier in Ref. [99] by a more complicated procedure and using the on-shell 
expressions. Wc have obtained it without using the field equations or any special boundary 
conditions thereby demonstrating the generality of the result. 

Using a suitably defined Fourier decomposition of the Ti and T2 in Eq. (184) wc can find 
the central charge and zero mode eigenvalue. We use a Fourier decomposition of Ti and T2 
given by: 



Tl — AmTm', T2 — BnTn , 



(185) 



with A* = A_„, B*, = B^rn and 



1 



T„i = — exp [im {at + g{x) + p.x±) 



(186) 



where a is a constant, g{x) is a function that is regular at the Killing horizon, p is an integer 
and x± are the (d — 2) tangential coordinates. Here the t — x plane defines the null surface. 
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Obviously, the result will depend on the choice made for a and we need to fix this to get a 
unique value for entropy. A natural choice, arising from the fact that near-horizon Rindler 
geometry exhibits periodicity in imaginary time with period 27r/K, is 

a^K (187) 

Substituting Eq. (186) in Eq. (184) and defining a quantity 

i = -i y Vhd'''^xp-'''^^lab^^,d , (188) 

which is proportional to the Wald entropy, we obtain, 

A 

K[U,Cn] ,5„+„,o . (189) 

Similarly, using the Fourier decomposition, Q[^] = A„iQ[^„i] in Eq. (182), we obtain, 

m 

Q[U] - ^-^".o . (190) 
and the algebra of charges turns out to be: 

A 

i[Qmi Qn] = (m - n)Q[^,„+„] + -——5n+,n,0 ■ (191) 

OTTLr 

This is the standard form of the Virasoro algebra in Eq. (173) with Q[£^rn+n] = Qm+n- We 
can identify the central charge and the zero mode eigenvalue as: 

S = 8^= ^[^"] = 8^ • (19^) 
The standard Cardy formula for the entropy is given by [79, 80], 



S^2^\j—; A^Qo-^ (193) 



which leads to 



(194) 

This exactly matches with the Wald entropy of the theory. 

The conceptual basis of this approach was concretely demonstrated (in the case of general 
relativity) from another point of view [106]. One can show that the Noether current that arises 
from the surface term of the action in general relativity, for diffeomorphisms which preserve 
the horizon structure, also possess a Virasoro algebra with similar properties. More precisely 
ref. [106] considers a class of diffeomorphisms which leave a particular form of the Rindler 
metric invariant and introduces the Noether current associated with the surface term in the 
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action functional for gravity (in contrast to previous approaches in the literature in which 
either the bulk term or the Einstein- Hilbert action was used to define the Noether current). 
Given the Noether current and diffeomorphism generators, there is a natural Virasoro algebra 
which can be associated with the horizon. The central extension of this Virasoro algebra again 
leads — via Cardy formula — to the correct entropy density of the horizon. Just to provide a 
contrast with the previous analysis, we will mention couple of technical points. In the earlier 
work, to obtain the correct entropy one had to either shift the zero mode energy (as done in. 
e.g., [82]) or choose a parameter contained in the Fourier modes of T (for instance a) as the 
surface gravity (as done in, e.g., [103]) or both [104] depending on the action of the theory we 
have started with. Here, interestingly, one did not require any such ad hoc prescription. This 
is because the parameter a in the expression of does not appear in the final expression 
of entropy in this approach. 

Using (a) the Noether current corresponding to the York-Gibbons-Hawking surface term 
and (b) choosing the vector fields by demanding invariance of the Rindler form of the metric, 
have not been attempted before and are new features of the work in Ref. [106]. This is 
important because it further strengthens the peculiar feature present in all gravitational 
actions, viz., that the same information is encoded in both bulk and boundary terms of the 
action. 

From a conceptual point of view, the analysis in Ref. [106] provides a nice, simple, phys- 
ical picture of the connection between horizon entropy and the degrees of freedom which 
contribute to it. In conventional physics, we are accustomed to thinking of degrees of free- 
dom of a system as absolute and independent of observer or the coordinate system used by 
the observer. In such a picture, the entropy — which is related to the logarithm of the 
degrees of freedom — will also be absolute and independent of the observer. We however 
know that horizon entropy, horizon temperature etc. must be treated as observer dependent 
notions because a freely falling observer through a black hole horizon and a static observer 
outside the black hole will attribute diff'erent thermodynamic properties to the horizon. It 
follows that any microscopic degrees of freedom which leads to horizon entropy must also be 
necessarily observer and coordinate dependent. So the question "what arc the degrees of 
freedom responsible for black hole entropy?" has no observer independent answer. We find 
that the notion of residual gauge symmetries which are perceived as real degrees of freedom 
under a restricted class of diffeomorphism, allows us to quantify this notion. An observer 
who perceives a horizon works with a theory having less diffeomorphism symmetry if she 
wants to retain the structure of the metric near the horizon. This necessarily upgrades some 
of the original gauge degrees of freedom to effectively true degrees of freedom as far as this 
particular observer is concerned. As a result, this particular class of observers will attribute 
an entropy to the horizon in an observer dependent manner. 

The result of [106] has already been generalised to m = 2 and m = 3 Lanczos-Lovelock 
terms in [107], using the corresponding surface terms for these lagrangians given in Sec. 2.6, 
and we expect it to hold for all Lanczos-Lovelock terms in a straightforward manner [108]. 
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4 Lanczos-Lovelock models and the emergent gravity 
paradigm 

Recent results indicate that gravity could be an emergent phenomenon, like fluid mechanics 
or elasticity, with the field equations governing gravitational dynamics having the same status 
as, say, the equations of fluid mechanics. This alternative perspective for classical gravity 
appears to be conceptually more satisfying than the usual perspective in certain aspects. (See 
for a review, Rcfs. [3, 4], [109].) and several features of Lanczos-Lovelock models strengthen 
this paradigm. In this section we shall describe these features. 

The back drop in which such an approach arises is the following: it is reasonable to 
assume that there are certain pre-geometric variables underlying the spacetime structure and 
that there exists a microscopic theory describing their dynamics. If we think of continuum 
spacetime as analogous to a continuum description of a fluid, the microscopic — quantum 
gravitational — description is analogous to statistical mechanics of the molecules of the 
fluid. (That is, the pre-geometric variables will act as the "atoms of spacetime".) The exact 
quantum theory of gravity, when it is available, will allow us to construct, in a coarse-grained 
long wavelength limit, a smooth spacetime and the effective degrees of freedom (like the 
metric tensor) in terms of pre-geometric variables. This procedure is analogous to obtaining 
variables like pressure, temperature etc. for a fluid in terms of microscopic variables. The 
dynamical equations of the microscopic theory will also lead to some effective description 
for the emergent degrees of freedom. For example, given the microscopic description of 
molecules moving randomly inside a container and colliding with the walls, one can obtain the 
ideal gas law PV = NksT governing the macroscopic ("emergent") variables like pressure, 
temperature etc. Similarly, once we have the correct theory of quantum gravity (in terms 
of pre-geometric variables), we should be able to obtain the field equations governing the 
geometry by a suitable approximate description. 

Such an approach, however, works only when we know the underlying microscopic theory. 
How do we proceed when we do not know the underlying theory, which is the current situ- 
ation as regards the spacetime? In the case of normal matter, this was indeed the situation 
couple of centuries back and physicists could successfully describe the behaviour of matter 
using thermodynamic concepts, even before we understood the molecular structure of matter. 
The emergent paradigm for gravity can be viewed as a similar attempt to describe gravity 
in a thermodynamic language, given the fact that we do not know the exact, microscopic 
description of spacetime. 

In the case of normal matter, the evidence for the existence of microscopic degrees of 
freedom is provided by the elementary observation that matter could be heated up. Boltz- 
mann emphasized the fact that matter could store and transfer energy in the form of heat 
only because it had microscopic degrees of freedom. If a fluid admits continuum description 
all the way and has no internal degrees of freedom, it cannot store heat or have a notion 
of temperature. Thus the existence of thermal phenomena provides strong internal evidence 
that continuum fluid mechanics misses some essential aspect of physics. This is clearly seen, 
in the case of an ideal gas, if we write 

PV 

NkB = — (195) 
The quantities on the right hand side of Eq. (195) are well defined and meaningful in ther- 
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modynamics but the N in the left hand side has no physical meaning in thermodynamics. It 
is proportional to the number of microscopic degrees of freedom of the gas and is irrelevant 
in the strictly continuum limit. With hindsight, we now know that the existence of the gas 
constant R = Nks for a mole of gas has a microscopic interpretation related to the internal 
degrees of freedom. This is a situation in which we could have guessed (as Boltzmann did) 
the existence of microscopic degrees of freedom — and even counted it, as Avogadro number 
for a mole of gas — from the existence of thermal phenomena in normal matter. 

In a similar vein, we can ask whether the continuum description of spacetime dynamics 
exhibits evidence for a more fundamental level of microscopic description in terms of "atoms 
of spacetime" . It would, of course, be impossible to describe the statistical mechanics of such 
microscopic degrees of freedom knowing only the continuum description of spacetime, just as 
one could not understand the constituents of matter in the days of Boltzmann. But it must 
be possible to provide internal evidence strongly suggesting the existence of such atoms of 
spacetime just by probing the continuum theories of gravity. 

The first task of the emergent paradigm will be to come up with such 'internal evidence' 
by closely examining classical theories of gravity. We will discuss several peculiar features 
in the structure of classical gravitational theories which have no deeper explanation within 
the standard framework and have to be accepted as just algebraic accidents. These peculiar 
features provide us with hints about the underlying microscopic theory. In that sense, they are 
conceptually similar to the equality of inertial and gravitational mass (which could have been 
thought of as an algebraic accident in Newtonian gravity but does find a deeper explanation 
when gravity is described as spacetime geometry) or the fact that matter can be heated 
up (which defied a fundamental explanation until Boltzmann postulated the existence of 
microscopic degrees of freedom) . These features strongly suggest interpreting classical gravity 
as an emergent phenomenon with its field equations having the same status as equations of, 
say, fiuid mechanics. Careful analysis of classical gravity (with a single quantum mechanical 
input, viz., the Davies-Unruh temperature [39] of local Rindler horizon) leads us to this 
conclusion. 

In the case of normal matter, a more formal link between microscopic and macroscopic 
description is established by specifying certain thermodynamic potentials like entropy, free- 
energy etc. as suitable functions of emergent variables like pressure, temperature etc. The 
functional form of these potentials can, in principle, be derived from the microscopic theory; 
but, when we do not know the microscopic theory, they are postulated phenomenologically 
from the known behaviour of the macroscopic systems. Similarly, we would expect such a 
thermodynamic approach to work in this case of spacetime as well. It should be possible to 
express, say, an entropy or free energy density for spacetime in terms of suitable variables, the 
extremum of which should lead to a consistency condition on the background spacetime — 
which will act as the equation of motion. Such a consistency condition arises even in normal 
thermodynamics, though it is not often described as one such. For a gas of N molecules in 
a volume V, we can express the kinetic energy and momentum transfer through collisions 
to the walls of the container per unit area and time, entirely in terms of the microscopic 
variables. Coarse graining these quantities we obtain the macroscopic variables T and P. 
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The equations of microscopic physics now demand the consistency condition 



mean momentum transfer to walls 
per unit time per unit area 



N 



P 



(196) 



(mean kinetic energy) 



V 



T 



between the coarse grained variables. We will show in Sec. 4.2 that such an approach is 
indeed possible to obtain the field equations of gravity. 

It is also true that, in the case of normal matter, the laws of thermodynamics hold uni- 
versally and do not depend on the kind of matter (ideal gas, liquid crystal, metal, ... ) one 
is studying. The information about the specific kind of matter which one is studying is pro- 
vided by the specific functional form of the thermodynamic potential lik e.g., the free energy 
F = F(T, V). Similarly, the thermodynamic framework can describe a wide class of possible 
gravitational field equations for the effective degrees of freedom. Which of these actually 
describe nature depends on the specific functional form of the thermodynamic potential, say, 
the entropy density of spacetime. It turns out that, for the Lanczos-Lovelock models which 
we will concentrate on, this information is encoded in the tensor P^^ (which we will call the 
entropy tensor) and the nature of the resulting theory depends on the dimension of spacetime. 
In particular, if Z? = 4, the thermodynamic paradigm selects Einstein's theory uniquely. 

This approach can be thought of as a "top-down" view (in real space) from classical gravity 
to quantum gravity. Specific quantum gravitational models which approach the problem 
'bottom-up' has to be consistent with the features of classical gravity described in the sequel. 
In particular, this approach makes precise the task of the microscopic quantum gravity model 
viz., that it should lead to a specific functional form for the entropy or free energy density 
of spacetime, just as microscopic statistical mechanics will lead to a specific entropy (or free 
energy) functional for a material system. In this sense, the emergent approach complements 
microscopic approaches based on models of quantum gravity. Such a perspective also implies 
that quantizing any classical gravitational field will be similar to quantizing, say the equations 
of fluid dynamics. Gravitons will have the same status as phonons in a solid. Neither will give 
us insights into the deeper microstructure — spacetime or atoms. (The alternative paradigm 
also has important implications for cosmology and, in fact, suggests linking cosmology with 
the emergence of space itself in a special manner. The details can be found in ref. [109]. ) 

4.1 Emergent aspects of the Lanczos-Lovelock equations of motion 

The first 'internal evidence' suggesting that Lanczos-Lovelock models describe an emergent 
structure of spacetime arises from the remarkable fact that Lanczos-Lovelock field equations 
reduce to a thermodynamics identity TdS = dE+PdV on the horizon in any static spacetime. 
We will begin by describing this result. Obviously, it will depend crucially on the near 
horizon properties of certain tensors associated with Lanczos-Lovelock field equations which 
we will describe in Sec. 4.1.1. These will be shown to lead to the thermodynamics identity 
TdS = dE + PdV on the horizon, in Sec. 4.1.2. 

4.1.1 Near-horizon symmetries of equation of motion 

To investigate the near-horizon symmetries, it is convenient to use a setup which is best 
suited to describe geometry in the neighbourhood of a static horizon. We will describe the 
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static case (viz., the existence of a timelike Killing vector at least in the neighbourhood of 
the horizon); the more general case of stationary horizons will only be mentioned when we 
require it. 

A coordinate system which facilitates taking the near horizon limit is given by the Gaus- 
sian normal coordinates, which, for static spacetimes, simplifies further and can be written 
as [110] 

ds^ = -N^dt^ + An^ + UABdy'^dy^ (197) 

where (TAB{n, y^) is the transverse metric, and the Killing horizon, generated by the timelike 
Killing vector field ^ = dt, is approached as — >■ 0,n — 0. We give below some geometric 
quantities associated with the above metric, which will be put to good use in subsequent 
analysis. 

1. Taylor expansion of metric components: Demanding the finiteness of curvature invari- 
ants on the horizon leads to the following Taylor series expansion for the lapse N{n, y) 
and the transverse metric aAB [HO]: 



N(n,y) = nn 



lAB = ['TH(2/)]^B + -[a2(2/)]^s"' + 0(n') (198) 



where we have collectively called the transverse coordinates as y, and R± is the Ricci 
scalar corresponding to y = constant surface. 

2. Domain of validity: We shall be interested in a small region in the neighborhood of 
the spacelike {D — 2)-surface; more precisely, we shall assume n <^ ' , and would 
be interested in the k — >■ oo. That is, we shall require the length scale set by k to be 
the smallest of all length scales in the problem. Note that the coordinate n gives the 
normal distance from the horizon. 

3. The null geodesic congruences: Consider the null vectors given by: 

= (-1,+A^~\0,0) 
h = {-l,-N-\0,0) (199) 

Explicit calculation shows that VjZ has only "y" -components, given by: 

Vil = ~ {2k-^) [a^^ dBR±] 9a + 0{n^) (200) 

which go to zero in the limit dR±/ — 0; therefore, in the limit we are interested in, it 
is easy to check that these vectors satisfy the geodesic equation in affinely parametrized 
form, that is; Vil = = Vfcfc. 

The affine parameter A, defined by I ■ VA ^ 1 can be found by using the above form of 
N{n,y); to the leading order, we find that, A ~ A^ -I- (l/2)Kn^ where A = \h is the 
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location of the horizon. Note that, N'^l — ;> ^\h, which imphes, 2k (A — A^f) I — > ^\h- In 
some part of the subsequent analysis, the differentials of various geometric quantities 
(such as entropy) defined on the horizon, which are directly involved in the statement 
of the first law of thermodynamics, would be interpreted as variations with respect to 
the affine parameter along the outgoing null geodesies, i.e., A. All of above comments 
hold (except for few sign changes) also for the ingoing null vector k. 

Comoving observers and their acceleration: Let us next look at the comoving observers 
in this spacetime, n =const, y =const, with 4-velocity u = N^^dt- It is easy to check 
that their acceleration a = VuM is given by: 

" = -^a'-'^^uN (201) 

The components can be worked out explicitly by using Taylor series expansions in 
Eqs. (198); these are given by [with C = -{l/2)R±{y- n = 0)] 



a± 

„2 



A R 

a aAas 



+ AC + 0{n'^) + 0{n^) 



so that 



Therefore, we see that, 



— = 1 



0(n3) ; ^ 
a 



a = [1 + 0{n^)] dn + 0{n^)dA 



(202) 



(203) 



This result indicates how the curvature i?^ of the t ^ n part of spacetime, which is a 
plane spanned by u and n, will in general tend to produce an acceleration in the "y" 
directions. We see that even in a curved spacetime, the unit vector along a points along 
n to an excellent approximation [0(n'^)]. 

Getting back to the symmetries of equations of motion, we can now use the above near- 
horizon form of the metric to show that the Lanczos-Lovelock equations of motion satisfy 
the following near-horizon symmetries, generalizing the result first derived in [110]: 
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(204) 



where E± is given by 

Ei_ 

cB 



0-2 



1 m 
167r2™ 

j-BAi...Bm- 

°ACi...D„,- 



CB 



^CAcB 
a t A 



1 



1 (-D-2) nCiDi (D-2) n' 



(205) 
(206) 



59 



In Sec. 4.1.2, we will use the above form when proving that Lanczos-Lovelock field equations 
reduce to TdS = dE + PdV. Here, we give a brief sketch of how the above, highly symmetric, 
form of Lanczos-Lovelock equations of motion arise near a static horizon as a result of the 
geometric structure of the Lanczos-Lovelock Lagrangian. More details can be found in the 
appendices of Ref. [111]. 

Let us first briefly outline the proof for Einstein tensor. We begin with the following 
expressions for the decomposition of the Ricmann tensor for the metric in Ec^. (197) (sec, for 
example. Section 21.5 and Exercise 21.9 of [26]). 





= 








— (3) DM 



(207) 
and 

''^^Rabcd ~ ^^^Rabcd — {KacKbd — KadKbc) 

^^^RnBCD = Kac-.b - Kab-.c (208) 

Here, | and : arc the covariant derivatives compatible with the induced metric on {t = 
constant} and {t — constant, n = constant} surfaces respectively, and Kab = — (l/2)9nO'AB 
is the extrinsic curvature of the {t = constant, n = constant} 2-surface as embedded in the 
{t = constant} surface. These expressions lead to 



jytA AC 

R tB = 



BCi 

N,CB - KcBdnN 



N 

Using the relevant Taylor series expansions (see Eqs. (198)), we obtain 



(209) 



5„trif |„=o = -^tr[(T2] (210) 



which gives, correct to O(n^), 



and 



which, in turn, implies 



R"^A^~7:tr[^2]^R'fA (211) 



i?^g = (^)i?^g + 0(n2) 



Ri^ = ii-ll + 0{n') (212) 



Finally, we use the following general expression for the Einstein tensor (see, for example, 
[26], section 14.2, pp. 344): 



— ^ \ ^ iiA + 7;^AB 



1 

2' 

g: = -{R'tA + UiB] (213) 
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Plugging in the above expressions for the Riemann tensor, and finally taking the limit n ^ 0, 

we immediately obtain part of Eq. (204). (Note that, with the up-down components, G~ = 

Gj.) One can go further and analyze, in the same way, the remaining components of the 
Einstein tensor; the final result is [110] 
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where Ef\H 
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gi\h. 



One can now proceed to do a similar analysis and obtain the near-horizon symmetries for a 
generic Lanczos-Lovelock Lagrangian. The full result, turns out to be the same as Eq. (204). 
We skip the proof which involves a bit of combinatorics, but otherwise straightforward to 
establish. Here, we shall only give a proof of the identity i?| = E'^ which would be 
directly relevant when we discuss thermodynamic structure of the field equations. 

Using the first equality in Eq. (30), it is easy to deduce the form of E^ (no summation 
over k). The alternating determinant simplifies upon using 5^ — 1, and the fact that none of 
the other indices can be k due to total anti-symmetry. Therefore, we are left with 



El = -\L,^{k] 



(215) 



where Lm {k^ denotes terms which do not contain the index k at all. We now specialize to 
the case when the index k corresponds to n. The RHS can be further split depending on the 
number of occurrences of the index t, as 



En = -\ [E,n {n, t} + L„, {n, 2t} + L,n {n,t}] 



(216) 



The first set on the RHS contains terms like R'^^g which arc identically zero (see first of 
Eqs. (207)), while the last set is the same as appears in E^. So we only have to prove that 
Ljn {n, 2t} = L„i {t, 2n}. The set Lm {n, 2t} will have two t^s appearing either on different 
factors of which would again vanish identically for the same reason as the first set. or 

it can have the two Vs appearing on the same factor, which would contribute 
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(217) 



Using Eqs. (209), we see that L„ {n, 2t} = L,„ {t, 2n} + 0{n'^). Therefore, E^l = Ef + 0{n'^), 

and the equality holds on the horizon. The form for E^^^ maybe deduced in the same way, 
but since it is not of much use in study of horizon thermodynamics, we do not give that 
analysis here. 
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4.1.2 Thermodynamic structure of Lanczos-Lovelock field equations 

To demonstrate the equivalence of field equations with the identity TdS = dE + PdV, we 
shall use the results derived above in Sec. 4.1.1. We begin with the t — t component of the 
field equations, and, following an analysis similar to the one leading to Eq. (30) of [112], we 
arrive at, 



where 



4 = ^i-^^i-^^cB-'^^l-lLi--'^ + oH 



— "ACi...r>,„_i ^AiBi ' ' ' ^A,„-iB^-i 



(218) 



(219) 



In Sec. 4.1.1, we proved the equality E^ — E~ on the horizon. Therefore, we have 

p?| _ pn I _ ^ ""^ \„ 1 rrCAcB 1 t{D-2) 

Using again the Taylor series from Sec. 4.1.1, and the affine parameter A introduced in that 
section, we also have 

Sx<JAB = — W2]ab + 0[{\ - XhY'^ 5\] (220) 
Using this, we arrive at the expression: 

2e1^5X ^ ri^-^y^^ 5,<yscV^- L\^-'^V^5X (221) 



5a S 



where we have introduced the Hawking temperature, T = k/{2tt). As the final step, we 
will now show that the factor multiplying T, as well as the second factor on the right hand 
side, can both be expressed as variations of quantities defined in terms of intrinsic horizon 
geometry, as is anticipated by the notations "i5aS"' and ''d\E" we have used for them. This is 
a remarkable fact, and is most probably true only for the Lanczos-Lovelock models. Moreover, 
the quantities S and E will turn out to be the entropy and energy associated with the horizon 
in Lanczos-Lovelock theory. 

To begin with, consider the quantity 

S = 47r?7i J dS iJ^Ji^^ (222) 

and note that, purely algebraically, the variation of the above expression must give equations 
of motion for the (m — 1)"^ order Lanczos-Lovelock term in {D — 2) dimensions. (The 
variation would also produce surface terms, which would not contribute when evaluated at 
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A = Xh because we will assume that the horizon is a compact surface with no boundary.) 
We have: 



SxS ^ -47rm / dS 8'^^ SxObc (223) 

JH 

where we have evaluated the variation on A = \h ■ Noting that the Lagrangian is 

j-{D~2) _ 1 1 j;AiBi...B„,-i (D-2) pC„_iL>„-i 

and using the first of Eqs. (30), we see that 

1 : 



£'c = AiLt^I^^E (224) 



Therefore, we obtain 



^>^S = l^i^ 6,aBc (225) 



which is precisely the integral of the factor multiplying T in Eq. (221). As mentioned above, 
S defined in Eq. (222) is a function of A, and its derivative with respect to A is well defined 
and finite on the horizon. The expression for S, evaluated at A = A//, 

5'(A = Ah) =47rm / dE L^^Z^'^ (226) 

JH 

is what we shall interpret as the entropy of the horizon. As mentioned before, S has no 
physical significance away from the horizon. (As an aside, we would also like to point out the 
following fact; when D = 2m, that is, when the corresponding Lanczos-Lovelock action is the 
Euler characteristic of the full Z?-dimensional manifold, then D — 2 = 2(m — 1), and we see 
from Eq. (226) that 5 is proportional to the Euler characteristic of the {D — 2)-dimensional 
horizon surface, determined by L^^S-i' therefore, S is just a constant number. For example, 
the Gauss-Bonnet term (m = 2) in = 4 will contribute a constant additive term to the 
standard Bekenstein-Hawking entropy.) 

Muhiplying Eq. (221) by d'^'^^y, integrating over the horizon surface, and taking the 
n — >■ limit, we now see that it can be written as 

T^xS - / dE 4f -^A = / dS tI 6X 

dE r| dX 

H 



dE P_i ^A (227) 



H 



where we have used the field equations E~ ~ {1/2)T~ in the first equality, and the relation 
E~ \h = EZ \h in the second equality. We have also used the interpretation of T~ as normal 
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pressure, P±, on the horizon. The above equation now has the desired form of the first law 
of thermodynamics, provided: (i) We identify the quantity 5, defined by Eq. (226) as the 
entropy of horizons in Lanczos-Lovelock gravity; indeed, exactly the same expression for 
entropy has been obtained in the hterature using independent methods, see e.g, Ref.[75]. 
(That the same expression is obtained using the Noether charge approach of Wald, was 
demonstrated in Ref. [113].) (ii) We also identify the second term on left-hand-side as 6xE; 
this leads to the definition of E to be 

E = J SxJdT. Li^-^'^ (228) 

where the A — > Xh limit must be taken after the integral is done (therefore, we would require 
the detailed form of Lm~'^'' as a function of A to calculate this explicitly). For D = 2{m+ 1), 
the integral over H above is related to the Eulcr characteristic of the horizon, in which case 
E (X Xh (upon setting the arbitrary integration constant to zero). For m = 2,D = 4, this 
reduces, of course, to the expression obtained earlier in the case of Einstein gravity. 

Let us briefly comment on the general form of E for spherically symmetric spacetimes for 
general Lanczos-Lovelock Lagrangians, with horizon at r = rn, and A = r. In this case, 
im* ~ (l/A^)™ and also ^/a ~ X^~^. The integrand in the energy expression therefore 
scales as \iD-2)-2m ^ Integrating the RHS of Eq. (228), and taking the X ^ Xh hmit after 
doing the integration, we see that E ^ A^ As mentioned above, for D ~ 2{m+ 1), 

E ~ Xh- In fact, in the case of spherically symmetric spacetimes in Lanczos-Lovelock theory, 
the above expression can be formally shown to be exactly equivalent to the one already known 
in the literature (see [112], and also Ref. [16] therein). As far as we are aware, no general 
expression for energy in Lanczos-Lovelock theory exists in the literature, and the one we have 
arrived at could be thought of as one definition which appears to be reasonable from physical 
point of view. In any case, the expression for energy clearly deserves further investigation. 

Putting everything together, we see that, for generic static spacetimes in Lanczos-Lovelock 
gravity, the field equations can be written as a thermodynamic identity: 

TSxS - SxE = F6X (229) 

where we have found it appropriate to interpret 

[ P±V^d^y (230) 

as the the average normal force over the horizon "surface" (somewhat in line with the spirit 
of the so called membrane paradigm, which we discuss in further detail in Sec. 3.3). The term 
F SX has the interpretation of virtual work done in displacing the horizon by an afhne distance 
(5A, in perfect analogy with the work term in the standard first law of thermodynamics. 

We have therefore established the equivalence of Lanczos-Lovelock field equations with the 
thermodyanamic identity, generalizing the result already known in Einstein gravity. In fact, 
as will be a recurring point in this review, we have achieved more than just a generalization of 
result from Einstein to all of Lanczos-Lovelock - the analysis here provides a more rigorous 
and robust geometric basis for the thermodynamic interpretation of field equations in Einstein 
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theory itself. To put the result in appropriate context, note that we have shown that the 
following relations 

E{n,n) = -E{u,u) & E{u,u) ^ (1/2) T{u,u) 

near— horizon symmetry field equations 

have a thermodynamic structure. Democracy of all observers then implies that E = (1/2)T 
also has a thermodynamic structure: 

Eiu,u) = ^Tiu,u) ^ E^\t 

L all u Z 

Further, in the limit — > 0, we have N'^k — > Using this and u = ^/-/V, we see that 
T{u,u) — > r(^,fc), which is the flow of energy for a virtual displacement of the horizon 
along k. Another way to interpret this result is as follows: the object g^f^ = 2^(afcb) is 
the transverse part of the induced metric, and hence Tabi°'l^ = (1/2) tr2[Tah], where tr2 is 
the trace with respect to the metric 2^(Qfct,). For spherically symmetric spacetimes, since 
To° = T; = P, we have (1/2) traiTat] = P, which is the result derived in Rcf. [30, 114]. 



4.1.3 Field equations from the Clausius relation 

The above result suggests a strong relation between gravitational dynamics and horizon ther- 
modynamics. This idea can be further strengthened by another result, viz. that gravitational 
field equations can be interpreted as the entropy balance condition on local Rindler horizons. 
In the most straight forward manner, this result can be obtained very simply from the fol- 
lowing argument [115]. 

Consider an infinitesimal displacement of a local patch of the stretched (local Rindler) 
horizon T-L along the direction of its normal r^, by an infinitesimal proper distance e, which will 
change the proper volume by dVprop = ey/^d^~^x where aab is the metric in the transverse 
space. The flux of matter energy through the surface will be T^^^Va (where be the 
approximate Killing vector corresponding to translation in the local Rindler time) and the 
corresponding entropy flux can be obtained by multiplying the energy flux by /3ioc = N j3. 
Hence the 'loss' of matter entropy perceived by the outside observer, because the virtual 
displacement of the horizon has engulfed some matter, is 

SSra = PxocSE ^ PlocT'^'^arjdVprop. (231) 

Interpreting PiocJ"" , where J° is the Noether current corresponding to the local Killing vector 
given by J° = 2Q^^'' + i^", as the gravitational entropy current, the change in the 
gravitational entropy is given by 

= l^loc'^aJ dVprop (232) 

(Note the occurence of the local, rcdshifted, temperature through f3\oc — N (3 in these expres- 
sions.) As the stretched horizon approaches the true horizon, Nr°- — )• and f3^°-(,aL 0. 
Hence we get, in this limit: (JS'grav = P^aJ°'dVprop = '2,/3Q"'^ ^a^jdVprop- Comparing (55grav 
and SSm we see that the field equations 2Q^ = can be interpreted as the entropy balance 
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condition SSgrav = SSmatt- This provides a direct thermodynamic interpretation of the field 
equations as local entropy balance in local Rindler frame. 

Though we work with entropy density, the factor /? = 2t:/k cancels out in this analy- 
sis — as it should, since the local Rindler observer with a specific k was introduced only 
for interpretational convenience — and the relation T5Sm — TdSgi-^v would have the same 
interpretation. The expression in the right hand side is the change in the horizon ('heat') 
energy Hgur = TS due to injection of matter energy. (For a more detailed discussion of this 
expression, see Ref. [116]). The context we consider corresponds to treating the local Rindler 
horizon as a physical system (like a hot metal plate) at a given temperature and possessing 
certain intrinsic degrees of freedom. Then integration SS = SE/T at constant T to relate 
change in horizon energy to injected matter energy. Any energy injected onto a null surface 
hovers just outside the horizon for a very long time as far as the local Rindler observer is 
concerned and thermalizes at the temperature of the horizon if it is assumed to have been 
held fixed. This is a local version of the well known phenomenon that, the energy dropped 
into a Schwarzschild black hole horizon appears to hover just outside R = 2M as far as an 
outside observer is concerned. In the case of a local Rindler frame, similar effects will occur 
as long as the acceleration is sufficiently high; that is, if k/k^ ^ 1. 

While the above results hold quite rigorously, they only show that the field equations 
lead to an entropy balance condition while what we are interested in is the converse. We 
want to prove that the entropy balance condition implies the field equations. There has been 
attempts in the literature to 'reverse-engineer' the above argument [117] but these ideas face 
a major stumbling block. As pointed out by Padmanabhan [118] the reverse engineering 
requires interpreting J° as an entropy current which in turn is justified only if we assume the 
field equations (e.g, Wald needed the field equations to show that Noether charge is physical 
entropy.) 

Clearly the idea works only if we can have an expression for entropy which can be inde- 
pendently justified. In Einstein's theory, entropy is proportional to area of the horizon which 
allowed Jacobson [119] to use the entropy balance condition to obtain the field equation in 
general relativity from Clausius relation. The result relied on three crucial inputs: (a) The 
entropy of the horizon is proportional to the area, (b) Raichaudhury equation relates the 
change in the area of cross-section of null congruence to the matter stress tensor, (c) Certain 
technical assumptions regarding the shear and expansion of the congruence. Given these, it 
is straight forward to obtain Einstein's equations from the Clausius relation. 

Various attempts have been made towards generalization of the original result of Jacobson 
to higher curvature gravity theories. Since entropy is no longer proportional to the area in a 
general Lanczos-Lovelock model, the Raichaudhury equation is pretty much useless as a tool 
and there is no way of directly generalizing the original idea of Jacobson. As a consequence, 
all the "generalizations" arc cither quite different from the original idea at a conceptual level 
(in that the Clausius relation is not employed in the way it was originally done), or else 
marred by technical errors and/or ad hoc assumptions tailor-made to get the desired result. 
The most immediate generalization given by Jacobson and collaborators themselves [120], 
for the f{R) theories, is distinct in concept in that additional dissipative terms had to be 
added to the heat flux to get the "correct" field equations. (Obviously, any field equation 
can be written in the form TdS — dE = (something) with the right hand side interpreted 
as giving some 'dissipation", so this result is conceptually vacuous.) More recent attempts 
towards generalization are, on the other hand, based on a judicious use of the expression 
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for the Noether potential Jab\ essentially, since diffeomorphism invariance of the action is a 
non-dynamical symmetry, the corresponding current J° is conserved off-shell, and using its 
expression for a Killing field, one can hope to derive the equations of motion from a relation 
like TdS = 5E. In these works, while the definition of matter flux is similar to Jacobson's, 
the Raychaudhuri equation is never invoked, thereby avoiding any need for some assumptions 
such as vanishing expansion etc. present in original Jacobson work. 

The situation, however, is far from clear since these works do not all agree with each other. 
For e.g., Parikh and Sarkar have pointed out issues with the Brustein and Hadad [117] paper 
whereas [118] has highlighted certain conceptual issues regarding interpretation of the results. 
Specifically, [118] has stressed some of the subtleties (alluded to above) in interpreting these 
results as a derivation of field equations from thermodynamics; these results, as Ref. [118] 
argues, are better viewed as interpretation of field equations as a thermodynamic relation, 
rather than a derivation of the former from the latter. More importantly, Brustein and Hadad 
derive their result using different definitions for various quantities as compared to the other 
two papers; this calls for a more detailed and critical look at the analysis therein. 

The most recent work in this direction is that in ref. [121] which we shall briefly describe. 
First, one starts with the following (very general) ansatz for horizon entropy 

S = [ s'^^'eabdA, (233) 

(where the entropy density s"'' is an antisymmetric tensor, e°'^ is the binormal to the slice and 
dA is the area element on the slice), and then calculate the change of this entropy between 
two horizon slices S and Eo using Stokes' theorem 

6S ^ (f s^'^eabdA = -2 / Vf.s"'' kadVdA (234) 

where V is a parameter along k (which is tangent to horizon generators) . Various assumptions 
about the slices and enclosed horizon region H can be found in [121]. Second, the Clausius 
relation, dS ~ SQ/T, is imposed with the heat flux given by 

SQ^ [ {-Tab)^'dHa (235) 

Matching the order on two sides of the above relation appropriately using behaviour of the 
local Killing vector (LKV) ^-^ gives 

-(fi/7r)V,s''^ ka = T^'^a ka + Oix^). (236) 

where 0{x^) refers to higher order terms in expansion of the LKV. Note that, at this stage, the 
relation derived above is nothing but the local entropy balance condition which was obtained 
earlier in [115] , except that the form of Sab has not yet been fixed. This third step is attempted 
as follows. One motivates a form for Sab in terms of the LKV and its covariant derivatives, 
and imposes the above equation to determine conditions on (and inter-relationships between) 
the various (tensorial) coefficients in this expression. This leads to the final form for Sab in 
terms of a 4 rank tensor with symmetries of Riemann, with certain integrability conditions 
imposed, which eventually identify this tensor with Pabcd- This basically summarizes the 
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key steps in [121] which attempts to generahze the Clausius relation maintaining as close a 
connection as possible with the original result of Jacobson. There are, however, key (and 
inevitable) differences from earlier work along these lines, and several other subtleties which 
arc discussed in detail in [121]. 

Finally, we end this section by clarifying an important point concerning distinction be- 
tween Sec. 4.1.3, where we discussed the equivalence of field equations with the first law of 
thermodynamics, and the discussion in the present section where one attempts to derive the 
field equations from the Clausius relation; that is, we wish to clarify the similarities and 
differences between dE + PdV and SQ, both of which are eventually equated to TdS. This 
was done in [122], where it was pointed out that: 

1. There is a precise sense in which the "heat flux" term of Jacobson, and the "PdV" 
term of Padmanabhan, arc different; this difference is connected with the difference in 
null vectors which govern the horizon deformation in the two cases. 

2. The change in area of a horizon cross-section is determined by the expansion (and not 
its first derivative) of the ingoing null congruence k normalized to have unit Killing 
energy. An explicit expression for the expansion of this congruence can be given in 
terms of combination of curvature tensor components as 

In = RabCk'' - Rabcdu-'n''u'n'' (237) 

dA 

(where A is the afffine parameter along fc, and u and n are unit timelike and spacelike 
vectors spanning the horizon, see [122]), which clearly shows that the area change 
involves not just Ricci, but also the Riemann tensor - a point which is of relevance in 
the context of deriving field equations from thermodynamics. 

3. From explicit (Z) — 2) + 1 + 1 decomposition of the Einstein tensor, one can show that 
there is also a term corresponding to quasi-local energy of the horizon, which must be 
separately accounted for when considering energy flow across the horizon. 

4. The Raychaudhuri equation with the prescribed null congruence has the 0(A) term 
which does not vindicate or necessitate setting the expansion to zero. The most curious 
aspect of the analysis, however, is the fact this term nevertheless yields the null-null 
component of Einstein equations, although the algebraic reason is completely different! 
In fact, in their recent generalization [121], Jacobson et. al. have indeed found that no 
conditions are imposed on their null congruence by the Clausius relation although this 
was invoked in the original paper of [119] and [120]. 

All these facts, along with the previous discussion in this section, shows that while in- 
terpreting the field equations as thermodynamic relations is mathematically rigorous, the 
derivation of field equations from themodynamic considerations has unresolved conceptual 
and algebraic issues when it is attempted in the above manner. We shall see later (see Sec. 5) 
that this is best done from a completely different approach, viz. from a thermodynamic 
variational principle. 
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4.2 Equipartition law for microscopic degrees of freedom 

Given the fact that spacetime appears to be hot, just hke a body of gas, we can apply the 
Boltzmann paradigm ("If you can heat it, it has microstructure" ) and study the nature of 
the microscopic degrees of freedom of the spacetime — exactly the way people studied gas 
dynamics before the atomic structure of matter was understood. One key relation in such 
an approach is the equipartition law AE = {l/2)kBTAN relating the number density AiV 
of microscopic degrees of freedom we need to store an energy AE at temperature T. (This 
number is closely related to the Avogadro number of a gas, which was known even before one 
figured out what it was counting!). 

If gravity is the thermodynamic limit of the underlying statistical mechanics, describing 
the 'atoms of spacetime', we should be able to relate E and T of a given spacetime and deter- 
mine the number density of microscopic degrees of freedom of the spacetime when everything 
is static. Remarkably enough, we can do this directly from the gravitational field equations 
[123, 124, 125]. The establishment of such a relation would obviously be important from the 
perspective of the emergent gravity paradigm as well. 

Since the basic idea here is quite simple (it is explained, for e.g., [123, 124, 125]), it is 
convenient to first describe the result for Einstein theory (which is essentially some well known 
relations interpreted differently) , and then describe the steps needed to derive it more formally 
for the entire class of Lanczos-Lovelock models using the Noether current and potential. We 
shall work with a static spacetime with metric components 500 ~ ~-^^i3oa = 0,9ap ~ ha^- 
An observer at rest in this spacetime with four- velocity u° = Sq /N will have an acceleration 

= (0,a^) where = {d^N/N). We consider a spatial 3-volume V in this spacetime, 
having a boundary dV. 



4.2.1 Equipartition law in Einstein theory 

To fix the ideas, we will begin with the simplest context in which one can determine the density 
of microscopic degrees of freedom, viz., a static spacetime in Einstein's general relativity. In 
a static spacetime, it is easy to show that 

RabU^u" = V,a' = N-^D^,{Na^') (238) 

where is the covariant derivative operator corresponding to the 3-space metric . Using 
the Einstein field equations in Eq. (238), we can relate the divergence of the acceleration to 
the source f^b = [Tab - il/2)gabT) by: 

D^{Na^) = SnNfabu'^u'' = 47rpKomar (239) 

where pKomar is the (so called) Komar mass-energy density. Integrating both sides of Eq. (239) 
over V and using Gauss theorem gives: 

E = / d^xVhpKon..r = 77 / ^^-TT- } (240) 

JV 2 Jgy Lp 1^ 27r J 

where a is the determinant of the induced metric on dV and is the spatial normal to dV . 
(We have temporarily restored G = Lp ^ 1 keeping h—c—ks — l-) 
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This result has a remarkable interpretation directly related to the equipartition law of 
conventional thermodynamics[l23, 124, 125]. To see this, choose dV to he a N = constant 
surface so that the normal is in the direction of the acceleration and a^n^ = \a\ is the 
magnitude of the acceleration. The local, redshiftcd, Davies-Unruh temperature [39] is then 
given by: T = NTioc = {Na''n^/2'K) = {N\a\/2T:). Therefore Eq. (240) is exactly in the form 
of the equipartition law (with fcs = 1) where 

E=U duT; An^y^ (241) 

That is, demanding the validity of Einstein equations, and interpreting Eq. (241) as a law of 
equipartition, we can determine the number of microscopic degrees of freedom in an element 
of area A A = ^/ad'^x on dV to be ^/ad'^x/Lp. If these microscopic degrees of freedom are 
in equilibrium at the temperature T, then the total equipartition energy contributed by all 
these degrees of freedom on dV is equal to the total energy contained in the bulk volume 
enclosed by the surface — which could be thought of as another occurence of the holographic 
principle. 



4.2.2 Equipartition law in Lanczos-Lovelock models 

Remarkably enough, the equipartition law exists for any diffcomorphism invariant theory of 
gravity and allows us to identify the corresponding surface density Anj ^faSx of microscopic 
states in a given theory. However, we restrict our attention to Lanczos-Lovelock models 
since the general case requires an additional assumption involving bifurcate Killing horizon 
which makes it's status a bit unclear. 

To demonstrate how the result generalizes, we follow precisely the same steps [124, 125] 
as in the Einstein case, but now giving a generalized version of the equations encountered 
before. For Lanczos-Lovelock actions, the object which replaces Ricci tensor of Einstein 
theory is defined by Eah + {l/2)Lgab, which, from Eq. (10), gives 

nab = Pa"'"'' Rbcde (242) 

(Note that TZab => (167r)~^i?Q(, in the Einstein case.) We shall also require the expression for 
Noether current J° associated with the timelike Killing vector which is quite simple and 
is given by J" = 27?,^^''. One is now set to rerun the steps outlined for Einstein theory above, 
as follows. 

Using the definition of J°'^ (or expression Eq. (41) directly for a Killing vector), we get: 

2nabu''u'' = VaiJ'^'ubN-^) (243) 

Further, in a static spacctimc, V^Q* = N^^DaiNQ"") for any static vector (with 9oQ' = 0) 
and hence we can write this as: 

2N'Rabu''u^ = D^i^J^^Ub) (244) 

which generalizes of Eq. (238) to an arbitrary theory of gravity. The source for gravity in a 
general theory (analogous to Komar energy density) is defined through 

p EE ANUabU^U^ (245) 
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For Lanczos-Lovelock models, p can be written as a function the matter stress-energy tensor 
Tab (generalizing the Tab used to define Komar energy in Einstein case above) using the 
field equations. This is possible for Lanczos-Lovelock models since trace of the equations 
of motion is proportional to the Lagrangian. However, the explicit form will not be needed 
here. Integrating both sides over a region V bounded by dV, we get 

E = I Vhd^-^x p 

V 

d^-^x^in^UbJ'") - / d'^-^xy/^{Nnc,.r°) (246) 

dV JdV 

which is the analogue of Eq. (240). In Einstein's theory, J^" = q^/Stt which will reduce 
Eq. (246) to Eq. (240). In a general theory, the expression for An is not just proportional to 
the area and quantifies this difference. 

To obtain the equipartition law in a more explicit form, we note that (from Eq. (41)), we 
get J'''' = 2P°^'=''VcCd and hence J"° = 4\a\P°^n'^. So 

E= f d^'-^x^ilQ'KP^^n^n^) { ^) ^ \ I dnT (247) 



27r / 2 



where T ~ N\a\/ 271 is the Davies-Unruh temperature as before. The surface density of 
microscopic degrees of freedom, when the field equations arc satisfied (step 2 above), can be 
now read-off as: 

^^ = 32.P-....- (248) 

where tab = (l/2)(Manb — UbUa) is the binormal to dV. The connection with Wald entropy 
(see section (3.1)) is evident. In fact, 

S={l/2)I3E (249) 

provides the relationship between entropy, energy and temperature in all Lanczos-Lovelock 
models including general relativity as a special case. (In Rcf. [123], the implications of this 
result for general relativity was discussed extensively.) It also follows that, in all these models, 
we have 4 A 5 = An. 

To sum up, the equipartition law allows one to identify the number density of microscopic 
degrees of freedom on a constant redshift surface. Using this one can define the entropy of the 
horizon in a general theory of gravity, which, as shown above, agrees with Wald entropy. These 
ideas will play a key role in our discussion in Sec. 5. We now discuss one specific context 
in which the equipartition idea presented here might have some fundamental implication. 
However, what follows is mathematically independent of the results derived above. 

4.2.3 Holography and quantization of gravitational entropy 

Most of the arguments in the emergent gravity paradigm are free of unnecessary speculations 
since the approach taken is "top-down" and relies on well-established principles. For example. 
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the key result of the previous section expressed in the form of Eq. (248), was derived from a 
classical theory with the only quantum mechanical input being the formula for Davies-Unruh 
temperature. The h enters the expression only through the formula fc^T = {h/c){a/2TT) which, 
along with G in classical gravitational field equations leads to Gh/c^, one of the central 
constants in quantum gravity. So one does not have to make any unwarranted speculations 
from the domain of quantum gravity to obtain Eq. (248). 

However, it is worth exploring where a result like Eq. (248) might emerge from a micro- 
scopic candidate theory of quantum gravity like for example, string theory or loop quantum 
gravity. To begin with, it is clear that our result in Eq. (248) should not depend on the details 
of the correct theory of quantum gravity. Any model of quantum gravity which has (i) the 
correct classical limit and (ii) is consistent with Davies-Unruh temperature, will lead to our 
Eq. (248). Since one would expect any quantum gravity model to satisfy these two conditions 
(i) and (ii) above our "top-down" approach will meet the "bottom-up" approach of such a 
theory in the overlap domain. So if string theory or loop-quantum-gravity or dynamical tri- 
angulations or any other candidate model for quantum gravity satisfies (i) and (ii), it will 
lead to the equipartition. (On the other hand, if a quantum gravity model does not satisfy (i) 
and (ii), such a model is probably wrong in any case.) Even though the thermodynamic limit 
contains far less information about the system than a microscopic description, one could still 
wonder about the possible route or mechanism by which a wide class of candidate models 
in quantum gravity (all satisfying (i) and (ii) above) can lead to a relation like Eq. (248) in 
the appropriate limit. In this section, we will indicate one route through which such results 
can emerge. (Although the motivation for our discussion comes from ideas in the emergent 
gravity paradigm, the results themselves are of independent interest since they generalize 
certain ideas about quantum mechanical spectrum of a black hole and its connection with 
the highly damped quasi-normal modes.) 

Let us first consider Einstein's theory in which the equipartition law assigns A/ degrees 
of freedom to a proper area element A so that A sa nLp for n ^ 1. (This corresponds to 
the semiclassical limit in which one can meaningfully talk about proper area elements etc. in 
the semiclassical metric.) This is, of course, nothing but area quantization in the asymptotic 
limit. Though LQG probably leads to such a result, in the limit we are interested in, it can 
be obtained from the Bohr-Sommerfeld quantization condition applied to horizons. In fact, 
Bekenstein conjectured [126] years back that, in a quantum theory, the black hole area would 
be represented by a quantum operator with a discrete spectrum of eigenvalues. (Bekenstein 
showed that the area of a classical black hole behaves like an adiabatic invariant, and hence, 
according to Ehrenfest's theorem, the corresponding quantum operator must have a discrete 
spectrum.) Extending these ideas to local Rindler horizons treated as the limit of a stretched 
horizon, one can understand how a result like quantization of any spatial area element can 
arise in a microscopic theory. 

The situation becomes more interesting when we go from the Einstein gravity to Lanczos- 
Lovelock models. In Einstein gravity, entropy of the horizon is proportional to its area and 
hence one could equivalently claim that it is the entropy which has an equidistant spectrum. 
But, in the Lanczos-Lovelock models, this proportionality between area and entropy breaks 
down. (Generalization of LQG to Lanczos-Lovelock models does not exist as yet.) The 
question then arises as to whether it is the quantum of area or quantum of entropy (if either) 
which arises in these models. This question was addressed in ref. [127] where it was shown 
that in the Lanczos-Lovelock models, it is indeed the entropy that is quantized with an 
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equidistant spectrum. This matches nicely with the fact that it is the quantity in the right 
hand side of Eq. (248) that takes integral values in the equipartition law. Thus one can 
alternatively interpret our result as entropy quantization. 

We now indicate a sufficiently general 'mechanism' by which any microscopic model of 
quantum gravity can lead to such a quantization condition, in terms of two features. First one 
is the peculiar 'holographic' structure of the action functionals in Lanczos-Lovelock models 

[128, 29, 21]. As noted in Sec , the action functionals in all these theories can be separated 

into a bulk and surface term and the surface term has the structure of an integral over d{pq). 
It can also be shown [3, 128, 29] that the same ^d{pqy structure emerges for the on-shell 
action in the Lanczos-Lovelock models; that is, 

^|o„^shcn« / d'^xd^iU'^'g^k)^ I dH.ir^^gjk (250) 

JV JdV 

where IT^'' is the suitably defined canonical momenta corresponding to gjk- The second 
ingredient is that, in all these theories, the above expression for the action leads to the 
Wald entropy of the horizon. (For example, in Einstein gravity, the surface term in action, 
evaluated on a horizon will give one quarter of area [128, 29].) In the semiclassical limit, Bohr- 
Sommerfeld quantization condition demands that the integral of d{pq) should be equal to 
2?™. Since the action in Eq. (250) has this 'd{pqy structure, the Bohr-Sommerfcld condition 
reduces to AJon-shcU = 27r7i. Since Ajon-shcii is also equal to Wald entropy, we obtain 

•Swald = ^jon-sholl = ZtTTI (251) 

The Bohr-Sommerfeld quantization condition, of course, was the same used originally by 
Bekenstein and others to argue for the area quantization of the horizon but in the more 
general context of Lanczos-Lovelock models it appears as entropy quantization [127]. (It 
is also possible to argue that in the semiclassical limit, the on-shell value of the action 
will be related to the phase of the semiclassical wave function cx cxp (iAjon-shcii)- If 
the semiclassical wave function describing the quantum geometry, relevant for a local Rindler 
observer, is obtained by integrating out the inaccessible degrees of freedom beyond the horizon 
then one can argue that [127], this phase should be 2?™ in the asymptotic limit leading to 
^jon-shoii = 2?™. Given the conceptual ambiguities related to interpretation of 'wave function 
of geometry', it is probably better to invoke Bohr-Sommerfeld condition.) To summarize, 
the three features: (i) the holographic structure of the gravitational action functional (ii) 
the equality of on-shell gravitational action functional and Wald entropy and (iii) the Bohr- 
Sommerfeld quantization condition, combine together to provide a context in which Eq. (248) 
can arise in any microscopic theory. 

While this gives a general result that in Lanczos-Lovelock theories it is the entropy of the 
horizon that is quantized, it would be useful to reinforce it by an explicit calculation in the 
standard context. Fortunately, this can be done for the m ~ 2 Lanczos-Lovelock term, that 
is, for the GB model, using arguments first suggested by Hod [129] based on quasi-normal 
modes of black hole oscillations. We give here a qualitative sketch of the arguments, and 
refer to [127] for technical details. 

We first briefly summarize Hod's argument. He started from Bekenstein's arguments 
regarding quantum area spectrum of a non extremal Kerr-Newman black hole, and showed 
that the spacing of area eigenvalues can be fixed by associating the classical limit of the quasi- 
normal mode frequencies, Wc, with the large n limit of the quantum area spectrum, in the 
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spirit of Bohr's correspondence principle (n being the quantum number). For a Schwarzschild 
black hole of mass M 'm D = A dimensions, the absorption of a quantum of energy ujc, would 
lead to change in the black hole area eigenvalues as, = An+i — An = {dA/dM)ujc and 
for entropy AS" = {dS/dM)ujc- For a Schwarzschild black hole the level spacing of both area 
and entropy eigenvalues were indeed found to be equidistant, allowing one to associate the 
notion of a minimum unit, a quantum, of area and entropy. 

To apply this idea to Lanczos-Lovclock models, we need an expression for the quasi- 
normal mode frequencies associated with black hole solutions in these theories. These are 
not known in general, but are known in the case of 5D Einstein-Gauss-Bonnet (GB) black 
holes. The form of the highly damped quasi-normal modes of these black holes suggest that 
it is the entropy which has a equally spaced spectrum [127]. The result essentially depends 
only on the fact that Ke ujc oi T leading to [dS / dM)ujc (x T{dS/dM) which is a constant. 
On the other hand, since in this case area is a non-linear function of entropy, A = F{S), the 
area spectrum is not equidistant. Qualitative comments concerning the precise value of the 
quantum of entropy (in particular, the condition under which it is 2n) can be found in [127]. 

The case we have made for the quantization of entropy can be analyzed in several different 
ways. If one applies our argument based on quasi-normal frequencies with the imaginary 
rather than real part of the frequencies (a strong case for doing so was suggested by Maggiore 
[130]), then the result for constant spacing of entropy eigenvalues follows simply from the 
fact that this imaginary part, being derivable from a scattering matrix approach [131, 132], is 
independent of the specific model of gravity and depends only on the existence of a horizon. 

The physical implications of the statement: gravitational entropy is quantized, are un- 
clear at present. As should have been evident from our discussion in earlier sections, there 
are various proposed explanations of horizon entropy in the literature, with widely differing 
interpretations. For Lanczos-Lovelock models, most of these reproduce Wald entropy or 
yield expressions closely resembling it (as in Sec. 3.2), and therefore the statement that Wald 
entropy is quantized with equidistant spacing will also generically have implications for the 
physical principles on which a particular approach to horizon entropy is based. Since at least 
some of these approaches remain close to well known physics (for instance, the approach 
described in Sec. 3.2), the statement of quantization of entropy might be of relevance at more 
fundamental level than has yet been appreciated. 

4.3 Emergent aspects of the Lanczos-Lovelock action 

We now turn to the thermodynamic relevance of results discussed in Sec. 2.7 about the 
peculiar structure of the action functional and related issues. We saw that the splitting of 
the total action into two terms led to a holographic relationship between the two. It turns 
out that, there exists another splitting of Lanczos-Lovelock actions which also displays 
similar holographic structure. This separation is a generalization of the split of the Einstcin- 
Hilbcrt Lagrangian as i? = 2{Gq — Rq). For static spacetimc with a Killing horizon, the 
0-th component is uniquely determined by the Killing time and hence unambiguous. In the 
case of Einstein's theory, the first term (Gg) in the above split is related to energy E while 
the second one (i?g) to entropy S, of a (static) horizon. More precisely, the first term is 
associated with the ADM Hamiltonian while the second one to projection of the Noether 
charge corresponding to diffeomorphism invariance of the theory, which gives the entropy of 
horizon in such theories. This gives the total action a thermodynamic interpretation as the 
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free energy (at least in spacetimes with horizon). Further, it can be shown that the new pair 
above also obeys the same holographic relation Eq. (77); the surface term in the pair arises from 
the Rq part of the split - a result which holds for static spacetimes [21]. Since the bulk and 
surface terms in this case are related to energy and entropy, the holographic relation between 
them may be thought of as analogous to inverting the expression for entropy 5* = S{E, . . .) 
to obtain the energy E = E(S, . . .) in terms of entropy in normal thermodynamic systems. 

It can be shown that the very same holographic relationship between the bulk and the 
surface term exists for an arbitrary splitting of the Lagrangian into Lgur and Lbuik, provided 
the trace of the equations of motion is proportional to the Lagrangian: g^-^Eat oc L. This is 
certainly true for Lanczos-Lovelock Lagrangians, although there is no rigorous proof as to 
whether there are any other Lagrangians satisfying it. (The proofs for the above statements 
can be found in [21].) We will now discuss these features. 



4.3.1 Lanczos-Lovelock action as free energy of spacetime 

We now present one more aspect of the gravitational action functional which adds strength 
to the thermodynamic interpretation. One can show that, in any static spacetime with a 
bifurcation horizon (so that the metric is periodic with a period f3 = 2tt/k in the Euclidean 
sector), the action functional for gravity can be interpreted as the free energy of spacetime for 
all Lanczos-Lovelock models. The Nocthcr current for Lanczos-Lovelock theory 

ja ^ (2£:°''eh + LC) ■ (252) 

will play a prominent role in our discussion, where we have assumed staticity (so that Sv"" = 0), 
and will work in the natural coordinate system in which the timelike Killing vector field 
= (1,0) explicitly exhibits the static nature. 

We start with the Noether current J° and integrate it on-shell (i.e., with 2Eab = Tab), 
over a constant-i hypersurface with the measure dSa = S';lNVhd^^^x where = and h 
is the determinant of the spatial metric. (We will work in the Euclidean sector.) Multiplying 
by the period /? of the imaginary time, we get 

Pjrdj:a = 13 j T^^^dl^a+P j LCdl^a 

= J{(3N)T^Uau''Vhd^-^x + j\tE J Ly/^d^-^x (253) 



where w° = (^'^/N = N^^Sq is the 4- velocity of observers moving along the orbits of 
and dSfi = Ua\/hd^^^x. The term involving the Lagrangian gives the Euclidean action for 
the theory, while in the term involving Tab, we note that /37V = /3ioc is the rcdshifted local 
temperature. If we define the (thermally averaged) energy E as 

j {J3N)T;^Uau''Vhd'^-^x = J PiocTl^UaU^'Vhd^-^x = f3E (254) 

we obtain 

A = I3 j J^d^a-I^E (255) 
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Upon using Stokes theorem, we see that the first term (with the Noether charge) gives horizon 
entropy. Therefore, we find that 



A = S-f3E = -I3F 



(256) 



where F is the free energy. (Usually, one defines the Euclidean action with an extra minus 
sign a.s A = —Ae so that the Euclidean action can be interpreted directly as the free energy.) 
One can also obtain from Eq. (253) the relation: 



S = I3 J''dT.a = 



-gd°x{p + L) 



(257) 



where p = Tabu^u". This gives the entropy in terms of matter energy density and the 
Euclidean action. Alternatively, if we assume that the Euclidean action can be interpreted 
as the free energy, then these relations provide an alternative argument for interpreting the 
Noether charge as the entropy. (Similar results have been obtained in Ref. [42, 43] by a more 
complicated procedure.) 

There is another way of expressing the above result by introducing a vector Ik = Vkt = 
dkt = (1, 0) which is (un-normalized) normal to t = constant liypersurfaces. The unit normal 



is h 



Nl, 



~Uk in the Lorcntzian sector. Using J" ~ \/i,J , wc can write 



2 E^^'' 



(258) 



and since h^aJ = Va(/fc J ), wc get 



'9 



-da 



(259) 



which is equivalent to the previous expression with J", but now expressed in a more formal 
way and using Noether potential (rather than current). 

The free-energy structure of the Lanczos-Lovelock action can be very explicitly demon- 
strated for black hole solutions which are spherically symmetric, in which case, the the action 
can be written as a total derivative. The details can be found in section 4 of Ref. [60] , and 
follows from an algebraic identity which holds for Lanczos-Lovelock Lagrangians of a given 
order m evaluated for spherically symmetric spacctimes with —.goo = f{^) = 9^^ (with horizon 
at /(r = a)= 0): 



Lr. 



d 

dr 



{D-2m- 1)! 



(1 - /)"V 



m^D — 27n — l 



d 

dr 



(D - 2m) 



-i^D- 



■2rn 



(260) 



The Euclidean action, given by dr J d'^x,/gELm[f(r)], when integrated over a spacetime 
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region bounded by the horizon at r = a and an arbitrary radius r = b{> a), becomes 



-A, 



167r 



16^ 



(D - 2m - i; 



{D-2)\ 



{D-2m)\ 



m/'(l - 



m-l„D-2m 



{D-2)\ 



,_D-2m-l 



(D- 2m- 1)! 



m(£-2)! „ 
^ ' {D-2m)\ 



(167r/Sl)£;„ 

-5,„ + Q[/(6),/'(6)] 



4S,„/n 



Q[f{h),nh)] 



where VL is the area of a unit (D — 2)-sphere, and we have used known expressions for 
horizon energy and entropy for sphericahy symmetric black holes in Lanczos-Lovelock theory. 
(The functional Q depends on behaviour of the metric near r — b, and is irrelevant for our 
argument.) The free energy structure is self-evident in the contribution to the Euclidean 
action on the horizon, thereby supporting the general arguments given earlier in this section. 



(261) 



5 Entropy density of spacetime and Lanczos-Lovelock 
models 

As already described in the previous section, one of the most important basis for treating 
gravity as an emergent phenomenon is the possibility of deducing field equations of gravity 
from conventional laws of thermodynamics applied to spacetime horizons perceived by rele- 
vant class of observers. Since horizons play a central role in such discussions, it is natural 
to consider an approach in which null surfaces provide the fundamental variables from which 
gravitational dynamics of the background spacetime emerges. Indeed, since null surfaces act 
as one-way membranes and can block information for a specific family of observers (time-like 
curves), it is natural to attribute an entropy to such surfaces. 

If we take this point of view seriously, then the deformations of spacetime (x* — x*) = 
associated with a vector field becomes analogous to deformations of a solid in the study 
of elasticity. By and large, such a spacetime deformation is not of much consequence except 
for the null surfaces. As we have described earlier, null surfaces can be thought of as acting 
as a local Rindler horizon to a suitable set of observers. The deformation of a local patch of 
a null surface will change the amount of information accessible to the local Rindler observer. 
Equivalently, such an observer will associate certain amount of entropy density with the 
deformation of a null patch with normal n". This suggests that extremizing the sum of 
gravitational and matter entropy associated with all null vector fields simultaneously, could 
lead to the equations obeyed by the background metric. 

Conceptually, this idea is similar to the manner in which we determine the influence of 
gravity on other matter fields. If we fill the spacetime with freely falling observers and insist 
that normal laws of special relativity should hold for all inertial observers simultaneously, we 
can arrive at the generally covariant versions of equation obeyed by matter in an arbitrary 
metric. This, in turn, allows us to determine the influence of gravity on matter fields and 
fixes the kinematics of gravity. To determine the dynamics, we do the same but now by fllling 
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the spacetime with local Rindler observers. Insisting that the local thermodynamics should 
lead to extremum of an entropy functional associated with every null vector in the spacetime, 
we should be able to obtain the equations which determine the background spacetime. 

There is no a priori assurance that such a program will succeed and hence it is a surprise 
that one can actually achieve this. Let us associate with every null vector field n°(x) in 
the spacetime a thermodynamic potential (say, entropy) which is quadratic in n° and 

given by: 

'^K] = ^grav K] + 3,„a« ["'^1 = - {iP^^V aU'^V - Tabu'^n'') , (262) 

where P^J^ and Tab are two tensors that play a role analogous to elastic constants in the theory 
of elastic deformations. If we extremize this expression with respect to n°-, one will get a 
differential equation for n'^ involving its second derivatives. We, however, want to demand 
that the extremum holds for all thereby constraining the background geometry. Further, 
the insistence on strictly local description of null-surface thermodynamics translates into the 
demand that the Euler derivative of the functional S5(n°) should not contain any derivatives 
of n'\ 

It is indeed possible to satisfy these conditions by the following choice: We take P°^ to be 
a tensor having the symmetries of curvature tensor and divergence-free in all its indices and 
take Tab to be a divergence-free symmetric tensor. (The conditions ^ aPcd ~ 0' ^ aT§ = 
can be thought of as a generalization of the notion of "constancy" of elastic constants of 
spacetime [133].) Once we get the field equations we can read off Tab as the matter energy- 
momentum tensor; our notation anticipates this result. We also know that the P"-^'^'^ with 
the assigned properties can be expressed as P^^ = dL/dR^^ where L is the Lanczos-Lovelock 
Lagrangian and Rabcd is the curvature tensor [3, 4]. (Again, the notation already anticipates 
the fact that it will turn out to be the same object as we have introduced in earlier parts 
of this review. However, since we are now hoping to arrive at the Lanczos-Lovelock actions 
in the context of emergent gravity paradigm, we hope to present a new perspective on the 
physical relevance of this tensor.) This choice in Eq. (262) will also ensure that the equations 
resulting from the entropy extremisation do not contain derivatives of the metric of higher 
order than second. 

We now demand that d'^/Sn'^ ~ for the variation of all null vectors n" with the condition 
rian" = imposed by a Lagrange multiplier X{x)gabn°'n^ added to S[n'']. An elementary 
calculation and use of generalized Bianchi identity and the condition V aT^ = leads us to 
[3, 4, 134, 135] the following construct on background geometry: 

Qb = n- \fbL = ]^Tl + K51 (263) 

where A is an integration constant. These are precisely the field equations for gravity in a 
theory with Lanczos-Lovelock Lagrangian L with an undetermined cosmological constant A 
arising as an integration constant. 

The thermodynamical potential corresponding to the density 3 can be obtained by in- 
tegrating the density 3[n°] over a region of space or a surface etc. based on the context. 
The matter part of the 5, proportional to Tabn°'Ti^, will pick out the contribution (p -I- p) 
for an ideal fluid, which is the enthalpy density. Multiplied by /3 = 1/T, this reduces to 
the entropy density because of Gibbs-Duhem relation. When the multiplication by /3 can 
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be reinterpreted in terms of integration over (0,^) of the time coordinate (in the Euchdean 
version of the local Rindler frame), the 3 can be interpreted as entropy and the integral over 
space coordinates alone can be interpreted as rate of generation of entropy. (This was the 
interpretation provided in Refs. [3, 4, 134, 135] but the final result is independent of this 
interpretation as long as suitable boundary conditions can be imposed). One can also think 
of ^[n°'] as an effective Lagrangian for the collective variables n° describing the deformations 
of null surfaces. 

The gravitational entropy density in Eq. (262) can be expressed in terms of the Killing 
{£q9ab ~ '^aQb + b^a) and isocutropic {Jab — '^a^b — bQa) parts of the deformations 
introduced in Ref. [5] as: 

-AP^^S/aq^ybq'' = P'''''{£qg^,){£ggbd) - {l/2)P-'-^JabJcd (264) 

So the gravitational entropy density has two parts: one contributed by the square of the 
Noether potential (which vanish for isocntropic deformations with qa = Va4>) and another 
which depends on the change in the metric under the deformation (which will vanish for the 
Killing deformations). When qj is a pure gradient, Jab will vanish and one can identify the 
first term in Eq. (264) with a structure like Tr{K'^) — (TrK)'^. On the other hand, when qa 
is a local Killing vector, the contribution from Sij to the entropy density vanishes and the 
entropy density is just the square of the Noether potential Jab- For a general null vector, 
both the terms contribute to the entropy density. Variation of entropy density with respect to 
either of the two contributions (after adding suitable Lagrange multiplier to ensure vanishing 
of the other term) will also lead to the gravitational field equations. 

In this approach, there arise several new features which are worth mentioning. 

First, the extremum value of the thermodynamic potential, when computed on-shcll for a 
solution with static horizon, leads to the Wald entropy[134] . This is a non-trivial consistency 
check on the approach because it was not designed to reproduce this result. When the field 
equations hold, the total entropy of a region V resides on its boundary dV which is yet another 
illustration of the holographic nature of gravity. 

Second, the entropy functional in Eq. (262) is invariant under the shift Tab Tab + 
PoQab which shifts the zero of the energy density. This symmetry allows any low energy 
cosmological constant, appearing as a parameter in the variational principle, to be removed 
thereby alleviating the cosmological constant problem to a great extent [135, 136]. As far as 
we know, this is the only way in which gravity can be made immune to the zero point level of 
energy density. It is again interesting that our approach leads to this result even though it is 
not designed for this purpose. This works because the cosmological constant, treated as an 
ideal fluid, has zero entropy because p + p = and thus cannot affect the dynamics in this 
perspective, in which gravity responds to the entropy density rather than energy density. 

Third, the algebraic reason for the idea to work is the identity: 

4P„tVc""Vrfn'' = 27^afc?^''n'' + V j4P„f n-^V^n^] (265) 

which shows that, except for a boundary term, we are extrcmising the integral of {2TZab — 
Tab)n°'n'' with respect n" subject to the constraint Uau'^ = 0. The algebra is trivial but not 
the underlying physical concept. 

Fourth, the gravitational entropy density — which is the term in the integrand '^grav oc 
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{-Pab'^cn^Vdn") in Eq. (262) — also obeys the relation: 



^^''^^ -P:^V,n'')^^{Wan^-S:W,n^) (266) 



where the second relation is for Einstein's theory. This term is analogous to = — (5^7^ 
(where Kab is the extrinsic curvature) that arises in the (1+3) separation of Einstein's equa- 
tions. (More precisely, the projection to 3-space leads to t^.) This term has the interpretation 
as the momentum conjugate to the spatial metric in (1+3) decomposition and Eq. (266) shows 
that the entropy density leads to a similar structure. That is, the canonical momentum con- 
jugate to metric in the conventional approach and the momentum conjugate to in 3giav 
are similar. 

One final question which must be asked is what justifies the starting expression for entropy 
functional as the "right" one. This has no answer in absence of a complete quantum gravity 
model. The next best thing one can do, and as we have done repeatedly in this review, is 
to compare the on-shell value of the expression near a Rindler horizon with Wald entropy. 
This can be easily done, and the entropy functional described here can be shown to reproduce 
correctly the known expression for Wald entropy on-shell. The proof of this is straightforward 
and can be found in [134]. 



6 Conclusions 

The Lanczos-Lovelock models provide a simple and elegant generalization of Einstein's the- 
ory of gravity to higher dimensions. If it turns out that we actually live in a spacetime 
with more than 4-dimcnsions, then Lanczos-Lovelock models probably will prove to be of 
practical utility in the study of nature. But even if this is not the case, Lanczos-Lovelock 
models are of great interest in providing a natural backdrop to test different conceptual and 
mathematical aspects of Einstein's theory. This has been the point of view taken up and 
illustrated throughout the review. On several occasions we have pointed out how looking 
at a generalized result in Lanczos-Lovelock theory allows us to obtain a better insight into 
the phenomena in Einstein's theory. Such a study also allows us to delineate special features 
and numerical coincidences which occur in the Z) = 4 Einstein gravity from results of much 
greater generality. 

One of the surprising developments in recent years is the effortless generalizations of many 
concepts related to thermodynamics of horizons from Einstein gravity to these models. Based 
on this one could also argue that the entire emergent paradigm of gravity possesses a much 
broader basis than one would have first guessed from the study of Einstein's theory. Because 
of its importance, we have used these topics as the focal point for this review. We hope it 
helps to attract more attention from the community to this fascinating subject. 
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Appendices 



A Algebraic aspects of certain tensors in higher curva- 
ture theories of gravity 

We here give a detailed analysis of certain tensors defined in Sec. 2.1 which appear naturally 
while varying a general Lagrangian which is a function of (and only of) metric and Ricmann 
tensor and clarify many of the issues listed in Sec. 2.1. We will prove some useful identities 
which are non-trivial in the sense that they do not emerge from the algebraic symmetry 
properties of tensors such as Pabcd defined in Sec. 2.1; instead they arise due to the diffeo- 
morphism covariance of the Lagrangian. This leads to the fact that, when scalar quantities 
are constructed from tensorial objects (like in the case of L built from Rabcd and g*-'), the 
derivatives of the scalar quantity with respect to these tensors must obey some identities. 
While a few of these results can be proved by other methods in special cases, the procedure 
outlined here, based on reference [7], is very powerful and readily adaptable to more general 
cases. 

While varying L — L[gab, Rabcd] to obtain the equations of motion, we come across two 
tensors in the following partial derivatives 

jabcd _ I \ . uab I dL 



KdRabcdJg^^ KdgabJn^^^^ 

The subscripts specify the variables being held fixed, which, of course is, what gives meaning 
to the partial derivatives. The P"-'"^'^ is what we called the entropy tensor of the theory 
[41, 3]. The derivatives appearing in Eq. (A.l) are defined within the subspace of infinitesimal 
deformations which preserve the algebraic symmetries of the independent variable. Therefore, 
P"'' is symmetric while P'^^'^'^ inherits the algebraic symmetries of the curvature tensor: 

pabcd pbacd pabdc. pabcd pcdab. pa[bcd] q 2^ 

One can construct several such tensors — which have the properties in Eq. (A. 2) — from the 
curvature tensor and the metric. But these properties alone do not guarantee that such a 
tensor can be expressed as the derivative of some scalar with respect to the curvature tensor. 
The fact that P"'"='^ has the form in Eq. (A.l) leads to several further properties and inter- 
relationships which form the main thrust of this appendix. In particular, we will prove the 
following results: 

(1) Let us construct the tensor 

^ab ^ p^^Vkpb^^^ (^ 3) 

which can be thought of generalized Ricci tensor in these theories. We will prove that 
P"'' ~ —272,°''. That is. the derivatives of an arbitrary scalar function with respect to the 
metric and the curvature are not independent but are connected by the relation: 
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It follows from Eq. (A. 4) that the tensor TZ.""^ is symmetric. (This was noted earlier in a 
different context in [137]). This result does not follow merely from the algebraic symmetries 
of P"-'"^'^ in Eq. (A. 2) above. The nature of the proof which we give below shows that it holds 
only because P"-^'^^. expressible as a derivative of a scalar with respect to the curvature 
tensor. 

(2) From Eq. (A. 4) we can obtain several other relations connecting the partial derivatives 
when different independent variables like {g°'^ , Rabcd) or {g"''' , R°'i,cd) or {g°'^,Rtd) ^^ed to 
describe the system. First, if we use the pair {g""^ ^ Rabcd) as independent variables, we get: 

gaigbm = 27Cim (A. 5) 



Instead, if we use the pair {g"''', R'^bcd) as independent variables, one can also show 

bed 

and, more interestingly, if we use the pair (5°^,-RI?^) as independent variables, we obtain 

(3) Using the fact that p^'^'^d- is anti-symmetric in b and c, one can show 

VfcVcP''""' = i {n"^ - 7^*) = (A.8) 

where the last equality follows from the symmetry of TZ"^^. Thus we get another surprising 
result that, for any scalar L built from curvature tensor and the metric, we have the identity 

V,V, (j^) = (A.9) 

KdRabcdJg^^ 

We shall now sketch a proof of these results. 



Proof of the identities: 



As wc mentioned in section (2.1), all the identities we want to prove follow from diffcomor- 
phism cinvariancc of the Lagrangian. So let us begin with an infinitesimal diffeomorphism 
— > + ^"(x) which changes L, gab and Rijki by infinitesimal amounts. The essential 
trick is to express the Lie derivative of L in two different ways and equate the results. First, 
because L is a scalar which depends on a:* only through gab{x) and Rijki{x), we have: 

£^L = f "V™L - rP'^^'^mgab + rP'^^^'^'^mRabcd = P'^^'^T^ mRabcd (A.IO) 

because Vmffafc = 0. On the other hand, if we think of the change 5L in L, due to small 
changes in 5gab = £^gab and SRijki = £^Rijkh we also have: 

£^L = P'^'£^gab + ?''"'£ eR,jki (A.ll) 
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We will now expand the right hand side of Eq. (A. 11) bringing it into a form similar to the 
right hand side of Eq. (A. 10) and equate the two expressions. That will lead to Eq. (A. 4). 
The first term on the RHS of Eq. (A. 11) is easily dealt with; using £i^gab = Va^fe + Vb^a and 
the symmetry of P°'^ wc obtain 

P'^'^l^Qab = 2P^''Wa^b (A.12) 

The term involving Lie derivative of curvature tensor gives (on contracting with py*^' and 
using Eq. (A. 2)) 

P'^''£^R,jkl - P'"'T^.nR^Jkl + 2P''^' [{'^^DRmJkl + {V kDRrjml] (A.13) 

With some further index manipulations, this can be written as 

P'^'*^' £^R,jki = P''''T^mR^Jkl + 4(V,^™)7e™ (A.14) 

where we have used the definition 7^"" = P^^'^'i?"^^;. Using Eq. (A.14) and Eq. (A.12) in 
Eq. (A. 11) we get 

£^L = P'^'^'CV„%fc/ + 2(V,U)[^'"" + 27^""] 

= £5L + 2(V,em)[P"" + 27^™■] (A.15) 

Since is arbitrary, it follows that the term within square brackets in the last expression 
has to vanish, leading to the condition in Eq. (A. 4): 

p^n^ ^ f ^ -2P'^''' R"\^,i = (A.16) 

Since the left-hand-side is symmetric, it follows that 72.™' must be symmetric. This symmetry 
is certainly not obvious from the definition of TZ"^ itself. 

Equation (A. 8) can be be obtained immediately from the above result, as follows: 

2VbV^P''''^ = [Vb,V^]P'"^'^ 

rib Tjkcid , tdc jjbkid , jji jjbckd , rid jjbcik 

— ^ kbc^ + ^ kbc^ + kbc^ + ^ kbc^ 

= -n'^'+TZ"^ = (A.17) 
where the anti-symmetry and pair-exchange symmetries of P^"^"^ have been used. 

Important Corollaries: 

One can derive a few corollaries from Eq. (A.16) which are useful in expressing the field 
equations obtained from our Lagrangian. Wc first note that, since Sgab = ^daigbjSg''-' ^ we 
have: 

Let us next consider how our results change if we use the pair (g*™, R'\cd) ^^e pair (5"'', R^^j) 
as the independent variables. Obviously that {dL/dg^"^) with Rabcd held fixed is not the same 
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as {dL/dg'^"^) with R'^^,^^ held fixed. To find their relationship, we note that 

' _d_ 



dL = (^) dg^^ + P^'"'dRatcd^ (^) dg"^ + P,''^'d{g^''Rkbcd) 



P^^'^^^dRmbcd + rfff"" <! {-^ ] + 7^„„ !> (A.19) 



bed 



from which it follows that 

dL \ f dL 



- n. 



= 2TZi,n — TZim = T^irn — Pi Rmbcd (A. 20) 

More interestingly, if we work with (g"'', R]jli) as the independent variables we get in place of 
Eq. (A.19) the reuh: 



^ ( ) dg"'' + P dRabcd = P'""'''dRabcd + dg 

^ ' Rabad I 

leading to 

dL \ f dL 




(A.22) 



This result shows that scalars built from g'^ , i?"^ are actually independent of the metric tensor 
when is held fixed. 

We conclude by proving another important corollary: The tensor VmV„P"'""'' is sym- 
metric in indices a and b. This result is important to establish directly the symmetry of 
Eab in Eq. (10). We have already proved that TZat is symmetric, which means the last term 
VmV„P°'""^ in Eq. (10) must also be symmetric in a and b. To show this explicitly, consider 
the anti-symmetric part of this tensor 

V,„V„ (F''"'"'' _ pbmna^ ^ y^^y^ (^pamnb ^ panbrr.^^ ^ _^ ^^^^ ^^pabvan 23) 

To obtain the first equality we have used the pair exchange symmetry and the anti-symmetry 
of the first two indices p^*^'; in arriving at the second equality we have used the cyclic relation 
pi\jkt\ _ Q follows that the last term in Eq. (10) is symmetric only if the right hand side 
of Eq. (A. 23) vanishes, which is assured by Eq. (A. 17) and the symmetry of P"-^'^'^ under pair 
exchange. 

Therefore, we see that each of the terms in Eq. (10) are individually symmetric. The 
results of this section are based on arguments which would apply to Lagrangians more gen- 
eral than the ones we have considered (and hence more general than the Lanczos-Lovelock 
models) , presumably giving more identities etc when, say, covariant derivatives of Riemann 
are also involved. In some sense, the main basis for deriving such identities lies in the very 
first equation Eq. (A. 10) which actually encapsulates the scalar nature of the Lagrangian L. 

A more detailed discussion of the results proved here, and some further implications of 
the same, can be found in [71. 
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B Construction of the Gauss-Bonnet lagrangian 



Let us consider building a Lagrangian L2 which is a hncar combination of quadratic terms 
in curvature , RatR"''' and RabcdR"^'"^'^ , The general form for L2 can be written as 



Lo = 



aR^ + PRabR"^ + iRabcdR" 



bed 



and its variation will have the form 

6 {V^ L2) = Ef Sgab + ^kSv'' 
A long but straightforward calculation gives 
1 



- ab T 

-g L2 



2aRR 



ab 



(2/3 + 47) RkcR 



kacb 



and 



2M 



cbak 



+ 



M-^'^'yjgab 
(^AaRg'^^'^g^'^^ - 



(B.l) 

(B.2) 

(B.3) 
(B.4) 

(B.5) 

(3 {2g^^R'"'' - g^-'R^" - g^^R^") + 47i?''"'^") + {a ^ blB.&) 



AjR^^^R', + 27i?L,i? 



{2a + /3 + 27) VVR - ( 2a + I ) g'^'OR - (^ + 47) DR'"' 



For Eab not to contain higher than second derivatives of the metric, one must demand 

2a + /3 + 27==0 (B.7) 
4a + /3 = (B.8) 
/3 + 47 = (B.9) 

The solution is a = C2 = 7, {3 — — 4c2, where C2 is some arbitrary constant. We therefore 
arrive at the following final form for the full second order Lagrangian 



L = (167rG)"^ [R+ {32Tr^y'^aL2] + boundary terms 



(B.IO) 



where [a] = (length)^ and we have set C2 — (327r^) ^ (this convention follows from the fact 
that Gauss-Bonnet term will then be the Euler character in AD) with 



L2 (X [R"^''^'^ Rabcd — AR'^^Rab + 



(B.ll) 



The Einstein-Gauss-Bonnet (EGB) system, with lagrangian given by a linear combination 
of Einstein-Hilbert lagrangian and L2 above is the simplest truncation of Lovelock model 
which has been well studied in the literature. For example, some consequences of the quasi- 
linearity of Einstein-Gauss-Bonnet (EGB) equations on the Cauchy problem as well as on 
"junction conditions" are discussed in [138], while [139] focusses on discussion of the conserved 
charges in EGB. 
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C Proof of identities in Eqs. (146) - (148) 



To prove the required identities, we will use the following two relations concerning Lanczos- 
Lovelock actions of order m in D dimensions (see Eqs. (30)) 



E 



E 



(m) - 



^ ciaibi . . .a„ 

2 167r2™ 3<^idi...c„ 
D-2m n 



1 aioi 



(C.l) 



where El^,^\ is the equation of motion tensor for Lanczos-Lovelock action; for e.g., for m = 1, 



E 



J(l) 



j{m) 

(167r)~^G' where G* is the Einstein tensor. Note that for maximally symmetric 



spacetime, which is the case we are considering for the (D — 2) sub-manifold corresponding 
to the {D — 2) dimensional horizon, we have G2 = G3 = • • • = G^l} due to isotropy. Hence 
we can write the second relation in Eq. (C.l) for equation of motion tensor for the {D — 2) 
dimensional maximally symmetric subspace as 



E 



A{D-2} 
B{rn) 



{D -2)- 2m 
D -2 



L 



D-2 



(C.2) 



Using the above identity it is easy to prove Eq. (146) as 
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(C.3) 



where to obtain the second equality we have used Eq. (117) and replaced Rq% with the 

intrinsic curvature R-cd defined completely in terms of the induced metric ^ab of the horizon. 
Such a replacement is valid since the corresponding extrinsic curvature vanishes due to our 
assumption of staticity. The last equality then follows from using Eq. (C.2). To prove the 
second identity of Eq. (147) note that 
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(C.4) 



where to obtain the second equality we have used = which is true due to (i) static 
nature of the background geometry and (ii) maximal symmetry of the (D — 2) dimensional 

subspace and we have again replaced R^f) with the intrinsic curvature Rq%. Further, in the 
last equality we have again exploited maximal symmetry to get a proportionality with a 
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proportionality constant Aim which we can determine by contracting the above equation 
with R^^ to get 

Lif^^'^ - {Q^E)m~"^ k% = M.i^-'^S^Sk'B-iMi^-'^^-'R (C.5) 

This completes the proof of identity in Eq. (147). To prove the third identity of Eq. (148) we 
again use the same procedure and write 

r[SODFl/i •••1^2^-4] 0)11)12 5A'2m-5M2m-4 _ c-[B D F Bi- ■ ■ - i] nAiA2 fyMm^sMm^i 

"[A0CEfll---fl2m-i] •^11^2 ' " ^!'2m-5f2m-4 ~ " [AC E Al ■ ■ ■ A2,„ - i] B i B2 " ' ^ B2„^ ^ & B2^ - i 

= nL^-'^S^SE (C.6) 

with the proportionality constant N^_i'^ determined by contracting both sides of the above 
expression by and then taking a trace to get 

(16.)2™-^^^Lf„I? = ^(f„f ^(16.)2ig-iLf- (C.7) 

This completes the proof of identity in Eq. (148). 
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